Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,948,689 G→A 24.2% intergenic (+19/+146) cmoB → / ← torZ tRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent/trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,948,6890GA24.2% 538.1 / 182.5 281intergenic (+19/+146)cmoB/torZtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (106/107);  new base A (28/40);  total (134/147)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

AAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAC  >  CP009273/1948540‑1948838
                                                                                                                                                     |                                                                                                                                                     
aagaGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCg                                                                                                                                                       <  1:3601367/150‑1 (MQ=255)
  gagCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGa                                                                                                                                                     <  1:2502109/150‑1 (MQ=255)
    gcagcGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATg                                                                                                                                                   <  1:5179464/150‑1 (MQ=255)
    gcagcGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATg                                                                                                                                                   <  1:3773980/150‑1 (MQ=255)
      agcGACGCACCGACTGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATgcg                                                                                                                                                 <  1:1112013/150‑1 (MQ=255)
      agcGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATgcg                                                                                                                                                 <  2:1428800/150‑1 (MQ=255)
      agcGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATGTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATgcg                                                                                                                                                 <  2:70157/150‑1 (MQ=255)
            gCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGa                                                                                                                                           >  2:1595039/1‑150 (MQ=255)
            gCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGa                                                                                                                                           >  1:1489359/1‑150 (MQ=255)
              aCCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAAc                                                                                                                                         <  2:987525/150‑1 (MQ=255)
                cGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGc                                                                                                                                       >  2:4810451/1‑150 (MQ=255)
                  aaTGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCt                                                                                                                                     >  2:2323332/1‑150 (MQ=255)
                  aaTGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCt                                                                                                                                     >  2:4246288/1‑150 (MQ=255)
                   atggatggTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCtt                                                                                                                                    >  2:2003438/1‑150 (MQ=255)
                    tggatggTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                   >  1:4847251/1‑150 (MQ=255)
                     ggatggTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                  <  2:2145288/150‑1 (MQ=255)
                      gatggTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATc                                                                                                                                 <  1:1867417/150‑1 (MQ=255)
                      gatggTGACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATc                                                                                                                                 <  2:4740649/150‑1 (MQ=255)
                             aCCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                          <  1:1329929/150‑1 (MQ=255)
                              ccGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                         >  2:3236453/1‑150 (MQ=255)
                               cGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                        <  2:151010/150‑1 (MQ=255)
                                    cGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAcc                                                                                                                   <  1:405892/150‑1 (MQ=255)
                                         gCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCt                                                                                                              >  2:4882412/1‑150 (MQ=255)
                                         gCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCt                                                                                                              <  1:2230/150‑1 (MQ=255)
                                            gATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCGAAt                                                                                                           >  2:3856287/1‑150 (MQ=255)
                                            gATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAt                                                                                                           >  1:2375635/1‑150 (MQ=255)
                                             aTTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                                                                          >  1:3290131/1‑150 (MQ=255)
                                             aTTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                                                                          <  1:3533677/150‑1 (MQ=255)
                                               ttCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                                                                        <  1:1613906/150‑1 (MQ=255)
                                               ttCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                                                                        <  2:3719702/150‑1 (MQ=255)
                                                 ccTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAat                                                                                                      <  2:1317936/150‑1 (MQ=255)
                                                   tCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTCAGTCAatat                                                                                                    >  2:4513991/1‑150 (MQ=255)
                                                   tCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatat                                                                                                    >  2:2601442/1‑150 (MQ=255)
                                                   tCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatat                                                                                                    >  2:2478119/1‑150 (MQ=255)
                                                    cGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatata                                                                                                   <  2:4481651/150‑1 (MQ=255)
                                                    cGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatata                                                                                                   <  2:2826206/150‑1 (MQ=255)
                                                       cccACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                                                                <  1:4319000/150‑1 (MQ=255)
                                                       cccACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                                                                <  2:3801876/150‑1 (MQ=255)
                                                          acaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAg                                                                                             >  1:4538731/1‑150 (MQ=255)
                                                          acaTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAg                                                                                             >  1:106318/1‑150 (MQ=255)
                                                            aTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAggg                                                                                           >  2:2910127/1‑150 (MQ=255)
                                                            aTGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAggg                                                                                           <  2:1598156/150‑1 (MQ=255)
                                                               aTCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACt                                                                                        >  1:1423444/1‑150 (MQ=255)
                                                                 ccGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAt                                                                                      <  1:203861/150‑1 (MQ=255)
                                                                  cGGGTACAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                                     <  2:5051778/150‑1 (MQ=255)
                                                                  cGGGTAAAACGGTGGAAGGTTATCCGGCACCTAACCGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                                     <  2:2532051/150‑1 (MQ=255)
                                                                  cGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                                     <  2:5110598/150‑1 (MQ=255)
                                                                     gTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAg                                                                                  <  2:1722968/150‑1 (MQ=255)
                                                                      tAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAgg                                                                                 >  1:958859/1‑150 (MQ=255)
                                                                        aaaCGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTCTTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                                               >  2:2239735/1‑150 (MQ=255)
                                                                        aaaCGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                                               <  1:3085281/150‑1 (MQ=255)
                                                                         aaCGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCt                                                                              >  1:4515918/1‑150 (MQ=255)
                                                                         aaCGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCt                                                                              >  1:2481616/1‑150 (MQ=255)
                                                                          aCGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTg                                                                             <  2:4680509/150‑1 (MQ=255)
                                                                           cGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGa                                                                            >  2:1880748/1‑150 (MQ=255)
                                                                           cGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGa                                                                            >  2:1904635/1‑150 (MQ=255)
                                                                              tGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATaa                                                                         >  1:3431599/1‑150 (MQ=255)
                                                                              tGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATaa                                                                         <  1:2281772/150‑1 (MQ=255)
                                                                              tGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATaa                                                                         >  1:823999/1‑150 (MQ=255)
                                                                               ggAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                                                                        >  1:4296534/1‑150 (MQ=255)
                                                                                 aaGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCa                                                                      >  1:3319161/1‑150 (MQ=255)
                                                                                  aGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCAt                                                                     >  2:2705581/1‑150 (MQ=255)
                                                                                  aGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCAt                                                                     >  2:2705615/1‑150 (MQ=255)
                                                                                      tATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                                                 <  1:768887/150‑1 (MQ=255)
                                                                                       aTCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcgc                                                                >  2:1118747/1‑150 (MQ=255)
                                                                                        tCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCa                                                               <  1:2558671/150‑1 (MQ=255)
                                                                                         ccGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                                              >  1:300951/1‑150 (MQ=255)
                                                                                             cACCTAACCGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAgg                                                          <  1:464677/150‑1 (MQ=255)
                                                                                               ccTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGACGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCa                                                        <  1:890397/150‑1 (MQ=255)
                                                                                                cTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAg                                                           >  1:1737796/1‑146 (MQ=255)
                                                                                                 tAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGc                                                      >  2:2390738/1‑150 (MQ=255)
                                                                                                   aaCGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                                                    >  2:2296561/1‑150 (MQ=255)
                                                                                                         cgGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGtt                                              <  1:3160516/150‑1 (MQ=255)
                                                                                                           gTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTg                                            <  1:4737927/150‑1 (MQ=255)
                                                                                                              cTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                         >  2:655901/1‑150 (MQ=255)
                                                                                                               tGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                                        <  1:1901406/150‑1 (MQ=255)
                                                                                                                  ttGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAcc                                     <  1:12313/150‑1 (MQ=255)
                                                                                                                          aaGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGcc                             >  2:3319134/1‑150 (MQ=255)
                                                                                                                           aGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc                            >  2:4414090/1‑150 (MQ=255)
                                                                                                                           aGCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc                            <  2:1831113/150‑1 (MQ=255)
                                                                                                                            gCCGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCg                           <  2:3011064/150‑1 (MQ=255)
                                                                                                                             ccGTAAAGGTCTGGTATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCATACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGt                          <  1:3519649/150‑1 (MQ=255)
                                                                                                                               gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGttt                        <  1:2971667/150‑1 (MQ=255)
                                                                                                                                 aaaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTcc                      <  2:867856/150‑1 (MQ=255)
                                                                                                                                     gTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGgc                  >  2:2936172/1‑150 (MQ=255)
                                                                                                                                         ggTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTa              <  2:601889/150‑1 (MQ=255)
                                                                                                                                           tATTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAtt            <  2:1248552/150‑1 (MQ=255)
                                                                                                                                            aTTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTc           >  2:3299960/1‑150 (MQ=255)
                                                                                                                                            aTTACTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTc           >  2:5023007/1‑150 (MQ=255)
                                                                                                                                             ttACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCa          <  1:2705581/150‑1 (MQ=255)
                                                                                                                                             ttACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTGTCCGGGCGCTATTCa          >  2:146859/1‑150 (MQ=255)
                                                                                                                                              tACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCAc         <  1:5024725/150‑1 (MQ=255)
                                                                                                                                              tACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCAc         <  2:512637/150‑1 (MQ=255)
                                                                                                                                                cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGt       >  2:3821153/1‑150 (MQ=255)
                                                                                                                                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACAACCTTGCTACGTAAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTGTGCATAACCTCAGCGCCCGGTTCCGGGCGCGAGTCACggc      >  2:4653696/1‑148 (MQ=21)
                                                                                                                                                   ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCtt    <  2:4956068/150‑1 (MQ=255)
                                                                                                                                                   ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCtt    <  2:3086685/150‑1 (MQ=255)
                                                                                                                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAc  >  2:3530712/1‑150 (MQ=255)
                                                                                                                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAc  >  1:3967667/1‑150 (MQ=255)
                                                                                                                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAc  >  2:1916840/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAC  >  CP009273/1948540‑1948838

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: