Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,006,260 G→A 18.9% intergenic (‑84/‑153) metR ← / → metE methionine biosynthesis regulon transcriptional regulator/5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,006,2600GA18.9% 621.8 / 126.3 249intergenic (‑84/‑153)metR/metEmethionine biosynthesis regulon transcriptional regulator/5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (98/104);  new base A (15/32);  total (113/136)
Fisher's exact test for biased strand distribution p-value = 5.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

AGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATAATTAGAGGAAGAA  >  CP009273/4006111‑4006409
                                                                                                                                                     |                                                                                                                                                     
aGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACg                                                                                                                                                       <  2:273421/150‑1 (MQ=255)
  cGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGc                                                                                                                                                     <  2:3227470/150‑1 (MQ=255)
   ggCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCGTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCt                                                                                                                                                    <  1:4167735/150‑1 (MQ=255)
    gCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTa                                                                                                                                                   >  1:1690968/1‑150 (MQ=255)
         gagcgagcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATc                                                                                                                                              <  1:2601695/150‑1 (MQ=255)
          agcgagcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATcc                                                                                                                                             >  1:480903/1‑150 (MQ=255)
           gcgagcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCt                                                                                                                                            <  1:921460/150‑1 (MQ=255)
           gcgagcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCt                                                                                                                                            <  1:1594655/150‑1 (MQ=255)
            cgagcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCtt                                                                                                                                           <  1:304382/150‑1 (MQ=255)
             gagcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTc                                                                                                                                          <  1:521092/150‑1 (MQ=255)
              agcCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACGTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCg                                                                                                                                         <  2:4677665/150‑1 (MQ=255)
                    aGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGGGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCtttt                                                                                                                                   >  2:1504289/1‑150 (MQ=255)
                    aGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCtttt                                                                                                                                   <  1:1344133/150‑1 (MQ=255)
                     gTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTa                                                                                                                                  >  2:3136375/1‑150 (MQ=255)
                     gTTCCGCAAAGCTTGTAGCGTTTTCAGGTTTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGTATATTCAGTAACCGCCAATTATGGATGTGTAAACATCTTGACAGCTAAAATCCTTCGTCTTTTa                                                                                                                                  <  2:4465300/150‑1 (MQ=255)
                      ttCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTaa                                                                                                                                 <  2:1825518/150‑1 (MQ=255)
                         cGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAAtt                                                                                                                              <  1:1744698/150‑1 (MQ=255)
                         cGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAAtt                                                                                                                              <  2:1872305/150‑1 (MQ=255)
                          gCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAAttt                                                                                                                             >  1:1168341/1‑150 (MQ=255)
                           cAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTa                                                                                                                            <  2:1360938/150‑1 (MQ=255)
                             aCGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATg                                                                                                                          <  2:3980888/150‑1 (MQ=255)
                                   gTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGt                                                                                                                    >  2:1596535/1‑150 (MQ=255)
                                     aGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTg                                                                                                                  <  2:3381002/150‑1 (MQ=255)
                                     aGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTg                                                                                                                  >  2:1240244/1‑150 (MQ=255)
                                        gTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCt                                                                                                               >  2:2860947/1‑150 (MQ=255)
                                          tttCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGc                                                                                                             >  2:901569/1‑150 (MQ=255)
                                          tttCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGc                                                                                                             <  1:1855360/150‑1 (MQ=255)
                                             cAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGtt                                                                                                          <  2:4076830/150‑1 (MQ=255)
                                               ggTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAAttt                                                                                                                             >  1:82060/1‑129 (MQ=255)
                                               ggTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAAttt                                                                                                                             <  2:82060/129‑1 (MQ=255)
                                                gtgtTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTtct                                                                                                       >  2:8138/1‑150 (MQ=255)
                                                 tgtTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTtctc                                                                                                      >  2:4171507/1‑150 (MQ=255)
                                                  gtcttACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTcc                                                                                                     >  2:1936582/4‑150 (MQ=255)
                                                     ttACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCAcc                                                                                                  >  2:1654008/1‑150 (MQ=255)
                                                       aCTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCAcccc                                                                                                >  2:1079007/1‑150 (MQ=255)
                                                          tCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGt                                                                                             <  2:3781894/150‑1 (MQ=255)
                                                           cGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTc                                                                                            <  1:3093228/150‑1 (MQ=255)
                                                           cGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTc                                                                                            >  2:3053516/1‑150 (MQ=255)
                                                            gATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCa                                                                                           <  2:1527729/150‑1 (MQ=255)
                                                              tCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCact                                                                                         <  1:1626302/150‑1 (MQ=255)
                                                               cATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCactt                                                                                        <  1:3533951/150‑1 (MQ=255)
                                                                aTGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCactta                                                                                       >  1:3385703/1‑150 (MQ=255)
                                                                aTGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCactta                                                                                       <  2:4215700/150‑1 (MQ=255)
                                                                aTGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCactta                                                                                       <  2:1168312/150‑1 (MQ=255)
                                                                 tGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCacttac                                                                                      >  1:3669597/1‑150 (MQ=255)
                                                                  gAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCacttact                                                                                     >  2:3668613/1‑150 (MQ=255)
                                                                     aGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCa                                                                                  >  1:3435468/1‑150 (MQ=255)
                                                                      gTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAg                                                                                 >  2:2698497/1‑150 (MQ=255)
                                                                             aCTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTc                                                                          <  1:872710/150‑1 (MQ=255)
                                                                                tCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCg                                                                       <  1:2466161/150‑1 (MQ=255)
                                                                                    cATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGa                                                                   >  2:2835667/1‑150 (MQ=255)
                                                                                    cATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGa                                                                   <  1:3134749/150‑1 (MQ=255)
                                                                                    cATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGa                                                                   <  1:3136375/150‑1 (MQ=255)
                                                                                    cATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGa                                                                   >  2:3514992/1‑150 (MQ=255)
                                                                                     aTGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGAt                                                                  >  1:4520657/1‑150 (MQ=255)
                                                                                      tGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATg                                                                 >  2:4407330/1‑150 (MQ=255)
                                                                                       gAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGa                                                                >  1:3340071/1‑150 (MQ=255)
                                                                                         ataatTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAAt                                                              >  2:2530213/1‑150 (MQ=255)
                                                                                         ataatTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAAt                                                              >  2:3360822/1‑150 (MQ=255)
                                                                                         ataatTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAAt                                                              <  2:345871/150‑1 (MQ=255)
                                                                                          taatTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTACCGGGGATGAATa                                                             >  1:2193839/1‑150 (MQ=255)
                                                                                              ttGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACt                                                         >  1:2180585/1‑150 (MQ=255)
                                                                                                 cgcTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTgc                                                      >  2:2508354/1‑150 (MQ=255)
                                                                                                  gcTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTgcc                                                     <  1:1988390/150‑1 (MQ=255)
                                                                                                    ttGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTgccgc                                                   >  2:3086281/1‑150 (MQ=255)
                                                                                                          aaTATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTGTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTccc                                             <  1:4223919/150‑1 (MQ=255)
                                                                                                           atatACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCt                                            <  1:4330030/150‑1 (MQ=255)
                                                                                                             ataCAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCtt                                                <  2:4602326/144‑1 (MQ=255)
                                                                                                             ataCAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCtt                                                >  1:4602326/1‑144 (MQ=255)
                                                                                                               aCAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAAt                                        <  1:3296484/150‑1 (MQ=255)
                                                                                                                cAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAAtt                                       <  1:994052/150‑1 (MQ=255)
                                                                                                                    aaCCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGACGCCTTCCCAAAATTCaaa                                   >  1:2659439/1‑150 (MQ=255)
                                                                                                                     aCCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCaaaa                                  >  2:2095903/1‑150 (MQ=255)
                                                                                                                     aCCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCaaaa                                  <  1:4408409/150‑1 (MQ=255)
                                                                                                                      ccGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAAt                                 >  1:1395367/1‑150 (MQ=255)
                                                                                                                      ccGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAAt                                 >  1:4642102/1‑150 (MQ=255)
                                                                                                                       cGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATc                                <  1:771207/150‑1 (MQ=255)
                                                                                                                       cGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATc                                <  2:3597064/150‑1 (MQ=255)
                                                                                                                       cGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATc                                >  1:2771380/1‑150 (MQ=255)
                                                                                                                         ccAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCa                              >  2:3997540/1‑150 (MQ=255)
                                                                                                                         ccAATTATGGATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCa                              <  1:2242191/150‑1 (MQ=255)
                                                                                                                          cAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCAt                             <  1:4506591/150‑1 (MQ=255)
                                                                                                                           aaTTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATa                            <  1:1194764/150‑1 (MQ=255)
                                                                                                                           aaTTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATa                            <  2:4205590/150‑1 (MQ=255)
                                                                                                                            aTTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAg                           >  1:3872512/1‑150 (MQ=255)
                                                                                                                              tATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGa                         <  1:1012438/150‑1 (MQ=255)
                                                                                                                                tGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGAtt                       >  1:1824770/1‑150 (MQ=255)
                                                                                                                                 ggATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGAttt                      >  2:1815382/1‑150 (MQ=255)
                                                                                                                                 ggATGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGAttt                      <  2:482806/150‑1 (MQ=255)
                                                                                                                                   aTGTGTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTAc                    <  1:2208099/150‑1 (MQ=255)
                                                                                                                                    tgtgTAAACATCTGGACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACa                   <  1:2440195/150‑1 (MQ=255)
                                                                                                                                      tgtAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACata                 <  2:2319277/150‑1 (MQ=255)
                                                                                                                                      tgtAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTAACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACata                 <  2:3933541/150‑1 (MQ=255)
                                                                                                                                      tgtAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGCCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACata                 <  2:480903/150‑1 (MQ=255)
                                                                                                                                         aaaCATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATaa              >  1:2102506/1‑150 (MQ=255)
                                                                                                                                             aTCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAATATCCATAGGATTTACATATAATTag          <  1:1056281/150‑1 (MQ=255)
                                                                                                                                             aTCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATCGGATTTACAGATAATTag          >  2:4575675/1‑150 (MQ=255)
                                                                                                                                             aTCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATAATTag          >  2:3722872/1‑150 (MQ=255)
                                                                                                                                               cTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATAATTagag        <  2:4751804/150‑1 (MQ=255)
                                                                                                                                                 ggACAGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGACTAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACCTATAATTAGAGga      >  1:3718711/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATAATTAGAGGAAGAA  >  CP009273/4006111‑4006409

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: