Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,567,248 G→T 20.9% R331S (CGC→AGC)  asd ← aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,567,2480GT20.9% 130.9 / 34.0 74R331S (CGC→AGC) asdaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base T (12/3);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 2.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

GCGGCATCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAAT  >  CP009273/3567102‑3567396
                                                                                                                                                  |                                                                                                                                                    
gcggcATCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATGTTCAGCTTACGCCggcggc                                                                                                                                                    >  1:378065/1‑149 (MQ=255)
     aTCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTAcc                                                                                                                                               >  2:436026/1‑149 (MQ=255)
     aTCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACCCTGGTCCCATATTCCGCTTACGACGGCCGCCTAcc                                                                                                                                               >  2:500658/1‑149 (MQ=255)
      tCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCg                                                                                                                                              <  2:1034507/149‑1 (MQ=255)
      tCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCg                                                                                                                                              >  1:881522/1‑149 (MQ=255)
             ccAGCTTTATTGGATAACGAGTACCGCAGTTGACGAAGCACCCGTCGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCTACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGt                                                                                                                                       <  1:84767/149‑1 (MQ=255)
              cAGATTTAATGAATACAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTGAAGGCTGACAGGAACTCTGGTCCCCTATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTc                                                                                                                                      >  2:1116510/1‑149 (MQ=255)
              cAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCCGCTTACGCAGGCGGCCTACCGGCGTGGTc                                                                                                                                      >  2:292893/1‑149 (MQ=255)
              cAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGGGGTc                                                                                                                                      >  1:334618/1‑149 (MQ=255)
              cAGATTTAATGAATAAAGACTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTc                                                                                                                                      >  1:437297/1‑149 (MQ=255)
                    tAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTg                                                                                                                                >  2:13867/1‑149 (MQ=255)
                      aTGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGcc                                                                                                                              <  1:705298/149‑1 (MQ=255)
                              aGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACgg                                                                                                                      <  2:693813/149‑1 (MQ=255)
                                  tACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGc                                                                                                                  <  1:391477/149‑1 (MQ=255)
                                    cGCCAGTTGACGAAGCATCCGACGCAGCGGCTACGCGGCCCCACACAGCAGCTGGTCGCCCCCGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTAAGCGTGCCGGTAACGGCGGCTg                                                                                                                >  1:471479/1‑149 (MQ=255)
                                     gCCAGTGGAAGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTgg                                                                                                               <  1:87776/149‑1 (MQ=255)
                                      ccAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAgagggg                                                                                                              >  1:13907/1‑144 (MQ=255)
                                      ccAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACCGTAAAGGCTGACCGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTCCCGGCGTGGTCAGCGTGCCGGTACCGGCCGCTggg                                                                                                              >  1:1065856/1‑149 (MQ=255)
                                          ttGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGCCAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCGGGGGTTa                                                                                                          >  2:650567/1‑149 (MQ=255)
                                          ttGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTa                                                                                                          >  2:247624/1‑149 (MQ=255)
                                             acgacgCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCt                                                                                                       <  2:690898/149‑1 (MQ=255)
                                             aCGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCt                                                                                                       <  2:713978/149‑1 (MQ=255)
                                                aaGCATCCGACGCAGCGGCTCCGCTGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCAc                                                                                                    >  1:920447/1‑149 (MQ=255)
                                                aaGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGCCGCGGTCAGCGTGCCGGTAACGGCCGCCGGGGTTAGCCCAc                                                                                                    >  2:60182/1‑149 (MQ=255)
                                                 aGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACg                                                                                                   <  1:436026/149‑1 (MQ=255)
                                                 aGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACg                                                                                                   <  1:500658/149‑1 (MQ=255)
                                                  gCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGTCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGACGGGGTTAGCTCACGc                                                                                                  >  1:875248/1‑149 (MQ=255)
                                                  gCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTCCCCGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGccaccgc                                                                                                  >  1:875303/1‑143 (MQ=255)
                                                  gCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTACCGGCCGCTCGGGTTCTCTCAccg                                                                                                  >  1:459008/1‑147 (MQ=255)
                                                                 acTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGa                                                                                   >  2:908958/2‑149 (MQ=255)
                                                                       cgGCCCCCGCCAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTc                                                                             <  2:936232/149‑1 (MQ=255)
                                                                       cgGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGGCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGCAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTc                                                                             <  2:736941/149‑1 (MQ=255)
                                                                          ccccccACAGCAGCTGGTCGCCCCCGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCTGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGa                                                                          <  1:311099/149‑1 (MQ=255)
                                                                          ccccccACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGa                                                                          <  2:288213/149‑1 (MQ=255)
                                                                          ccccccACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGa                                                                          <  1:992484/149‑1 (MQ=255)
                                                                           cccccACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGACCTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGTGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGaa                                                                         <  1:748979/149‑1 (MQ=255)
                                                                           cccccACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACCCATAGTGACTTCCCGGTCGTTCGGag                                                                         >  2:122316/1‑148 (MQ=255)
                                                                              ccACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGa                                                                      <  1:892289/149‑1 (MQ=255)
                                                                                    cagcCGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCg                                                                <  1:186354/149‑1 (MQ=255)
                                                                                         ggTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGACCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCAc                                                           >  2:1209470/1‑149 (MQ=255)
                                                                                            cGCCCACGGTAAAGGCTGACAGAAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCCCGCATAGTGATTTCCCGATCGTTCGGAACGACCTTCGCCCACggg                                                        >  1:943564/1‑148 (MQ=255)
                                                                                               ccACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATGGCGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATtg                                                     >  1:955857/1‑149 (MQ=255)
                                                                                                  cggacaaGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTgc                                                  <  2:1106109/144‑1 (MQ=255)
                                                                                                  cGGTAAAGGCTGACAGGACCTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTgc                                                  <  1:122443/149‑1 (MQ=255)
                                                                                                     tAAAGGCTGCCAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGca                                               <  2:967783/149‑1 (MQ=255)
                                                                                                     tAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGca                                               >  1:554938/1‑149 (MQ=255)
                                                                                                     tAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGca                                               >  1:1169886/1‑149 (MQ=255)
                                                                                                     tAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCACACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGca                                               >  1:1211692/1‑149 (MQ=255)
                                                                                                     gaacggCGGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCGCGACCTCTCAGAACGACTTTCGCCCACGGATTGTGCGca                                               <  2:881522/145‑1 (MQ=255)
                                                                                                       aaGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGGTGGGGTTAGCGCACGCATAGTGATTTCCCGGTCGTTCGGAAAGACTTTCGCCCACGGATTGTGCGcagc                                             >  1:290547/1‑149 (MQ=255)
                                                                                                       aaGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGcagc                                             >  1:850329/1‑149 (MQ=255)
                                                                                                       aaGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGcagc                                             <  1:1116377/149‑1 (MQ=255)
                                                                                                        aGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGcagcc                                            >  2:319827/1‑149 (MQ=255)
                                                                                                          gccgACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGcagccag                                          <  1:500631/146‑1 (MQ=255)
                                                                                                          gCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCGACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGcagccag                                          <  2:145426/149‑1 (MQ=255)
                                                                                                            tGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCcagca                                        >  1:1167662/1‑149 (MQ=255)
                                                                                                            tGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCcagca                                        <  2:1127294/149‑1 (MQ=255)
                                                                                                             gACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCCCGGATTGTGCGcagccggcgg                                       >  1:863173/1‑147 (MQ=255)
                                                                                                             gACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCCCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGCTTGACCGCAACcagaag                                       >  1:158914/1‑149 (MQ=255)
                                                                                                              aCAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGt                                      <  1:980621/149‑1 (MQ=255)
                                                                                                              aCAGGAACTCTGGTCCCATATTCAGCTAACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGGTAGCTCACGCATAGGGATTTCCCGATCGTTCGGAACGAGTTTCGCCCACGGATTGTGCGCAGCACGCAGt                                      >  2:271464/1‑149 (MQ=255)
                                                                                                                 ggAACTCTGGTCCCATATTCAGCTTACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCCCCGCTTGTGCGCAGCCAGCGGttct                                   >  2:360461/1‑149 (MQ=255)
                                                                                                                      tctGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCAc                              >  2:430126/1‑149 (MQ=255)
                                                                                                                       ctGGTCCCATATTCCGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACg                             <  1:1046663/149‑1 (MQ=255)
                                                                                                                          gTCCCTTATTTAGCTGACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTc                          <  2:982529/149‑1 (MQ=255)
                                                                                                                             ccATATTCAGCTTACGCAGGCTGCCTACCAGCTCATTCAACGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAAGTCTTCCTCgcacggg                       >  2:1192594/1‑143 (MQ=255)
                                                                                                                              cATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGaa                      <  1:610148/149‑1 (MQ=255)
                                                                                                                                  ttCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGa                  <  1:319827/149‑1 (MQ=255)
                                                                                                                                     atcatACGTAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCATTTCTTCCACGGTCGGAATAGacac               <  1:972937/145‑1 (MQ=255)
                                                                                                                                     aGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTACCGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGacac               >  1:1123017/1‑149 (MQ=255)
                                                                                                                                         tACGCAGGCTGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCt           <  1:13867/149‑1 (MQ=255)
                                                                                                                                         tACGCAGGCGGCCTACCGGCGTGGTGAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCCGCCAGCAGTTCTTCCACGGTCGGGATAGACACATCt           >  1:742396/1‑149 (MQ=255)
                                                                                                                                         tACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGGTTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACCCATCt           >  1:574657/1‑149 (MQ=255)
                                                                                                                                              aGGCGGCCTACCGGCGTGGTCCGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCtttttt      <  1:1184625/149‑1 (MQ=255)
                                                                                                                                              aGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCtttttt      <  1:1184722/149‑1 (MQ=255)
                                                                                                                                                  ggcCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCTAGCTGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAAt  <  2:207674/149‑1 (MQ=255)
                                                                                                                                                  ggcCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACGTCTTTTTTCAAt  <  1:893586/149‑1 (MQ=255)
                                                                                                                                                  ggcCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAAt  <  1:1067928/149‑1 (MQ=255)
                                                                                                                                                  ggcCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAAt  <  1:989668/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGGCATCGCGCCCCAGATTTAATGAATAAAGATTACGCCAGTTGACGAAGCATCCGACGCAGCGGCTCCGCGGCCCCCCACAGCAGCTGGTCGCCCACGGTAAAGGCTGACAGGAACTCTGGTCCCATATTCAGCTTACGCAGGCGGCCTACCGGCGTGGTCAGCGTGCCGGTAACGGCAGCTGGGGTTAGCTCACGCATAGTGATTTCCCGATCGTTCGGAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAAT  >  CP009273/3567102‑3567396

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑