Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,437,156 G→T 17.1% P638P (CCG→CCT mnmC → fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,437,1560GT17.1% 528.6 / 63.0 179P638P (CCG→CCTmnmCfused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base G (63/85);  new base T (11/20);  total (74/105)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

ATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTAC  >  CP009273/2437009‑2437301
                                                                                                                                                   |                                                                                                                                                 
 tatGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGa                                                                                                                                                  <  1:1016135/148‑1 (MQ=255)
  atGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGAt                                                                                                                                                 <  1:599463/148‑1 (MQ=255)
    gCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTc                                                                                                                                               <  1:1770611/148‑1 (MQ=255)
     cgtcgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGGCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTcc                                                                                                                                              <  1:2551004/148‑1 (MQ=255)
     cgtcgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTcc                                                                                                                                              <  1:2551114/148‑1 (MQ=255)
       tcgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCga                                                                                                                                            <  1:2224925/148‑1 (MQ=255)
        cgtTGGCGGAGCAGAAAGATGAGGCGGTACGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgat                                                                                                                                           >  1:2145295/1‑148 (MQ=255)
        cgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgat                                                                                                                                           <  2:2035410/148‑1 (MQ=255)
        cgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgat                                                                                                                                           >  1:3209102/1‑148 (MQ=255)
        cgtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgat                                                                                                                                           <  2:1643659/148‑1 (MQ=255)
         gtTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatg                                                                                                                                          <  2:2084717/148‑1 (MQ=255)
          ttGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatgg                                                                                                                                         <  2:2089027/148‑1 (MQ=255)
           tGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGTGACGAACCGATTCCgatgga                                                                                                                                        <  2:2756984/148‑1 (MQ=255)
           tGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCgatgga                                                                                                                                        <  1:2436324/148‑1 (MQ=255)
           gggCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatgga                                                                                                                                        <  2:1814449/147‑1 (MQ=255)
                   cAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTa                                                                                                                                <  2:1521477/148‑1 (MQ=255)
                     gAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGATCCGATTCCGATGGATGCCAGTACg                                                                                                                              <  1:3251766/148‑1 (MQ=255)
                      aaaGATGAGGCGGTACGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGc                                                                                                                             >  2:3112571/1‑148 (MQ=255)
                      aaaGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGc                                                                                                                             <  2:2694702/148‑1 (MQ=255)
                       aaGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCt                                                                                                                            <  2:3240847/148‑1 (MQ=255)
                        aGATGAGGCGGTACGCGCGCCGGTTTTTGACGATCTCGTTATGTTTGCGGCTTTAGGTTCTCGCGGGTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTg                                                                                                                           >  1:1256576/1‑148 (MQ=255)
                        aGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTg                                                                                                                           >  1:1379273/1‑148 (MQ=255)
                        aGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTg                                                                                                                           <  1:994130/148‑1 (MQ=255)
                        aGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTATGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTg                                                                                                                           >  2:1656122/1‑148 (MQ=255)
                          aTGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTggc                                                                                                                         >  1:1053778/1‑148 (MQ=255)
                             aGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTggcggc                                                                                                                      >  2:2550707/1‑148 (MQ=255)
                               gCGGTAAGCGCGCCGGTTTTTTACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATTAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGt                                                                                                                    >  2:2783101/1‑148 (MQ=255)
                                cGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGtt                                                                                                                   <  1:2790437/148‑1 (MQ=255)
                                 ggTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTa                                                                                                                  <  2:1478310/148‑1 (MQ=255)
                                 ggTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTa                                                                                                                  <  2:1478403/148‑1 (MQ=255)
                                   tAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTaaa                                                                                                                >  2:3119060/1‑148 (MQ=255)
                                   tAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTaaa                                                                                                                >  2:1554572/1‑148 (MQ=255)
                                    aaGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAAc                                                                                                               >  2:1179013/1‑148 (MQ=255)
                                      gcgcgcCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAAccc                                                                                                             >  2:2584565/1‑148 (MQ=255)
                                        gcgcCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGa                                                                                                           <  2:1210614/148‑1 (MQ=255)
                                          gcCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAAt                                                                                                         <  1:642450/148‑1 (MQ=255)
                                          gcCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAAt                                                                                                         <  1:2682713/148‑1 (MQ=255)
                                           ccGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATc                                                                                                        >  2:889577/1‑148 (MQ=255)
                                           ccGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATc                                                                                                        >  1:1303996/1‑148 (MQ=255)
                                           ccGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATc                                                                                                        >  1:1771198/1‑148 (MQ=255)
                                              gTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGt                                                                                                     <  2:107002/148‑1 (MQ=255)
                                               tttttGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGtt                                                                                                    <  2:1345736/148‑1 (MQ=255)
                                                  ttGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATg                                                                                                 <  2:1348462/148‑1 (MQ=255)
                                                  ttGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATg                                                                                                 <  1:1081052/148‑1 (MQ=255)
                                                    gACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATggg                                                                                               <  2:1446271/148‑1 (MQ=255)
                                                        aTCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGGTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTAGGGGTGCg                                                                                           >  1:3169565/1‑148 (MQ=255)
                                                            ctTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTAGGGGTGCGGaaa                                                                                       >  2:2808254/1‑148 (MQ=255)
                                                             tttATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAAt                                                                                      <  1:2720790/148‑1 (MQ=255)
                                                             tttATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAAt                                                                                      >  1:3027462/1‑148 (MQ=255)
                                                                aTGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAAttgt                                                                                   >  1:1973226/1‑148 (MQ=255)
                                                                aTGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAAttgt                                                                                   <  2:1863308/148‑1 (MQ=255)
                                                                aTGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGGTATGGGTGCGGAAAttgt                                                                                   >  2:2688140/1‑148 (MQ=255)
                                                                    ttGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGaa                                                                               >  1:372241/1‑148 (MQ=255)
                                                                    ttGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGaa                                                                               >  1:30307/1‑148 (MQ=255)
                                                                    ttGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGaa                                                                               >  2:2811414/1‑148 (MQ=255)
                                                                       cGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCTAATCGGTTATGGGTGCGGAAATTGTTGAAggg                                                                            <  2:193931/148‑1 (MQ=255)
                                                                          cTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTaa                                                                         <  1:1762451/148‑1 (MQ=255)
                                                                           tttAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTaaa                                                                        <  2:189074/148‑1 (MQ=255)
                                                                            ttAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAg                                                                       <  2:96866/148‑1 (MQ=255)
                                                                            ttAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAg                                                                       <  1:3148935/148‑1 (MQ=255)
                                                                             tAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGc                                                                      <  2:526429/148‑1 (MQ=255)
                                                                                gTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGTAGGGTAAAGCGGt                                                                   >  1:1403947/1‑148 (MQ=255)
                                                                                gTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGt                                                                   <  2:898751/148‑1 (MQ=255)
                                                                                 ttCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGtt                                                                  >  1:550017/1‑148 (MQ=255)
                                                                                   ctcGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTaa                                                                <  1:1462897/148‑1 (MQ=255)
                                                                                       cgGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCg                                                            >  1:3240093/1‑148 (MQ=255)
                                                                                       cgGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCg                                                            >  1:2813649/1‑148 (MQ=255)
                                                                                        ggTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCgg                                                           <  1:3122854/148‑1 (MQ=255)
                                                                                        ggTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCgg                                                           <  2:2231541/148‑1 (MQ=255)
                                                                                            tgtgTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTa                                                       <  1:659192/148‑1 (MQ=255)
                                                                                            tgtgTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTa                                                       <  1:2319791/148‑1 (MQ=255)
                                                                                                 tCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTg                                                  >  2:2087539/1‑148 (MQ=255)
                                                                                                 tCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTg                                                  <  2:2900200/148‑1 (MQ=255)
                                                                                                  cTGCCCCGCGGCGTGCCGAGATTCTGGCGGCGCAGCTGAGCGACGAACCTATTCCGCTGGATGCCAGTACGCTGGCGGCGTTACACCCGACTCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGc                                                 <  2:1775992/148‑1 (MQ=255)
                                                                                                        cGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGc                                           >  1:936036/1‑148 (MQ=255)
                                                                                                             tgtgCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAgtgt                                      >  1:1171927/1‑148 (MQ=255)
                                                                                                                  cGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGcc                                 <  2:2522856/148‑1 (MQ=255)
                                                                                                                    agaTTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCgg                               >  1:1905058/1‑148 (MQ=255)
                                                                                                                     gaTTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGa                              <  1:854571/148‑1 (MQ=255)
                                                                                                                     gaTTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGa                              >  2:2375100/1‑148 (MQ=255)
                                                                                                                     gaTTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGa                              <  1:3234080/148‑1 (MQ=255)
                                                                                                                      aTTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGAt                             >  1:739891/1‑148 (MQ=255)
                                                                                                                        tCTGGCGGCGCAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATgc                           <  1:557027/148‑1 (MQ=255)
                                                                                                                          tGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATgcgg                         <  2:1192319/148‑1 (MQ=255)
                                                                                                                           ggcggcGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGGTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATgcggc                        >  1:847628/1‑148 (MQ=255)
                                                                                                                                cgcAGATGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGaa                   >  2:443015/1‑148 (MQ=255)
                                                                                                                                    gATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGcc               <  1:2704523/148‑1 (MQ=255)
                                                                                                                                    gATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGcc               <  1:2395689/148‑1 (MQ=255)
                                                                                                                                    gATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGcc               <  1:2704662/148‑1 (MQ=255)
                                                                                                                                    gATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGcc               <  1:1554574/148‑1 (MQ=255)
                                                                                                                                     aTGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCt              <  1:2517380/148‑1 (MQ=255)
                                                                                                                                     aTGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCt              <  1:2912938/148‑1 (MQ=255)
                                                                                                                                      tGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCtt             <  1:2759733/148‑1 (MQ=255)
                                                                                                                                      tGAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCtt             <  2:36018/148‑1 (MQ=255)
                                                                                                                                      tGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCtt             <  2:3209102/148‑1 (MQ=255)
                                                                                                                                       gAGCGACGAACCTATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTa            <  2:2991951/148‑1 (MQ=255)
                                                                                                                                        aGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTAt           <  2:1997186/148‑1 (MQ=255)
                                                                                                                                                aCCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTa   >  1:3304448/1‑148 (MQ=255)
                                                                                                                                                aCCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTa   <  2:1598697/148‑1 (MQ=255)
                                                                                                                                                 ccGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGGTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTAc  >  1:3061173/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTAC  >  CP009273/2437009‑2437301

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: