Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,759 T→G 7.0% V232G (GTT→GGT)  metK → S‑adenosylmethionine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,080,7590TG 7.0% 586.7 / 7.0 237V232G (GTT→GGT) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (88/125);  new base G (21/2);  total (110/127)
Fisher's exact test for biased strand distribution p-value = 2.93e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATG  >  CP009273/3080610‑3080908
                                                                                                                                                     |                                                                                                                                                       
tGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgt                                                                                                                                                         <  1:3410304/150‑1 (MQ=255)
    gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc                                                                                                                                                     <  2:1246589/150‑1 (MQ=255)
    gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc                                                                                                                                                     <  2:3307286/150‑1 (MQ=255)
    gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc                                                                                                                                                     <  2:4024065/150‑1 (MQ=255)
       gagCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGt                                                                                                                                                 <  2:1162361/148‑1 (MQ=255)
       gagCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGt                                                                                                                                                 >  1:4199041/3‑150 (MQ=255)
        tGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTg                                                                                                                                               >  1:4259067/1‑150 (MQ=255)
          cTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTCAGGGAAGAGATCATCAAGCCACTTCTGCCCGCTGACTGGCTGACTTCTGCCACCAAAGTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGc                                                                                                                                             >  2:348863/1‑150 (MQ=255)
           tttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGcc                                                                                                                                            >  2:1204446/1‑150 (MQ=255)
             tCCACTCATCACTCTGAAGAGATCGACCATAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCa                                                                                                                                          <  2:3014587/150‑1 (MQ=255)
                 cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAGTTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATgg                                                                                                                                      <  2:1281454/150‑1 (MQ=255)
                  tCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATggg                                                                                                                                     >  1:2327622/1‑150 (MQ=255)
                   cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGt                                                                                                                                    <  2:878481/150‑1 (MQ=255)
                    aGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAACTTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGGGGCCCACAGGGTg                                                                                                                                   >  1:2493068/1‑150 (MQ=255)
                     gCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATAGGTGGCCCACAGGGTgg                                                                                                                                  >  1:3907383/1‑149 (MQ=255)
                       aCTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACt                                                                                                                                >  2:4238032/1‑150 (MQ=255)
                       aCTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTTTCTGTTGCCCCAAGGGTGAcg                                                                                                                                >  1:1190611/1‑149 (MQ=255)
                        ctctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTg                                                                                                                               <  1:4355436/150‑1 (MQ=255)
                         tctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGc                                                                                                                              <  1:401035/150‑1 (MQ=255)
                          ctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCg                                                                                                                             <  1:1721568/150‑1 (MQ=255)
                           tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCATGGGTGACTGCgg                                                                                                                            >  2:218250/1‑150 (MQ=255)
                            gAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                           <  2:4264881/150‑1 (MQ=255)
                              agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                         <  1:3458500/150‑1 (MQ=255)
                               gagaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTg                                                                                                                        <  2:3326047/150‑1 (MQ=255)
                               gagaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTg                                                                                                                        <  2:4097775/150‑1 (MQ=255)
                                 gaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGAc                                                                                                                      <  1:2108810/150‑1 (MQ=255)
                                  aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTTTCGGTGGGCCCATGGGGGACTGCGGTCTGGCt                                                                                                                     >  2:1617887/1‑150 (MQ=255)
                                    cGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGTCTGCGGTGTGTCTgg                                                                                                                   >  1:1110669/1‑150 (MQ=255)
                                      aCCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTc                                                                                                                 <  2:2820496/150‑1 (MQ=255)
                                       ccAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCg                                                                                                                <  2:4521891/150‑1 (MQ=255)
                                          gAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTTTCGGGGGCCCAATGGGTGACTGCGGTCTGGCTGGTCGTaa                                                                                                             >  1:1091898/1‑150 (MQ=255)
                                           aaaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGGGGCCCCAAGGGGGGCTGCGGGCTGGCGGGGCGTcaa                                                                                                            >  1:4551504/1‑150 (MQ=255)
                                            aaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATAGGTGGCCCCAAGGGTGGCTGCGGGCTGACTGGGGGTaaaa                                                                                                           >  1:3673520/1‑150 (MQ=255)
                                              tCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCAGACTGGTCGGAAAAtt                                                                                                         >  2:601998/1‑150 (MQ=255)
                                               cGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGCCCCAAGGGGGTCGGCGGTGTGGGGGGTGGTAAAAATa                                                                                                        >  1:848290/1‑150 (MQ=255)
                                                gCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTAt                                                                                                       <  1:1901894/150‑1 (MQ=255)
                                                    cAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGtt                                                                                                   <  2:3608384/150‑1 (MQ=255)
                                                     aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCATGGGTGACGGCGGGATGACTGGTCGTAAAAATATAGttt                                                                                                  >  2:2407932/1‑149 (MQ=255)
                                                     aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTg                                                                                                  >  2:712062/1‑150 (MQ=255)
                                                     aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTg                                                                                                  <  1:4757107/150‑1 (MQ=255)
                                                        aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAAAGGGTGACGGCGGGCTGACTGGGCGTAAAATTATAGTTGATa                                                                                               >  1:3499537/1‑150 (MQ=255)
                                                           cGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCt                                                                                            <  2:2910577/150‑1 (MQ=255)
                                                           cGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGGTGATAACt                                                                                            >  2:4662038/1‑150 (MQ=255)
                                                              tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGTCTGCGGTCTGACTGGTCGGAAAAATATAGTTTTTACCTAcg                                                                                         >  1:4722310/1‑150 (MQ=255)
                                                               aaTGGAAGTGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGCTACCTAcgg                                                                                        >  2:3143135/1‑150 (MQ=255)
                                                                  ggAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcggcgg                                                                                     <  1:4430101/150‑1 (MQ=255)
                                                                    aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCa                                                                                   >  2:1720551/1‑150 (MQ=255)
                                                                      gagaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATg                                                                                 <  1:1747041/150‑1 (MQ=255)
                                                                       agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATgg                                                                                <  2:1000768/150‑1 (MQ=255)
                                                                       agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATgg                                                                                <  1:4292746/150‑1 (MQ=255)
                                                                       agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTTTAGGTTGGCCCAAGGGTGTCTGCGGTGGGTCTGGTGGTTAAATTATAGTTGTTAACAAAGGCGGCCTgg                                                                                >  1:1633832/1‑150 (MQ=255)
                                                                        gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGGCGTAAAATTAACGTTGATACCTACGGCGGCATGgc                                                                               >  2:1138433/1‑150 (MQ=255)
                                                                          tcatcaAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGGGGGTTTAGGTGGGCCCAATGGTGGCGGCGGTCTGGCGGGTGGGGAAATTATAGGTGTTAACAAACGCGGCAAGgcgc                                                                             >  1:3985244/1‑150 (MQ=255)
                                                                             tcaAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTc                                                                          <  1:621435/150‑1 (MQ=255)
                                                                              caAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCa                                                                         <  2:1878727/150‑1 (MQ=255)
                                                                               aaGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCAc                                                                        >  2:1956890/1‑150 (MQ=255)
                                                                                 gCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACgg                                                                      <  1:2941142/150‑1 (MQ=255)
                                                                                    aaTTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTgg                                                                   <  1:3933667/150‑1 (MQ=255)
                                                                                    aaTTCTGCCCGCTGAATGGCTGACTTCGGCCACCACCGTCTTCATCAACCCGACCGGTCGTGTCGTGATCGGGGGCCCACGGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTGGATACCTACGGCGGCATGGCGCGTCACGGTgc                                                                   <  1:3384610/150‑2 (MQ=255)
                                                                                     aTTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGc                                                                  <  2:1893776/150‑1 (MQ=255)
                                                                                      ttCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCg                                                                 <  2:257174/150‑1 (MQ=255)
                                                                                       tCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTTTAGGTGGGCCAATGGGTGACTGCGGTGTGGCTGGTGGTAAAATTATAGTTTTTACCTAAGGCGGGAATGGGCGTCACGGTGgggg                                                                >  2:2891558/1‑150 (MQ=255)
                                                                                       tCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTATAGGTGGCCCCAAGGGTGGGTGCGGTGTGAAGGGTCGTAAAAATATAGTTTTTAACAAAGGCGGCCAGTCGCGTCACGGTGgggg                                                                >  1:619143/1‑150 (MQ=255)
                                                                                        cTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGt                                                               <  1:2180400/150‑1 (MQ=255)
                                                                                        cTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGGGGCCCAATGGGTGTCTGCGGTCTGACGGGTCGTAAAATTATAGGTGATACCTACGGCGGCAAGGCGCGGCACGGTGGCGGt                                                               >  2:3417423/1‑150 (MQ=255)
                                                                                         tGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAAGGCGGCATGGCGCGTCACGGTGGCGGTg                                                              >  2:3898950/1‑150 (MQ=255)
                                                                                           cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCAAGGGTGTCTGCGGTCTGTCTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                            >  1:1961547/1‑150 (MQ=255)
                                                                                               cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtct                                                        <  1:2220968/150‑1 (MQ=255)
                                                                                               cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtct                                                        <  1:704742/150‑1 (MQ=255)
                                                                                                tGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGACTGCGGTCTGAATGGTCGTAAAATTATAGTTGATAACAAAGGCGGCATGGCGCGTCACGGTGGCGGTGCATtctc                                                       >  1:2619924/1‑150 (MQ=255)
                                                                                                tGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTAGTTATAGGTGGGCCCCATGGTGTCTGCGGTCTGCCTGGTCGTAAAAATATAGTTTATAACAACGGCGGCAAGGCGCGTCACGGTGGCGGTGCATtctc                                                       >  2:1577329/1‑150 (MQ=255)
                                                                                                tGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGGCCCAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGTTACCTAAGGCGGCATGGCGCGTCACCGTGGCGGTGCAttata                                                       >  1:3764096/1‑147 (MQ=255)
                                                                                                    tgacTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTGATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGt                                                   <  2:316163/147‑1 (MQ=255)
                                                                                                         gACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTTTAGGTGGGCCAAGTGGTGTATGCGGTGTGGAGGGTCGGTAAAATATCGTTTTTACCTACGGGGGGATGGCGCGTCACGGTGGGGGTGCATTATATGGTTAAGa                                              >  2:3987978/1‑150 (MQ=255)
                                                                                                           cTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAGGTGGCCCAATGGGGGACAGCGGTGTTTCATGGCGTAAAATTATAGTTGATAAAAAACGCGGCCATGCGCGTCACCGGGGCGGTGCATTATCTTGTtaataaa                                            >  1:2980680/1‑146 (MQ=255)
                                                                                                            ttCTGCCAACAAATTCTTCATAAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTAGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATcc                                           <  1:3214319/150‑1 (MQ=255)
                                                                                                             tCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAAAGGGTGTCTGCGGTCTGACTGGTCGTTAAATTATCGTTGATACCAAACGCCGCAAGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAATATACa                                          >  2:1385353/1‑150 (MQ=255)
                                                                                                             tCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCTGGTCGTAAAATTATAGTTGATACCTAAGGCGGCAATGCGCGTCACGGTGGCGGTGCATTCTATGGTAAAGATCCa                                          >  2:250886/1‑150 (MQ=255)
                                                                                                              cTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTAGTTTTAGGTGGCCCAATGGGTGACTGCGGTCTTACTGGTAGTAAAATTATATTTGTTAACAAACGCGGCATGGCGCGTGACCGTGGCGGTGCATTATCTTGTTAatataaat                                         >  1:3099242/1‑146 (MQ=255)
                                                                                                               tGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAAATGGTGTCTGCGGTCTGACTGGTCGTAAAAATATAGTTGATAACAACGGCGGCATTGCGCGTCACCGTGGCGGTGCATTCTATGGTAAatataaaaa                                        >  1:3115893/1‑145 (MQ=255)
                                                                                                                 ccaccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTTGCCCAAAGGGTGACTGCGGTCTTACTGGTCGTAAAAATATAGTTTTTACCAAAGGCGGCAATGCGCGTCCAGGTTGCGGTGGATTATCTGGTAAAAATACCTCaa                                      >  2:51863/1‑150 (MQ=255)
                                                                                                                     caAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                                  <  2:3555548/150‑1 (MQ=255)
                                                                                                                           cttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCg                            <  2:2021415/150‑1 (MQ=255)
                                                                                                                           cttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCg                            <  2:4425757/150‑1 (MQ=255)
                                                                                                                              catcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                         <  1:3603524/150‑1 (MQ=255)
                                                                                                                               atcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTcc                        <  2:4355231/150‑1 (MQ=255)
                                                                                                                                tcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCg                       <  2:2128947/150‑1 (MQ=255)
                                                                                                                                  aaCCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCa                     <  1:2794/150‑1 (MQ=255)
                                                                                                                                     ccGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGcc                  <  1:1566904/150‑1 (MQ=255)
                                                                                                                                     ccGACCGGTCGTGTCGGTATAGGTGGCCCCAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTTATACCTACGGCGGCCAGGCGCGTCACGGTGGCGGTGCAATCACTGGTAAACATACCAAAAAAGTGGACCCTTCCCCAGcc                  >  2:1067777/1‑150 (MQ=255)
                                                                                                                                      cGACCGGTCGTTTCGTGATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                 <  1:2894273/150‑1 (MQ=255)
                                                                                                                                          cGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgc             >  2:2625175/1‑150 (MQ=255)
                                                                                                                                           ggTCGTTTCGTTATCGGTGGCCCAAAGGGTGAATGCGGTCTGACTGGTCGTAAAAATATAGTTTATAACAACGGCGGCATTGCGCGTAACCGTGGCGGTGCAATATCTGGTAAATATACAAAAAAATTGTACCCTTACGCAGccaaacca            >  1:3013427/1‑144 (MQ=255)
                                                                                                                                            gTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgcag           <  1:4246277/150‑1 (MQ=255)
                                                                                                                                                tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACg       >  2:1301913/1‑150 (MQ=255)
                                                                                                                                                tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACg       <  1:1518454/150‑1 (MQ=255)
                                                                                                                                                 ttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt      >  1:101444/1‑150 (MQ=255)
                                                                                                                                                    gttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTAt   <  1:1065203/150‑1 (MQ=255)
                                                                                                                                                     ttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATg  <  2:3673520/150‑1 (MQ=255)
                                                                                                                                                     ttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATg  >  2:862892/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATG  >  CP009273/3080610‑3080908

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