Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,596,487 T→A 33.7% *344L (TAA→TTA)  rsmC ← 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,596,4870TA33.7% 265.6 / 226.6 217*344L (TAA→TTA) rsmC16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (69/75);  new base A (26/47);  total (95/122)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 220 total aligned reads displayed.

GCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGC  >  CP009273/4596342‑4596634
                                                                                                                                                 |                                                                                                                                                   
gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTaa                                                                                                                                                   <  1:1468047/148‑1 (MQ=255)
 cGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAATTCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAAc                                                                                                                                                  >  1:2073189/1‑148 (MQ=255)
   tGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCt                                                                                                                                                >  2:629388/1‑148 (MQ=255)
    gCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCtt                                                                                                                                               <  2:3336366/148‑1 (MQ=255)
    gCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCtt                                                                                                                                               <  2:1900835/148‑1 (MQ=255)
    gCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCtt                                                                                                                                               <  1:481136/148‑1 (MQ=255)
    gCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCtt                                                                                                                                               >  2:1839244/1‑148 (MQ=255)
       ggTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttct                                                                                                                                            >  1:3441267/1‑148 (MQ=255)
           tGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGc                                                                                                                                        <  1:1794560/148‑1 (MQ=255)
            ggTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGcc                                                                                                                                       <  2:2406804/148‑1 (MQ=255)
             gTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCt                                                                                                                                      <  1:332183/148‑1 (MQ=255)
             gTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCt                                                                                                                                      <  1:2434774/148‑1 (MQ=255)
             gTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTGATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCt                                                                                                                                      <  2:1937791/148‑1 (MQ=255)
                tgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGc                                                                                                                                   <  2:2979495/148‑1 (MQ=255)
                     gcGGGGGCTGATTCCAGAGGCGCCGCTGGCGCGTCTTATCTGGCCTACGACGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTc                                                                                                                              <  1:2732854/148‑1 (MQ=39)
                     gcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTc                                                                                                                              <  2:2398637/148‑1 (MQ=255)
                        ggggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATa                                                                                                                           <  1:777140/148‑1 (MQ=255)
                        ggggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATa                                                                                                                           <  1:2564452/148‑1 (MQ=255)
                         gggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATaa                                                                                                                          <  2:1211474/148‑1 (MQ=255)
                          ggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAAt                                                                                                                         <  1:1332171/148‑1 (MQ=255)
                          ggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAAt                                                                                                                         <  1:615184/148‑1 (MQ=255)
                            cagatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATgg                                                                                                                       <  2:2627917/148‑1 (MQ=255)
                               atgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                                                    >  1:977063/1‑148 (MQ=255)
                               atgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                                                    >  2:918642/1‑148 (MQ=255)
                                tgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGcgcg                                                                                                                   <  1:2841978/148‑1 (MQ=255)
                                 gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGa                                                                                                                  >  1:3369458/1‑148 (MQ=255)
                                 gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGa                                                                                                                  <  2:121039/148‑1 (MQ=255)
                                    agatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATac                                                                                                               >  1:763252/1‑148 (MQ=255)
                                     gatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATaca                                                                                                              >  1:1738149/1‑148 (MQ=255)
                                      atgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATacac                                                                                                             <  2:3389489/148‑1 (MQ=255)
                                      atgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATacac                                                                                                             <  2:2832635/148‑1 (MQ=255)
                                       tgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACAcc                                                                                                            <  2:1548872/148‑1 (MQ=255)
                                       tgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACAcc                                                                                                            <  2:3155149/148‑1 (MQ=255)
                                            cGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGaa                                                                                                       >  2:1312643/1‑148 (MQ=255)
                                            cGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGaa                                                                                                       >  1:2502662/1‑148 (MQ=255)
                                              cTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAgc                                                                                                     >  1:393728/1‑148 (MQ=255)
                                              cTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAgc                                                                                                     >  1:2184198/1‑148 (MQ=255)
                                                  cgcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGccc                                                                                                 <  2:3043014/148‑1 (MQ=255)
                                                  cgcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGccc                                                                                                 <  2:1520846/148‑1 (MQ=255)
                                                    cgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTg                                                                                               <  1:3521715/148‑1 (MQ=255)
                                                        ttATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTg                                                                                           >  2:2029213/1‑148 (MQ=255)
                                                              ggCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGAt                                                                                     >  1:2493837/1‑148 (MQ=255)
                                                               gCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATc                                                                                    >  1:1945153/1‑148 (MQ=255)
                                                                   aCGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACtt                                                                                >  2:823496/1‑148 (MQ=255)
                                                                      aaGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGt                                                                             >  1:821834/1‑148 (MQ=255)
                                                                      aaGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTAGCGCGATCACTTCGt                                                                             >  2:68272/1‑148 (MQ=255)
                                                                       aGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTg                                                                            >  2:1769209/1‑148 (MQ=255)
                                                                        gggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTgg                                                                           >  2:772988/1‑148 (MQ=255)
                                                                        gggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTgg                                                                           >  1:3192530/1‑148 (MQ=255)
                                                                         ggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGa                                                                          >  1:3048644/1‑148 (MQ=255)
                                                                         ggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGa                                                                          <  1:1624521/148‑1 (MQ=255)
                                                                         ggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGa                                                                          >  1:1537294/1‑148 (MQ=255)
                                                                           cTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAg                                                                        >  1:3395724/1‑148 (MQ=255)
                                                                            tAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGc                                                                       <  2:3550369/148‑1 (MQ=255)
                                                                             aaCGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGcc                                                                      <  1:863049/148‑1 (MQ=255)
                                                                             aaCGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGcc                                                                      >  1:3004320/1‑148 (MQ=255)
                                                                               cGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCaa                                                                    <  1:1039600/148‑1 (MQ=255)
                                                                                gTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCaaa                                                                   <  2:2276010/148‑1 (MQ=255)
                                                                                  gCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATg                                                                 >  2:1787134/1‑148 (MQ=255)
                                                                                   cAGGTTTTGTAGGTCGGATAATGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGt                                                                <  2:2931078/148‑1 (MQ=255)
                                                                                   cAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGt                                                                <  2:3303860/148‑1 (MQ=255)
                                                                                     ggTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtct                                                              <  1:1771188/148‑1 (MQ=255)
                                                                                      gTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtctc                                                             >  2:2353910/1‑148 (MQ=255)
                                                                                       ttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCa                                                            >  1:3235512/1‑148 (MQ=255)
                                                                                       ttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCa                                                            >  1:1109510/1‑148 (MQ=255)
                                                                                        tttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCAt                                                           <  1:2896806/148‑1 (MQ=255)
                                                                                            tAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAg                                                       >  1:2112300/1‑148 (MQ=255)
                                                                                               gTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCAc                                                    >  1:6739/1‑148 (MQ=255)
                                                                                                  ggATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTc                                                 >  2:776468/1‑148 (MQ=255)
                                                                                                   gATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTCCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTcc                                                >  2:2197211/1‑148 (MQ=255)
                                                                                                     tAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCgg                                              <  2:417541/148‑1 (MQ=255)
                                                                                                      aaGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGa                                             <  1:2596548/148‑1 (MQ=255)
                                                                                                         gCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATaa                                          <  1:1391281/148‑1 (MQ=255)
                                                                                                         gCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATaa                                          >  2:572446/1‑148 (MQ=255)
                                                                                                            ttCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGc                                       <  2:1168977/148‑1 (MQ=255)
                                                                                                              cACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAg                                     <  2:1952622/148‑1 (MQ=255)
                                                                                                                 gccgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGaa                                  >  1:1840209/1‑148 (MQ=255)
                                                                                                                 gccgcATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGaa                                  <  2:551660/148‑1 (MQ=255)
                                                                                                                   cgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAgg                                <  2:1243995/148‑1 (MQ=255)
                                                                                                                    gcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGc                               <  2:98137/148‑1 (MQ=255)
                                                                                                                    gcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGc                               <  2:3515859/148‑1 (MQ=255)
                                                                                                                     cATCCGACACGGTATTCGGCTAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCg                              <  1:3341855/148‑1 (MQ=255)
                                                                                                                     cATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCg                              <  2:2754000/148‑1 (MQ=255)
                                                                                                                        ccGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTc                           >  1:828818/1‑148 (MQ=255)
                                                                                                                           acacGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCt                        >  1:1301764/1‑148 (MQ=255)
                                                                                                                            cacGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTa                       >  2:290279/1‑148 (MQ=255)
                                                                                                                              cGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACa                     <  2:1895877/148‑1 (MQ=255)
                                                                                                                               ggTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACaa                    >  2:507632/1‑148 (MQ=255)
                                                                                                                                 tATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAAtt                  <  2:3106633/148‑1 (MQ=255)
                                                                                                                                  aTTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTc                 <  1:2340214/148‑1 (MQ=255)
                                                                                                                                  aTTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTc                 <  1:1646041/148‑1 (MQ=255)
                                                                                                                                  aTTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTc                 >  2:2787343/1‑148 (MQ=255)
                                                                                                                                     cGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCTATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCa              >  1:1186525/1‑148 (MQ=255)
                                                                                                                                       gCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGc            <  2:398735/148‑1 (MQ=255)
                                                                                                                                       gCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGc            <  2:2528933/148‑1 (MQ=255)
                                                                                                                                        cGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCt           <  1:581567/148‑1 (MQ=255)
                                                                                                                                           gaTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTcgc        >  2:356817/1‑148 (MQ=255)
                                                                                                                                           gaTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTcgc        >  1:3321059/1‑148 (MQ=255)
                                                                                                                                            aTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTcgcc       >  2:2420193/1‑148 (MQ=255)
                                                                                                                                                 tAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCTCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTcgccgccgc  <  1:3153154/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGC  >  CP009273/4596342‑4596634

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: