Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,940,829 A→T 9.2% intergenic (+7/‑27) metZ → / → metW tRNA‑Met/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,940,8290AT 9.2% 441.2 / 4.8 167intergenic (+7/‑27)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (52/94);  new base T (6/12);  total (59/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

TCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAG  >  CP009273/2940681‑2940977
                                                                                                                                                    |                                                                                                                                                      
tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                                       <  1:3439586/150‑1 (MQ=255)
tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                                       >  1:4745828/1‑150 (MQ=255)
  cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                                     >  1:787097/1‑150 (MQ=255)
   gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTtg                                                                                                                                                    >  2:2317273/1‑150 (MQ=255)
     gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                                  <  1:36688/150‑1 (MQ=255)
     gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                                  <  1:4566945/150‑1 (MQ=255)
     gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTtgat                                                                                                                                                  >  1:2745658/1‑150 (MQ=255)
         ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAg                                                                                                                                             <  1:2023295/150‑1 (MQ=255)
         ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAg                                                                                                                                             <  1:420135/150‑1 (MQ=255)
         ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAg                                                                                                                                             <  1:1395823/150‑1 (MQ=255)
           ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑Ta                                                                                                                                          <  2:568671/150‑1 (MQ=255)
            gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGGAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGTTG‑AAG‑Taa                                                                                                                                         >  2:3351313/1‑150 (MQ=255)
            gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑Taa                                                                                                                                         <  2:408339/150‑1 (MQ=255)
            gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑Taa                                                                                                                                         <  2:409095/150‑1 (MQ=255)
             tAAAAAAAAGGTTGAATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑Taaa                                                                                                                                        <  1:2698950/150‑1 (MQ=255)
              aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATG‑AAG‑Taaaag                                                                                                                                      <  2:1809152/150‑3 (MQ=255)
                aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCa                                                                                                                                     >  2:383451/1‑150 (MQ=255)
                aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCa                                                                                                                                     <  1:4391418/150‑1 (MQ=255)
                 aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAg                                                                                                                                    <  2:3317633/150‑1 (MQ=255)
                 aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTTCACGCAAG‑AAAAGCAg                                                                                                                                    >  2:3290428/1‑150 (MQ=18)
                  aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTGGTTG‑ATT‑AAAAGCAGt                                                                                                                                   >  2:171257/1‑150 (MQ=18)
                   aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTa                                                                                                                                  <  2:3926071/150‑1 (MQ=255)
                   aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCCAAGCTCGTCTGATCCAATCCCCCACCCGCAACCAATGAAATTGGTATT‑AGG‑AAAAGCAGTa                                                                                                                                  >  1:1052879/1‑150 (MQ=18)
                   aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACAACACTGTCGCCGCTTCAATTCAGACCCACGCAACCAATTCAAATATCAGG‑AGG‑TAGTGCAGTa                                                                                                                                  >  2:2719661/1‑150 (MQ=21)
                    aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCACTTAAATTTTGATG‑AAG‑TAAAGCAGTAc                                                                                                                                 >  1:3096278/1‑150 (MQ=255)
                    aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGG‑AAG‑TAAAGCAGTAc                                                                                                                                 >  2:1476377/1‑150 (MQ=255)
                     aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCACATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCCCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACg                                                                                                                                <  1:3673607/150‑1 (MQ=255)
                       gTTTCATGAAAACGCGAGCGGAGTATAGTGCTCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGATCATAACCAGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTTGATG‑AAG‑TAAAGCAGTACGGt                                                                                                                              <  2:256436/150‑1 (MQ=18)
                        ttGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTg                                                                                                                             <  1:4784514/150‑1 (MQ=255)
                                       aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGatca                                                                                                              >  2:1609069/1‑150 (MQ=255)
                                          ggaggaTAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTACAAAC‑AAG‑CACGTCAGTACGGGGACCGGGGGTGGAGCAAcc                                                                                                           >  1:2794703/1‑150 (MQ=14)
                                             gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                                        >  1:2719180/1‑150 (MQ=255)
                                              tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                                       <  1:1494815/150‑1 (MQ=255)
                                              tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                                       <  1:3895704/150‑1 (MQ=255)
                                               ataGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                                      <  2:3282341/150‑1 (MQ=255)
                                                taGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                                                     <  1:3282871/150‑1 (MQ=255)
                                                      gcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                                               <  2:473895/150‑1 (MQ=255)
                                                       cATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAGTTTCAGG‑AAG‑AATGGCATTACGGTGTAGAGGGGTGGAGCAGACAGGTAGCtctacg                                                                                              >  2:4005284/1‑146 (MQ=14)
                                                        aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTTTTTG‑TAGTTAATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg                                                                                             >  1:2140821/1‑150 (MQ=11)
                                                        aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCACATAACCCGAAGGTCGTAGGTGCAAATCCCGACCCCGCAACCACTAAAAAGTTGATG‑AAG‑GAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTCGCTCGTCgg                                                                                             >  1:3940570/1‑150 (MQ=14)
                                                            aCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTACCGATCCGGCCCCCGAAACCAATACAATCAACCTG‑GCG‑TAAAGCAGTACGGTCACGCGGGGGGGAGCAGCCTGGTAGCTCGGCGGGcgc                                                                                         >  1:3173979/1‑148 (MQ=11)
                                                             cGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                                                        <  2:4326203/150‑1 (MQ=255)
                                                             cGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGAAGCCTGGTAGCTCGTCGGGCTCa                                                                                        >  1:253061/1‑150 (MQ=255)
                                                               gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATa                                                                                      <  2:558769/150‑1 (MQ=40)
                                                                 cgcgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAc                                                                                    <  2:1647204/150‑1 (MQ=40)
                                                                    ggggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGTTG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCg                                                                                 >  2:2008792/1‑150 (MQ=255)
                                                                       gTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                                              <  1:2633829/150‑1 (MQ=40)
                                                                            gcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgt                                                                         <  1:1950489/150‑1 (MQ=40)
                                                                             cagcCTTGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgtc                                                                        <  1:1899190/150‑1 (MQ=255)
                                                                             cagcCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgtc                                                                        <  1:2434603/150‑1 (MQ=255)
                                                                                  tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTc                                                                   <  2:680382/150‑1 (MQ=40)
                                                                                  tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTc                                                                   <  1:3216156/150‑1 (MQ=40)
                                                                                   ggTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCa                                                                  <  2:3060197/150‑1 (MQ=40)
                                                                                   ggTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCa                                                                  <  1:2263380/150‑1 (MQ=40)
                                                                                       gCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATc                                                              <  1:3450259/150‑1 (MQ=40)
                                                                                        cTCGTCTGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCACTTATATTTTGCTG‑AAG‑TAACGCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATcc                                                             <  2:2047328/150‑1 (MQ=11)
                                                                                         tcgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCg                                                            <  1:575822/150‑1 (MQ=40)
                                                                                          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TACAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCgg                                                           >  2:888032/1‑150 (MQ=255)
                                                                                             cgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccc                                                        <  2:1971845/150‑1 (MQ=40)
                                                                                              gggCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGcccc                                                       <  1:3396707/150‑1 (MQ=40)
                                                                                                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCg                                                     <  2:3749671/150‑1 (MQ=40)
                                                                                                     tAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGTTG‑AAG‑TCAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAcc                                                <  2:1033455/150‑1 (MQ=255)
                                                                                                      aaCCCGAAGGTCGTCGGTTCATATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCa                                               <  2:767075/150‑1 (MQ=255)
                                                                                                      aaCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCa                                               <  1:2599502/150‑1 (MQ=40)
                                                                                                      aaCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTAGAGG‑AAG‑TAAAGCATTAATGTGAAGAGGGGTGGAGAAGCATTGAAGCTCGTCGGGCGCATAACACGAAGGACGACGGATCAAATCCGCCCCCCGCAACCa                                               >  1:3135094/1‑150 (MQ=14)
                                                                                                        cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt                                             >  1:503423/1‑150 (MQ=40)
                                                                                                         ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATc                                            >  1:2787378/1‑150 (MQ=255)
                                                                                                          cGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATCa                                           >  2:241246/1‑150 (MQ=255)
                                                                                                          cGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑AAATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACACGAAGGGCGACGGTTCAAATCCGGCCCCCGCACCCCATCa                                           >  2:3441405/1‑150 (MQ=14)
                                                                                                           gAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCCGCCCCCGCAACCaaaaaa                                          >  2:4547035/1‑146 (MQ=255)
                                                                                                              ggTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtt                                       <  2:3330708/150‑1 (MQ=255)
                                                                                                               gtcgtcGGTTCAAATCCGTCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                      <  1:4515512/150‑1 (MQ=255)
                                                                                                                   tcgGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatg                                  <  1:383451/150‑1 (MQ=255)
                                                                                                                     ggTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGGGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                <  1:1173764/150‑1 (MQ=255)
                                                                                                                     ggTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                >  1:513682/1‑150 (MQ=255)
                                                                                                                     ggTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                >  1:2160860/1‑150 (MQ=255)
                                                                                                                         cAAATCCGGCCCCCGCAACCAATTAAATTTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                            <  2:1220629/150‑1 (MQ=255)
                                                                                                                         cAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                            <  2:1220460/150‑1 (MQ=255)
                                                                                                                         cAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCGGCCCCCGCAACCAATCAAATAGGAGGAAGTac                            >  2:3776740/1‑149 (MQ=255)
                                                                                                                          ataTCCTGACCCCTCAAACAATTATAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                           <  1:1004472/148‑1 (MQ=17)
                                                                                                                            aTCCTGCACCCGCAACCAATTAAATTGTGTTG‑TAG‑TAAAGGAGTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAg                         <  2:4602248/150‑1 (MQ=17)
                                                                                                                            aTCCTGACCCCGCAACCAATTATATTTGGTTG‑TAT‑TCTCGCATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAg                         <  2:2152718/150‑1 (MQ=17)
                                                                                                                            aTCCGTCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGGAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAg                         <  2:2152717/150‑1 (MQ=255)
                                                                                                                            aTCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAg                         <  1:2008792/150‑1 (MQ=255)
                                                                                                                             tCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGc                        <  2:513563/150‑1 (MQ=255)
                                                                                                                                ttaccccGCAACCAATTATATTTTGTTG‑AAG‑TATTGGAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                     <  1:144326/147‑1 (MQ=18)
                                                                                                                                ggCCCCCGCAACCAATTATACTTTGATG‑AAG‑TAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                     <  2:2871487/150‑1 (MQ=255)
                                                                                                                                ggCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCACCCGCAACCAATCAAAAATGATGAAGAAAAAGCAGt                     >  1:4236835/1‑150 (MQ=255)
                                                                                                                                ggCCCCCGCAACCAAATCTATATTGTTG‑AAG‑TCATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                     <  2:2871537/150‑1 (MQ=17)
                                                                                                                                 gCCCCCTCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                    <  2:3091980/150‑1 (MQ=255)
                                                                                                                                 gCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                    >  1:141943/1‑150 (MQ=255)
                                                                                                                                  cccccGCAACCAATTATAATTTGATG‑AAG‑TAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                   <  2:4745828/150‑1 (MQ=255)
                                                                                                                                    cccGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTCATAGCAGTAAgg                 >  2:4457676/1‑150 (MQ=255)
                                                                                                                                    cccGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg                 <  1:2317273/150‑1 (MQ=255)
                                                                                                                                     ccGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAATAGCAGTACGGt                >  1:3920707/1‑150 (MQ=255)
                                                                                                                                        cAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc             >  1:458712/1‑150 (MQ=255)
                                                                                                                                         aaCCAATTAAATTTTGATG‑TAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg            <  2:4870886/150‑1 (MQ=255)
                                                                                                                                         aaCCAATTAAATTTTGATG‑AAG‑TTTTGGAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg            <  2:4871958/150‑1 (MQ=18)
                                                                                                                                                   aaTTTTATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGAAGAAGTAAAAGCAGTAAGGTGACGCGGGGTGGag  >  2:1065677/1‑150 (MQ=255)
                                                                                                                                                   aaTTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGag  >  1:1785696/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATG‑AAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAG  >  CP009273/2940681‑2940977

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