Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,528,711 T→C 5.3% L351P (CTG→CCG)  ftsA → cell division protein FtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,528,7110TC 5.3% 279.7 / 4.1 94L351P (CTG→CCG) ftsAcell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base T (39/50);  new base C (2/3);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

CGAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGC  >  NC_002947/1528565‑1528856
                                                                                                                                                  |                                                                                                                                                 
cgaGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCTGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTgg                                                                                                                                                 <  6:19650/149‑1 (MQ=255)
cgaGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTgg                                                                                                                                                 <  4:367911/149‑1 (MQ=255)
      ctTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGc                                                                                                                                           >  1:402980/1‑149 (MQ=255)
       ttCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTg                                                                                                                                            <  1:3872/147‑1 (MQ=255)
         cACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACTCCTGGGCGTGc                                                                                                                                           <  4:404023/146‑1 (MQ=255)
         cACCATTATCCAGGCCGAACTGCGCCGCAGTGTCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCg                                                                                                                                         <  3:103317/148‑1 (MQ=255)
            ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                     <  6:399153/149‑1 (MQ=255)
            ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                     >  4:5361/1‑149 (MQ=255)
             cTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAg                                                                                                                                    >  3:360640/1‑149 (MQ=255)
             cTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCggcatg                                                                                                                                    <  4:191570/149‑6 (MQ=255)
              ttATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                   >  4:296093/1‑149 (MQ=255)
               tATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGTAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCCCATGCCGGTACGACTGGGCGTGCCGCACAGCg                                                                                                                                  >  8:101106/1‑149 (MQ=255)
                aTCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGt                                                                                                                                 >  6:103218/1‑149 (MQ=255)
                   cAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTAcc                                                                                                                                          <  1:425592/137‑1 (MQ=255)
                   cAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTAcc                                                                                                                                          >  2:425592/1‑137 (MQ=255)
                    aGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCgg                                                                                                                             <  3:5361/149‑1 (MQ=255)
                      gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCgggg                                                                                                                           <  2:226155/149‑1 (MQ=255)
                      gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCAGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCgggg                                                                                                                           >  6:9658/1‑149 (MQ=255)
                       ccGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCggggg                                                                                                                          <  6:298063/149‑1 (MQ=255)
                       ccGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCggggg                                                                                                                          >  1:356998/1‑149 (MQ=255)
                        cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGc                                                                                                                         <  6:19950/149‑1 (MQ=255)
                         gAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCt                                                                                                                        >  3:118484/1‑149 (MQ=255)
                                   gcAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtg                                                                                                              <  7:65071/149‑1 (MQ=255)
                                    cAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtgg                                                                                                             >  3:282671/1‑149 (MQ=255)
                                     aGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtggt                                                                                                            <  3:296093/149‑1 (MQ=255)
                                      gTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtggtg                                                                                                           <  8:207935/149‑1 (MQ=255)
                                       tGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                          <  5:194650/149‑1 (MQ=255)
                                       tGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                          <  4:387910/149‑1 (MQ=255)
                                       tGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                          <  5:261057/149‑1 (MQ=255)
                                       tGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                          >  8:248258/1‑149 (MQ=255)
                                         gCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgcg                                                                                                         >  6:356908/1‑148 (MQ=255)
                                         gCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgcg                                                                                                         <  7:248258/148‑1 (MQ=255)
                                           tACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGCCGTGGTGCGCaa                                                                                                      >  5:397829/1‑149 (MQ=255)
                                                aGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGCCGCCGTAGAACTGGCACAGGAAATCTTCCACATGCCGGAACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGACACCCg                                                                                                  >  6:343824/1‑148 (MQ=255)
                                                       gTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTAtt                                                                                           <  5:103218/148‑1 (MQ=255)
                                                       gTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGACCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTc                                                                                          <  1:410140/149‑1 (MQ=255)
                                                       gTGCCCGCCGGCATCGGCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTAt                                                                                            <  6:255764/147‑1 (MQ=255)
                                                            cgccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTAGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTg                                                                                     >  1:332869/1‑149 (MQ=255)
                                                                 gCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGtgtg                                                                                >  8:410942/1‑149 (MQ=255)
                                                                  cATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTgg                                                                               <  6:116329/149‑1 (MQ=255)
                                                                   aTCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTTGTGCGCAACCCGATCTATTCCACTg                                                                                     <  3:142630/142‑1 (MQ=255)
                                                                   aTCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTg                                                                                     >  4:142630/1‑142 (MQ=255)
                                                                   aTCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTggg                                                                              >  3:413322/1‑149 (MQ=255)
                                                                       tCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTg                                                                          <  4:239931/149‑1 (MQ=255)
                                                                       tCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTg                                                                          >  3:377776/1‑149 (MQ=255)
                                                                          tCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCt                                                                        >  3:226446/1‑148 (MQ=255)
                                                                                gcggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTAc                                                                 >  8:335670/1‑149 (MQ=255)
                                                                                 gggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACg                                                                <  3:155610/148‑1 (MQ=255)
                                                                                 cggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACg                                                                <  5:356908/149‑1 (MQ=255)
                                                                                 cggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACg                                                                <  2:402980/149‑1 (MQ=255)
                                                                                 cggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACg                                                                <  8:13508/149‑1 (MQ=255)
                                                                                  ggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACgg                                                               <  2:411436/149‑1 (MQ=255)
                                                                                       cgcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgc                                                          <  2:356998/149‑1 (MQ=255)
                                                                                           aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgctg                                                        >  6:71630/1‑147 (MQ=255)
                                                                                            aGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAg                                                     <  1:257525/149‑1 (MQ=255)
                                                                                            aGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAg                                                     <  5:220875/149‑1 (MQ=255)
                                                                                             gATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAg                                                     >  3:60384/1‑148 (MQ=255)
                                                                                              aTGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAg                                                     >  2:63846/1‑147 (MQ=255)
                                                                                               tGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGcc                                                             <  3:262340/138‑1 (MQ=255)
                                                                                               tGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGcc                                                             >  4:262340/1‑138 (MQ=255)
                                                                                                    ggCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                             <  4:226446/149‑1 (MQ=255)
                                                                                                    ggCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                             <  8:218686/149‑1 (MQ=255)
                                                                                                    ggCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                             <  5:71630/149‑1 (MQ=255)
                                                                                                    ggCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACACGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                             <  4:244144/149‑1 (MQ=255)
                                                                                                      cgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAgg                                           <  3:423211/149‑1 (MQ=255)
                                                                                                      cgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAgg                                           <  4:360640/149‑1 (MQ=255)
                                                                                                        ccGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                         <  4:377776/149‑1 (MQ=255)
                                                                                                              aaCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                                   <  7:410942/149‑1 (MQ=255)
                                                                                                              aaCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                                   <  5:9658/149‑1 (MQ=255)
                                                                                                              aaCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                                   <  1:63846/149‑1 (MQ=255)
                                                                                                                   gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGAc                              >  8:359265/1‑149 (MQ=255)
                                                                                                                     cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                            >  7:211291/1‑149 (MQ=255)
                                                                                                                     cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                            >  5:318785/1‑149 (MQ=255)
                                                                                                                     cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                            >  8:62590/1‑149 (MQ=255)
                                                                                                                      gAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGccc                                  >  8:305590/1‑142 (MQ=255)
                                                                                                                      gAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGccc                                  <  7:305590/142‑1 (MQ=255)
                                                                                                                      gAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                           <  2:298467/149‑1 (MQ=255)
                                                                                                                                 ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagca                   <  8:156860/146‑1 (MQ=255)
                                                                                                                                 ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagc                    <  6:56012/145‑1 (MQ=255)
                                                                                                                                 ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagc                    >  5:56012/1‑145 (MQ=255)
                                                                                                                                   gctgtgcGGCACGCCCAGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaacag               >  1:341146/7‑148 (MQ=255)
                                                                                                                                   acaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaacagc              <  7:335670/149‑1 (MQ=255)
                                                                                                                                   acaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaacag               >  5:36827/1‑148 (MQ=255)
                                                                                                                                     aTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAcagcag            <  4:60384/149‑1 (MQ=255)
                                                                                                                                     aTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAcagcag            >  5:35724/1‑149 (MQ=255)
                                                                                                                                       gCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                         >  2:20944/1‑118 (MQ=255)
                                                                                                                                       gCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                         <  1:20944/118‑1 (MQ=255)
                                                                                                                                       gCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAcagcag            >  4:250572/1‑147 (MQ=255)
                                                                                                                                       gCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAcagcag            >  3:105549/1‑147 (MQ=255)
                                                                                                                                        ccGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAcagcagcag         <  4:413322/149‑1 (MQ=255)
                                                                                                                                            tACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAAc     <  4:105549/149‑1 (MQ=255)
                                                                                                                                               gCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGTAGCAACAGc  >  6:77064/1‑149 (MQ=255)
                                                                                                                                               gCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGc  >  1:870/1‑149 (MQ=255)
                                                                                                                                               gCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGc  >  8:354825/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGC  >  NC_002947/1528565‑1528856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: