Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,406,553 C→A 100% R12L (CGG→CTG)  bsmA ← bioflm peroxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,406,5530CA98.9% 891.7 / 0.1 272R12L (CGG→CTG) bsmAbioflm peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (110/159);  total (110/162)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 272 total aligned reads displayed.

TGCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATG  >  CP009273/4406407‑4406698
                                                                                                                                                  |                                                                                                                                                 
tGCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAgg                                                                                                                                                 <  2:2926273/149‑1 (MQ=255)
 gCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGt                                                                                                                                                <  1:799105/149‑1 (MQ=255)
 gCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAAACAACTGGCAAACAGGt                                                                                                                                                >  1:3056258/1‑149 (MQ=255)
       gCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATaa                                                                                                                                          <  1:2096309/149‑1 (MQ=255)
        cTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAAc                                                                                                                                         <  2:186986/149‑1 (MQ=255)
         tACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACg                                                                                                                                        <  2:2202641/149‑1 (MQ=255)
          aCCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATTGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGAGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGc                                                                                                                                       <  2:61016/149‑1 (MQ=255)
          aCCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGc                                                                                                                                       <  1:3363564/149‑1 (MQ=255)
          aCCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGc                                                                                                                                       <  1:4053351/149‑1 (MQ=255)
          aCCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGc                                                                                                                                       <  1:2379326/149‑1 (MQ=255)
           ccTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCt                                                                                                                                      >  1:1090572/1‑149 (MQ=255)
             tATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAt                                                                                                                                    <  1:3981676/149‑1 (MQ=255)
             tATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCGTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAt                                                                                                                                    <  1:705644/149‑1 (MQ=255)
              aTTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAtt                                                                                                                                   >  2:3046152/1‑149 (MQ=255)
              aTTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAtt                                                                                                                                   >  2:2830440/1‑149 (MQ=255)
              aTTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAtt                                                                                                                                   <  2:1733746/149‑1 (MQ=255)
              aTTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAtt                                                                                                                                   >  1:2815944/1‑149 (MQ=255)
                   cTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGtt                                                                                                                              <  1:804962/149‑1 (MQ=255)
                   cTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGtt                                                                                                                              <  1:1433248/149‑1 (MQ=255)
                   cTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGtt                                                                                                                              >  1:2167763/1‑149 (MQ=255)
                     gTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTc                                                                                                                            <  1:307533/149‑1 (MQ=255)
                      tAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTcc                                                                                                                           <  1:993290/149‑1 (MQ=255)
                      tAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTcc                                                                                                                           <  1:2212928/149‑1 (MQ=255)
                      tAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTcc                                                                                                                           <  2:3732715/149‑1 (MQ=255)
                          ccTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCt                                                                                                                       <  2:225418/149‑1 (MQ=255)
                          ccTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCt                                                                                                                       <  2:225319/149‑1 (MQ=255)
                           cTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTa                                                                                                                      >  1:3595375/1‑149 (MQ=255)
                             tGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAAc                                                                                                                    <  1:1484468/149‑1 (MQ=255)
                              gCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAAcc                                                                                                                   <  1:1135017/149‑1 (MQ=255)
                               cGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCa                                                                                                                  <  1:2316059/149‑1 (MQ=255)
                               cGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCa                                                                                                                  <  2:3066097/149‑1 (MQ=255)
                                gTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                                                                                 <  2:3479379/149‑1 (MQ=255)
                                gTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                                                                                 <  1:133288/149‑1 (MQ=255)
                                gTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                                                                                 >  1:202917/1‑149 (MQ=255)
                                   tGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt                                                                                                              >  2:2533202/1‑149 (MQ=255)
                                   tGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt                                                                                                              >  1:1593556/1‑149 (MQ=255)
                                    ggTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtg                                                                                                             >  1:353989/1‑149 (MQ=255)
                                    ggTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtg                                                                                                             >  1:1073425/1‑149 (MQ=255)
                                     gTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgt                                                                                                            >  1:1316069/1‑149 (MQ=255)
                                     gTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgt                                                                                                            <  1:1766623/149‑1 (MQ=255)
                                        gcggcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtagt                                                                                                         >  2:1197867/1‑149 (MQ=255)
                                        gcggcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtagt                                                                                                         <  2:755675/149‑1 (MQ=255)
                                          ggcgAATCTCTTGCGGATGGTCCGTTACTGGCTGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc                                                                                                       >  1:264796/1‑149 (MQ=255)
                                          ggcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc                                                                                                       <  2:1377789/149‑1 (MQ=255)
                                           gcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGcc                                                                                                      >  2:3210926/1‑149 (MQ=255)
                                             gAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATACACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtt                                                                                                    <  1:3954382/149‑1 (MQ=255)
                                              aaTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta                                                                                                   <  2:3355205/149‑1 (MQ=255)
                                               aTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttat                                                                                                  >  1:1362756/1‑149 (MQ=255)
                                                    ttGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAtt                                                                                             <  1:3085011/149‑1 (MQ=255)
                                                    ttGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAtt                                                                                             <  1:3110184/149‑1 (MQ=255)
                                                       cGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTaa                                                                                          <  1:931844/149‑1 (MQ=255)
                                                        ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                         <  1:2689031/149‑1 (MQ=255)
                                                        ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                         <  1:285035/149‑1 (MQ=255)
                                                        ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                         <  2:2477732/149‑1 (MQ=255)
                                                          aTGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACtt                                                                                       <  1:2527257/149‑1 (MQ=255)
                                                           tGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACttt                                                                                      >  1:2982909/1‑149 (MQ=255)
                                                             gTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTgg                                                                                    >  1:306498/1‑149 (MQ=255)
                                                             gTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTgg                                                                                    <  2:663962/149‑1 (MQ=255)
                                                              tCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTggg                                                                                   <  2:1323395/149‑1 (MQ=255)
                                                                  ttACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGaa                                                                               <  1:347455/149‑1 (MQ=255)
                                                                     cTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGt                                                                            <  1:3158910/149‑1 (MQ=255)
                                                                             gctggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAAt                                                                    <  1:3319659/149‑1 (MQ=255)
                                                                             gctggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAAt                                                                    >  1:1206689/1‑149 (MQ=255)
                                                                                 gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGtt                                                                >  1:1206023/1‑149 (MQ=255)
                                                                                 gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGtt                                                                >  1:1206191/1‑149 (MQ=255)
                                                                                      ggCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTg                                                           <  2:2547300/149‑1 (MQ=255)
                                                                                         gTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTa                                                        <  2:3999463/149‑1 (MQ=255)
                                                                                          tACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAc                                                       <  1:3584471/149‑1 (MQ=255)
                                                                                            ccTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAcgc                                                     <  2:3713816/149‑1 (MQ=255)
                                                                                                gCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAtt                                                 <  2:2378606/149‑1 (MQ=255)
                                                                                                 cAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAttt                                                >  2:1708626/1‑149 (MQ=255)
                                                                                                 cAGTGCGATACAGGCACTTAACATCAACACCAATAATAAACTGGCAATCAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACTCGAttt                                                <  1:1755821/149‑1 (MQ=255)
                                                                                                  aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAtttt                                               >  2:3649317/1‑149 (MQ=255)
                                                                                                       gcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAAt                                          <  2:3565574/149‑1 (MQ=255)
                                                                                                        cTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATg                                         <  2:1986802/149‑1 (MQ=255)
                                                                                                         tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGa                                        <  2:30679/149‑1 (MQ=255)
                                                                                                               cACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGt                                  >  1:615015/1‑149 (MQ=255)
                                                                                                                 cTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGa                                >  1:3343533/1‑149 (MQ=255)
                                                                                                                    aaCATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGa                             <  2:1090572/149‑1 (MQ=255)
                                                                                                                     aCATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGaa                            >  1:2414937/1‑149 (MQ=255)
                                                                                                                      cATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGaaa                           >  2:1526318/1‑149 (MQ=255)
                                                                                                                       aTCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAg                          <  2:1889273/149‑1 (MQ=255)
                                                                                                                       aTCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAg                          <  2:734309/149‑1 (MQ=255)
                                                                                                                          aaCACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGt                       <  1:250597/149‑1 (MQ=255)
                                                                                                                           acacCAATAATAAACTGGCAAACAGTTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTTACGAAAGCGTa                      <  2:397644/149‑1 (MQ=255)
                                                                                                                                 ataataAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcac                <  2:2595578/149‑1 (MQ=255)
                                                                                                                                  taataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcaca               <  1:565921/149‑1 (MQ=255)
                                                                                                                                  taataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcaca               <  2:1063024/149‑1 (MQ=255)
                                                                                                                                  taataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcaca               <  1:2631967/149‑1 (MQ=255)
                                                                                                                                   aataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAt              <  2:2008493/149‑1 (MQ=255)
                                                                                                                                    ataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc             >  1:1426498/1‑149 (MQ=255)
                                                                                                                                    ataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc             <  1:420983/149‑1 (MQ=255)
                                                                                                                                    ataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc             <  2:535037/149‑1 (MQ=255)
                                                                                                                                    ataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc             >  1:3854755/1‑149 (MQ=255)
                                                                                                                                       aaCTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCa          >  1:1302574/1‑149 (MQ=255)
                                                                                                                                       aaCTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCa          >  1:3244712/1‑149 (MQ=255)
                                                                                                                                        aCTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCac         >  2:3808606/1‑148 (MQ=255)
                                                                                                                                           ggCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACAt      <  2:2167763/149‑1 (MQ=255)
                                                                                                                                              aaaCAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTAt   <  2:2815944/149‑1 (MQ=255)
                                                                                                                                               aaCAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATg  <  2:433433/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATG  >  CP009273/4406407‑4406698

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