Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,576,057 G→A 8.6% V152M (GTG→ATG)  ybhP‑II → phosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,576,0570GA 8.6% 234.1 / 8.2 81V152M (GTG→ATG) ybhP‑IIphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (32/42);  new base A (4/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

GCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGCCTGCA  >  NC_002947/4575914‑4576205
                                                                                                                                               |                                                                                                                                                    
gCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCAc                                                                                                                                                 <  6:325068/149‑1 (MQ=255)
gCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCAc                                                                                                                                                 <  5:28691/149‑1 (MQ=255)
gCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCAc                                                                                                                                                 <  4:147555/149‑1 (MQ=255)
gCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAACGACAGGTGCACGCAGACCGCATGCAc                                                                                                                                                 <  6:313393/149‑1 (MQ=255)
     cacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTc                                                                                                                                            >  8:43539/1‑149 (MQ=255)
       cGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTggcct                                                                                                                                          <  2:42175/149‑5 (MQ=255)
       cGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGt                                                                                                                                          <  6:306549/149‑1 (MQ=255)
          cAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGc                                                                                                                                       <  7:43539/149‑1 (MQ=255)
           aaTGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCt                                                                                                                                      >  2:428690/1‑149 (MQ=255)
                gTTGCTGTCCAAATACCCGATCATCGAGCACGGCAACCTGGACGGGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                                 <  4:383771/149‑1 (MQ=255)
                gTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTg                                                                                                                                     <  1:325054/145‑1 (MQ=255)
                gTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTg                                                                                                                                     >  2:325054/1‑145 (MQ=255)
                gTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                                 <  4:55410/149‑1 (MQ=255)
                    cTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCa                                                                                                                             >  1:149112/1‑149 (MQ=255)
                      gTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                            <  2:186976/148‑1 (MQ=255)
                      gTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                            >  1:186976/1‑148 (MQ=255)
                        ccAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCgg                                                                                                                         >  4:2681/1‑149 (MQ=255)
                         cAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCgg                                                                                                                         >  4:233931/1‑148 (MQ=255)
                         cAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGc                                                                                                                        <  8:383159/149‑1 (MQ=255)
                            aTACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCaaa                                                                                                                     <  8:79925/149‑1 (MQ=255)
                               cccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTCCTGGAAAGCCACCGGCAAAag                                                                                                                   <  8:152763/148‑1 (MQ=255)
                                ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagca                                                                                                                 <  7:144926/149‑1 (MQ=255)
                                ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAag                                                                                                                   >  6:332977/1‑147 (MQ=255)
                                  gATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcagc                                                                                                               >  8:196937/1‑149 (MQ=255)
                                       tcGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAg                                                                                                          >  8:33933/1‑149 (MQ=255)
                                              ccGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctgc                                                                                                   <  2:121770/149‑1 (MQ=255)
                                              ccGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctgc                                                                                                   <  3:2681/149‑1 (MQ=255)
                                                gCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGGCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGt                                                                                                 <  6:343262/149‑1 (MQ=255)
                                                    ccTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAACGACAGGTGCACGCAGACCGCATGCACCTGGCCTTg                                                                                                                                        >  6:361019/1‑106 (MQ=21)
                                                    ccTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAACGACAGGTGCACGCAGACCGCATGCACCTGGCCTTg                                                                                                                                        <  5:361019/106‑1 (MQ=21)
                                                          cGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCATGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTgg                                                                                       <  2:163563/149‑1 (MQ=255)
                                                           gtgtCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGa                                                                                      <  2:164750/149‑1 (MQ=255)
                                                           gtgtCAATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGa                                                                                      >  2:46529/1‑149 (MQ=255)
                                                             gtCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACg                                                                                    >  6:159984/1‑149 (MQ=255)
                                                                gATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGcc                                                                                  <  5:159984/148‑1 (MQ=255)
                                                                 aTCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCATTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCt                                                                                >  5:236631/1‑149 (MQ=255)
                                                                    aCCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGcc                                                                             >  7:82598/1‑149 (MQ=255)
                                                                           ccGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAAc                                                                      >  3:375848/1‑149 (MQ=255)
                                                                               aCGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGccc                                                                  <  7:213130/149‑1 (MQ=255)
                                                                                cgccgcGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGcccc                                                                 >  3:235442/1‑149 (MQ=255)
                                                                                cgccgcGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCATTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGcccc                                                                 <  6:236631/149‑1 (MQ=255)
                                                                                    gcgGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTg                                                             <  5:322999/149‑1 (MQ=255)
                                                                                       gCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatc                                                          <  2:8392/149‑1 (MQ=255)
                                                                                            ttGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGc                                                     >  4:189408/1‑149 (MQ=255)
                                                                                            ttGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGc                                                     <  1:428690/149‑1 (MQ=255)
                                                                                              gCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACAGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCCCCg                                                   >  5:397808/1‑149 (MQ=255)
                                                                                               cACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCgg                                                  <  1:46529/149‑1 (MQ=255)
                                                                                                  tgtgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGa                                               >  5:347748/1‑149 (MQ=255)
                                                                                                    tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                              <  3:189408/148‑1 (MQ=255)
                                                                                                    tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                              <  8:82598/148‑1 (MQ=255)
                                                                                                    tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACt                                             >  1:409375/1‑149 (MQ=255)
                                                                                                    tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACt                                             <  5:332977/149‑1 (MQ=255)
                                                                                                      tgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCTGCCAACGCCCCGGTGATCATCGCCGGCGACtt                                            <  4:168060/148‑1 (MQ=255)
                                                                                                       gCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGcc                                                    >  6:76726/1‑139 (MQ=255)
                                                                                                       gCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGcc                                                    <  5:76726/139‑1 (MQ=255)
                                                                                                        cTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCaa                                         <  7:33933/149‑1 (MQ=255)
                                                                                                           gAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGa                                      >  4:195433/1‑149 (MQ=255)
                                                                                                                gCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTgg                                  >  3:259515/1‑148 (MQ=255)
                                                                                                                 ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGa                                 >  8:306023/1‑148 (MQ=255)
                                                                                                                 ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGa                                 >  6:281113/1‑148 (MQ=255)
                                                                                                                 ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGAACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGAAGCCCTGCAGCCCAACGCCCCGGTTATCATCGCCGGCGACTTCAACGACTGGa                                 >  4:61052/1‑148 (MQ=255)
                                                                                                                   ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                              <  8:49011/149‑1 (MQ=255)
                                                                                                                     ggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAACGc                            >  8:323445/1‑149 (MQ=255)
                                                                                                                     ggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGc                            <  1:149864/149‑1 (MQ=255)
                                                                                                                      gCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCTGCGACTTCAACGACTGGAAAAGc                            >  2:338097/1‑148 (MQ=255)
                                                                                                                      gCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGa                                      >  3:214480/1‑138 (MQ=255)
                                                                                                                      gCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGa                                      <  4:214480/138‑1 (MQ=255)
                                                                                                                         aGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATg                        <  4:375848/149‑1 (MQ=255)
                                                                                                                            acAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                              <  4:45376/123‑1 (MQ=255)
                                                                                                                            acAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                              >  3:45376/2‑124 (MQ=255)
                                                                                                                              aGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAc                   >  3:388431/1‑149 (MQ=255)
                                                                                                                               ggTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCCTCGCCGGCGACTTCCACGACTGGAACAGCCATGGCAccc                  >  8:231934/1‑149 (MQ=255)
                                                                                                                                gTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAACAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCAGGCCAACGCCCCGGTGATCATCGCCGGAGACTTCAAAGACTGGAAAAGCCATGGCacgcg                 >  5:208497/1‑145 (MQ=255)
                                                                                                                                  gCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCa               >  3:32853/1‑149 (MQ=255)
                                                                                                                                    aTGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAcc             <  2:149112/149‑1 (MQ=255)
                                                                                                                                    aTGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAACGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGAAGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAAAGACTGGAAAAGCCATGGCACCCGCAcc             >  1:398856/1‑149 (MQ=255)
                                                                                                                                     tGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAccc            <  4:259515/149‑1 (MQ=255)
                                                                                                                                          tCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGc       <  3:195433/149‑1 (MQ=255)
                                                                                                                                          tCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGc       <  2:409375/149‑1 (MQ=255)
                                                                                                                                              cGTGCACCTGTGGTTGCTGGAAAGCAACGGGCAAATGCAGCTGCTGTTGCTGCGTAAACTGCTGTACGCCCTGCCGGCCAACGCCCGGGTTAGCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGCCTg    <  1:338097/148‑1 (MQ=255)
                                                                                                                                              cGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGCCTg    <  4:32853/148‑1 (MQ=255)
                                                                                                                                               gTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGCCTGCa  >  3:224127/1‑149 (MQ=255)
                                                                                                                                               gTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAACCATGGCAACCGCACCCTTGGCCTGCa  <  4:235442/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGCCTGCA  >  NC_002947/4575914‑4576205

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: ‑