Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,823,835 A→C 7.3% W53G (TGG→GGG)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,823,8350AC 7.3% 647.8 / 2.2 284W53G (TGG→GGG) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (161/86);  new base C (3/27);  total (164/120)
Fisher's exact test for biased strand distribution p-value = 4.65e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

CATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCAT  >  CP009273/1823689‑1823983
                                                                                                                                                    |                                                                                                                                                    
  tcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcgc                                                                                                                                                  >  1:3317757/1‑149 (MQ=255)
    gCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCgg                                                                                                                                                >  2:3259692/1‑149 (MQ=255)
     ccGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGa                                                                                                                                               >  1:3074837/1‑149 (MQ=255)
     ccGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGa                                                                                                                                               <  1:575612/149‑1 (MQ=255)
      cGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGaa                                                                                                                                              >  1:3738366/1‑149 (MQ=255)
       gTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGaaa                                                                                                                                             >  2:4681397/1‑149 (MQ=255)
        tCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAc                                                                                                                                            <  2:2978138/149‑1 (MQ=255)
         cACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcg                                                                                                                                           >  1:2611866/1‑149 (MQ=255)
          aCTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgc                                                                                                                                          <  2:3809641/149‑1 (MQ=255)
          aCTTCGGTTGCCGCTGCCCCGCGCGGCTTACCCGTTTCTCGGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACACCGGCATGACGTTCCGCAAATCAAAGCCGCGCCCCGCGCGGAAAcgc                                                                                                                                          <  1:2152867/149‑1 (MQ=255)
           cTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCAGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgcc                                                                                                                                         >  2:4709788/1‑149 (MQ=255)
              cGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCa                                                                                                                                      >  1:467704/1‑149 (MQ=255)
                gTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACg                                                                                                                                    >  2:866212/1‑149 (MQ=255)
                gTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACg                                                                                                                                    >  1:1256993/1‑149 (MQ=255)
                         tCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCg                                                                                                                           >  2:2823608/1‑149 (MQ=255)
                         tCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCg                                                                                                                           >  1:3923909/1‑149 (MQ=255)
                         tCCCAGCGCGGCTAACCCGTTTCTCTGGAAATAATCGCGTTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGATTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCg                                                                                                                           <  2:571900/149‑1 (MQ=255)
                          cccAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGa                                                                                                                          >  1:2666577/1‑149 (MQ=255)
                          cccAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGa                                                                                                                          >  1:2994856/1‑149 (MQ=255)
                          cccAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGa                                                                                                                          >  1:1547727/1‑149 (MQ=255)
                                gcgGATAACCCTTTTCTCTGGCAATAATCGGGGTTATTTCGGCTTTATTGCTTTCCTGCTGTGCGGCAAAGCGTTCGAAAACGGCATTACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGcc                                                                                                                    <  1:4366188/149‑1 (MQ=255)
                                  ggCTTACCCGTTTCTCTGGCAATAATCGCGTTTAATTCGGCTTTATTGCTTTCCAGCTGTGCGGCAACGGGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTg                                                                                                                  <  2:1004454/149‑1 (MQ=255)
                                       aCCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGa                                                                                                             >  1:3684519/1‑149 (MQ=255)
                                        cccGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                                            >  1:2029492/1‑149 (MQ=255)
                                            tttCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGa                                                                                                        >  1:4131238/1‑149 (MQ=255)
                                            tttCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGa                                                                                                        <  2:2769950/149‑1 (MQ=255)
                                                 ctGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggc                                                                                                   >  1:1652362/1‑149 (MQ=255)
                                                 ctGGAAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAACGCGTTCGACAACGGAATGACGTTCAGAAAATTAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggc                                                                                                   <  1:1188689/149‑1 (MQ=255)
                                                   ggatataatCGCGTTTATTTGGTCTTTTTTGCTGTCCAGCAGTGCGGAAATGGGTTGGACAACGGATTGACTTTCAGAATTTAAAAGCGGCGCCCTGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcat                                                                                                 <  2:3998310/145‑1 (MQ=255)
                                                   ggCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcat                                                                                                 >  2:974937/1‑149 (MQ=255)
                                                   ggCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcat                                                                                                 >  1:775719/1‑149 (MQ=255)
                                                    gCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatc                                                                                                >  1:1339210/1‑149 (MQ=255)
                                                       ataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggc                                                                                             >  2:3202216/1‑149 (MQ=255)
                                                           gcgcgGTTAATGCGGCTTTTGTGCTTGCCAGCGGTGCGGCAACGGGTTGTCCACCGGCATGAGTTACCCAAAATTAACCCCGCGCCCCGCGGGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGcatca                                                                                         <  1:637570/148‑1 (MQ=255)
                                                             ggggTTAATGCGGCTTTTTTGGTTTCCAGCTGTGCGGCAAAGGTTTCGACACGGGAATGACTTACACCAAAGAAAAGCCGCGCCCCGCGGGCAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                                                                       <  1:4191413/147‑1 (MQ=255)
                                                               ggTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgc                                                                                     >  1:2747956/1‑149 (MQ=255)
                                                                tttATTTCGGTTTTTTTGCTTTCCTGCTGTGCGCAAACGGTTGCGACACCGGAATGACGTCCACCAATTCAACCCCGCCCCCCCCGCGCAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcc                                                                                    <  1:3722989/148‑1 (MQ=255)
                                                                gTTATTTGGGCTTTTTGGCTTCCCCGCTGTGCGGAAACGCTTTCGAAACGGCATTGACGTTCAGCAAATCAAAGCCGCGCCCCGCGCGCACACGCCGCACTGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcc                                                                                    <  1:2270583/149‑1 (MQ=255)
                                                                   aaTTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttg                                                                                 >  1:4244639/1‑149 (MQ=255)
                                                                     ttCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttgcc                                                                               >  1:4450550/1‑149 (MQ=255)
                                                                       cGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCAt                                                                             >  2:1678102/1‑149 (MQ=255)
                                                                         gCTTTCTTGCTTTCCCGCAGTGCGGCAACGCGTTCGACACCGGCATGACTTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                                           <  2:4449799/149‑1 (MQ=255)
                                                                         gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                                           >  1:2342835/1‑149 (MQ=255)
                                                                         gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                                           >  2:4182636/1‑149 (MQ=255)
                                                                         gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                                           >  2:3723153/1‑149 (MQ=255)
                                                                           tttATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAaca                                                                         >  1:2418082/1‑149 (MQ=255)
                                                                            tttttGCTTTCCAGCTGGGCGGCAAAGCGTTCGACAACGGCTTGACGTTCCGCAAATGAAAGCCGCGCCCCGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTAGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAacac                                                                        <  2:1337742/146‑1 (MQ=255)
                                                                             tATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACAcc                                                                       >  2:219495/1‑149 (MQ=255)
                                                                               ttGCTTTCCAGCAGTGCGGAAAAGCGTACGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGTAATCGGCATCATTGCCTTGCCATAACACCTc                                                                     <  2:1190770/149‑1 (MQ=255)
                                                                                   tttCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGccc                                                                 >  1:817692/1‑149 (MQ=255)
                                                                                   tttCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGccc                                                                 >  2:2664833/1‑149 (MQ=255)
                                                                                      ccAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGAt                                                              >  2:738181/1‑149 (MQ=255)
                                                                                       cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                                                             >  2:4210091/1‑149 (MQ=255)
                                                                                       cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                                                             >  2:1819445/1‑149 (MQ=255)
                                                                                       cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTAGCCATAACACCTCGCCCGATa                                                             >  2:2993721/1‑149 (MQ=255)
                                                                                        agcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc                                                            >  1:882420/1‑149 (MQ=255)
                                                                                        agcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc                                                            <  2:1991551/149‑1 (MQ=255)
                                                                                         gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAcc                                                           <  2:3452476/149‑1 (MQ=255)
                                                                                         gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAcc                                                           >  1:1784410/1‑149 (MQ=255)
                                                                                         gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAcc                                                           <  2:2567490/149‑1 (MQ=255)
                                                                                            gTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGa                                                        >  2:2162363/1‑149 (MQ=255)
                                                                                             tGCGGCAATGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAt                                                       <  2:3702356/149‑1 (MQ=255)
                                                                                                ggCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTAc                                                    >  2:2094098/1‑149 (MQ=255)
                                                                                                ggCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTAc                                                    <  2:1540676/149‑1 (MQ=255)
                                                                                                  aaaaGCTTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGc                                                  <  2:3871560/148‑1 (MQ=255)
                                                                                                    aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCtt                                                >  2:4146399/1‑149 (MQ=255)
                                                                                                     aGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc                                               >  2:4099089/1‑149 (MQ=255)
                                                                                                        gTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcaca                                            >  1:2356808/1‑149 (MQ=255)
                                                                                                         ttCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCacac                                           <  2:4193361/149‑1 (MQ=255)
                                                                                                         ttCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCacac                                           >  2:3132160/1‑149 (MQ=255)
                                                                                                          tCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAGTACGCTTCACAcg                                          >  1:47661/1‑149 (MQ=255)
                                                                                                               cAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGAt                                      <  2:1256993/149‑1 (MQ=255)
                                                                                                               cAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGAt                                      >  1:4587292/1‑149 (MQ=255)
                                                                                                               cAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGAt                                      >  2:3291170/1‑149 (MQ=255)
                                                                                                                accGGAATGACGTTAAGAATATTAAAGCGGCGCCCCGGGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATg                                     <  1:4699229/147‑1 (MQ=255)
                                                                                                                 ccGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATgc                                    <  2:4750618/148‑1 (MQ=255)
                                                                                                                 aCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATgc                                    >  1:3410313/1‑149 (MQ=255)
                                                                                                                      aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCt                               <  1:1116519/149‑1 (MQ=255)
                                                                                                                       ttGACGTTCAGCAAATGAAAGCCGCGCCCTGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCt                               <  1:3503807/148‑1 (MQ=255)
                                                                                                                         gACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTgg                             <  2:4362222/149‑1 (MQ=255)
                                                                                                                              tCAGCAAATGAAAGCCGCGCCCTGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGt                        <  2:3228362/149‑1 (MQ=255)
                                                                                                                                aGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTa                      <  2:4391466/149‑1 (MQ=255)
                                                                                                                                 gCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTAt                     >  2:2377373/1‑149 (MQ=255)
                                                                                                                                    aaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCa                  >  1:4029739/1‑149 (MQ=255)
                                                                                                                                    aaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCa                  >  1:482180/1‑149 (MQ=255)
                                                                                                                                     aTGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAg                 >  2:2939908/1‑149 (MQ=255)
                                                                                                                                      tGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGt                <  1:1109278/149‑1 (MQ=255)
                                                                                                                                      tGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGt                >  1:2801701/1‑149 (MQ=255)
                                                                                                                                      tGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGt                >  2:1072397/1‑149 (MQ=255)
                                                                                                                                       gAAAGCCGCGCCCAGCGCGTACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTc               <  1:3259550/149‑1 (MQ=255)
                                                                                                                                       gAAAGCCGCGCCCAGCGCGTAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTc               <  1:3259681/149‑1 (MQ=255)
                                                                                                                                         aaGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCAc             <  1:1837563/149‑1 (MQ=255)
                                                                                                                                          aGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCAcc            >  2:2029928/1‑149 (MQ=255)
                                                                                                                                           gCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACCAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCg           <  2:3551601/149‑1 (MQ=255)
                                                                                                                                           gCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCg           <  2:195963/149‑1 (MQ=255)
                                                                                                                                           gCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCg           <  2:499808/149‑1 (MQ=255)
                                                                                                                                           gCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCg           <  2:3551638/149‑1 (MQ=255)
                                                                                                                                               cgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa       <  2:4518195/149‑1 (MQ=255)
                                                                                                                                               cgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATAATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa       <  2:3235995/149‑1 (MQ=255)
                                                                                                                                                  ccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTCATcc    >  2:1303272/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCAT  >  CP009273/1823689‑1823983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: