Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 8,705 C→G T156T (ACC→ACG)  talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,7050CG79.7% 1408.4 / inf 548T156T (ACC→ACG) talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (56/55);  new base G (240/196);  total (297/251)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 555 total aligned reads displayed.

GTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTG  >  CP009273/8558‑8853
                                                                                                                                                   |                                                                                                                                                    
 tctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctg                                                                                                                                                   >  1:644561/7‑150 (MQ=255)
 tctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctg                                                                                                                                                   >  2:1399836/7‑150 (MQ=255)
   aaTTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCctgct                                                                                                                                                 <  1:322518/150‑1 (MQ=255)
     tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACgctgct                                                                                                                                                 <  2:661897/146‑1 (MQ=255)
     tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATC‑AACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt                                                                                                                                              >  2:739347/3‑150 (MQ=255)
      cgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt                                                                                                                                              >  1:2211659/2‑150 (MQ=255)
       gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt                                                                                                                                              <  1:2317693/149‑1 (MQ=255)
        cGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTGCCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtct                                                                                                                                            <  2:1028215/150‑1 (MQ=255)
           aaGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCt                                                                                                                                         <  2:1500245/150‑1 (MQ=255)
                   cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                                                                                                                                 <  1:2102751/144‑1 (MQ=255)
                   cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                                                                                                                                 <  1:1406424/144‑1 (MQ=255)
                     catcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGc                                                                                                                               >  2:557322/5‑150 (MQ=255)
                        gATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGAAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCg                                                                                                                            >  2:1190423/1‑150 (MQ=255)
                         aTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCgc                                                                                                                           >  2:2369711/1‑149 (MQ=255)
                         aTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCgc                                                                                                                           >  1:1992933/1‑149 (MQ=255)
                         aTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCg                                                                                                                            <  2:1621524/149‑1 (MQ=255)
                         aTCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTATCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTc                                                                                                                             <  2:2010740/148‑1 (MQ=255)
                         aTCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCg                                                                                                                            <  1:401991/149‑1 (MQ=255)
                            aagctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCt                                                                                                                        <  2:1746246/144‑1 (MQ=255)
                             aaCTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCtt                                                                                                                       <  1:373567/150‑1 (MQ=255)
                                tgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgcg                                                                                                                    >  1:1343437/3‑148 (MQ=255)
                                  tACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgcg                                                                                                                    <  2:818595/148‑3 (MQ=255)
                                  tACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgcc                                                                                                                  <  1:1524095/150‑2 (MQ=255)
                                    cAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgcgccg                                                                                                                 >  2:1983015/1‑144 (MQ=255)
                                    cAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgccg                                                                                                                 <  2:1799236/149‑3 (MQ=255)
                                            cTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCg                                                                                                        >  2:780788/1‑150 (MQ=255)
                                              ggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgtc                                                                                                      <  2:144516/144‑2 (MQ=255)
                                                 atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                                                                                   >  1:1177192/4‑150 (MQ=255)
                                                    aGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCt                                                                                                 <  1:270654/149‑1 (MQ=255)
                                                     gCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGa                                                                                               >  1:2178479/1‑150 (MQ=255)
                                                      cAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAt                                                                                              <  2:2152553/150‑1 (MQ=255)
                                                       aaCGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAtt                                                                                             <  2:2223144/150‑2 (MQ=255)
                                                        aCGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGAtct                                                                                            >  2:1063760/1‑150 (MQ=255)
                                                           aTCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGc                                                                                         >  1:1269618/1‑150 (MQ=255)
                                                              gcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGt                                                                                      <  2:444302/148‑1 (MQ=255)
                                                              gTATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCg                                                                                       <  1:1936227/149‑1 (MQ=255)
                                                               caTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGtt                                                                                     >  2:613942/2‑150 (MQ=255)
                                                                aTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGttc                                                                                    >  1:1495553/1‑149 (MQ=255)
                                                                   cTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTttgtt                                                                                 >  2:1931543/1‑150 (MQ=255)
                                                                      aTCAAACTGGCTTCTACCTGACAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTgg                                                                               <  1:2293748/149‑1 (MQ=255)
                                                                      aTCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTtcgttg                                                                                <  2:154056/148‑1 (MQ=255)
                                                                          aaCTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTa                                                                          >  1:932540/1‑150 (MQ=255)
                                                                           aCTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCcgttcgtt                                                                                 >  2:156356/1‑142 (MQ=255)
                                                                             tGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTAtcc                                                                       >  2:1130305/1‑148 (MQ=255)
                                                                              ggCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCt                                                                      >  1:1690203/1‑150 (MQ=255)
                                                                               gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                     <  1:599959/144‑2 (MQ=255)
                                                                               gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                     <  2:516714/144‑2 (MQ=255)
                                                                               gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                     <  1:172451/144‑2 (MQ=255)
                                                                               gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                     <  2:487115/144‑2 (MQ=255)
                                                                                 ttCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGa                                                                   >  1:2349386/1‑150 (MQ=255)
                                                                                  tccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAc                                                                  >  1:883550/4‑150 (MQ=255)
                                                                                    caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAc                                                                  <  2:2109403/147‑1 (MQ=255)
                                                                                    caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTg                                                                <  2:384255/149‑1 (MQ=255)
                                                                                    caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTg                                                                >  2:703672/2‑150 (MQ=255)
                                                                                    caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTg                                                                >  2:283273/2‑150 (MQ=255)
                                                                                     aCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTgg                                                               >  2:24537/1‑150 (MQ=255)
                                                                                     aCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTgg                                                               >  1:42289/1‑150 (MQ=255)
                                                                                     aCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTgg                                                               <  1:1020905/150‑1 (MQ=255)
                                                                                      ccTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTATTGACTg                                                                <  2:1859153/148‑1 (MQ=255)
                                                                                       cTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTa                                                             <  1:674004/150‑1 (MQ=255)
                                                                                           cAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACaaa                                                         <  2:205951/150‑1 (MQ=255)
                                                                                                      cgccgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGa                                              >  2:2208003/8‑150 (MQ=255)
                                                                                                       gccgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGa                                              >  2:1331967/7‑149 (MQ=255)
                                                                                                        ccgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGt                                                              <  1:1860830/127‑1 (MQ=255)
                                                                                                         cgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGac                                             >  2:2284240/5‑147 (MQ=255)
                                                                                                          gctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGac                                             >  1:622889/4‑146 (MQ=255)
                                                                                                           ctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACaag                                          >  2:408896/3‑149 (MQ=255)
                                                                                                              aaCAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGATTGGTACAAAGCCAACACCGACaaga                                         >  2:1295322/1‑147 (MQ=255)
                                                                                                               aCAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGGGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACaagaaaga                                     <  2:2335872/150‑1 (MQ=255)
                                                                                                               aCAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACaagaaaga                                     <  2:904029/150‑1 (MQ=255)
                                                                                                                cAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAagag                                    <  1:1546203/150‑1 (MQ=255)
                                                                                                                cAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAagag                                    >  1:2240473/1‑150 (MQ=255)
                                                                                                                cAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAagag                                    <  1:1779001/150‑1 (MQ=255)
                                                                                                                   cTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTAc                                 <  1:1082121/150‑1 (MQ=255)
                                                                                                                   cTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTAc                                 <  1:2407334/150‑1 (MQ=255)
                                                                                                                      gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACg                                <  2:115448/148‑1 (MQ=255)
                                                                                                                      gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                              >  1:513451/1‑149 (MQ=255)
                                                                                                                        aaaaGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTcc                            <  2:1257104/150‑1 (MQ=255)
                                                                                                                         aaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCg                           >  1:190805/1‑150 (MQ=255)
                                                                                                                           agaagGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCgg                          >  1:1235417/1‑149 (MQ=255)
                                                                                                                             aagGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                       <  1:334931/150‑1 (MQ=255)
                                                                                                                                gCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                       >  1:7591/1‑147 (MQ=255)
                                                                                                                                gCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                       >  1:1310550/1‑147 (MQ=255)
                                                                                                                                gCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCag                       >  2:1875229/1‑147 (MQ=255)
                                                                                                                                 cATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                       >  1:1826491/1‑146 (MQ=255)
                                                                                                                                 cATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                       >  2:1361104/1‑146 (MQ=255)
                                                                                                                                  attaaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCGTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGAt                  <  1:116886/147‑1 (MQ=255)
                                                                                                                                    taaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATcc                >  1:824841/2‑150 (MQ=255)
                                                                                                                                    taaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATcc                >  1:971778/2‑150 (MQ=255)
                                                                                                                                    taaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATcc                >  1:1455304/2‑150 (MQ=255)
                                                                                                                                     aaCTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATcc                <  1:477290/149‑1 (MQ=255)
                                                                                                                                       cTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCggg             >  1:1588128/1‑150 (MQ=255)
                                                                                                                                          tAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGt          >  2:310870/1‑150 (MQ=255)
                                                                                                                                           aaCCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTg         >  2:1667496/1‑150 (MQ=255)
                                                                                                                                            aCCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTg         <  2:7591/149‑1 (MQ=255)
                                                                                                                                            aCCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTg         <  1:1651872/149‑1 (MQ=255)
                                                                                                                                             ccTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGt       <  2:1175874/150‑1 (MQ=255)
                                                                                                                                               tGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGttt     >  2:856632/1‑150 (MQ=255)
                                                                                                                                               tGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGttt     >  2:470578/1‑150 (MQ=255)
                                                                                                                                                  cccTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTg  >  1:1380956/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTG  >  CP009273/8558‑8853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑