Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,840 G→A 5.1% P249P (CCG→CCA)  bepA → OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,616,8400GA 5.1% 252.2 / 2.9 78P249P (CCG→CCA) bepAOM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base A (2/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

ATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACCAGCTCACCAGTGATTTGCTGGAT  >  NC_000913/2616692‑2616987
                                                                                                                                                    |                                                                                                                                                   
aTGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTTCAACGCTCTGGATTCGAACCTAAGGCGATTCCAACATTCATCGAAAAATTACTCGATAAGGCGCGTTACTCCTCGCGCCAGCCGGAAATTTTATTGACTCACCCGTTCCCg                                                                                                                                                     >  3:24986/1‑149 (MQ=255)
aTGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCg                                                                                                                                                     <  1:261425/149‑1 (MQ=255)
  gAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCgg                                                                                                                                                    >  7:148467/1‑148 (MQ=255)
   aaCAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGATGCCGGa                                                                                                                                                   >  7:198967/1‑148 (MQ=255)
    aCAGGAATCGGACCGCATTTGTATTCATGTGATGCAACGCTCGGGATTCGATCCGCAGGAGATGCCAACCTTACTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCATCAGTAAATTTTATTGACTCACCCTTTGCCGGaaa                                                                                                                                                 <  8:6291/149‑1 (MQ=255)
    aCAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGaaa                                                                                                                                                 <  5:395261/149‑1 (MQ=255)
               aCCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAACATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCa                                                                                                                                      <  2:275557/149‑1 (MQ=255)
               aCCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCa                                                                                                                                      <  2:269880/149‑1 (MQ=255)
                 cGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGa                                                                                                                                    <  2:46743/149‑1 (MQ=255)
                    aTTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGc                                                                                                                                 >  5:347261/1‑149 (MQ=255)
                    aTTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGc                                                                                                                                 >  7:60173/1‑149 (MQ=255)
                    aTTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGc                                                                                                                                 >  4:319315/1‑149 (MQ=255)
                    aTTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCAGAAAGTCGTCTGGCAGATGc                                                                                                                                 >  8:255789/1‑149 (MQ=255)
                    aTTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCAGAAAGTCGTCTGGCAGATGc                                                                                                                                 >  8:249597/1‑149 (MQ=255)
                       ggTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGccc                                                                                                                               <  8:198967/148‑1 (MQ=255)
                       ggTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGccc                                                                                                                               <  3:319315/148‑1 (MQ=255)
                       ggTATTAAGGTGCTGCAACGCTCGGGATTCGATACGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGccc                                                                                                                               <  5:21868/148‑1 (MQ=255)
                        gTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCg                                                                                                                              >  3:85522/1‑148 (MQ=255)
                         tATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGc                                                                                                                             >  7:118849/1‑148 (MQ=255)
                             cAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAAcc                                                                                                                         <  6:58428/148‑1 (MQ=255)
                               ggTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGaaa                                                                                                                                                 >  5:158911/1‑122 (MQ=255)
                               ggTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGaaa                                                                                                                                                 <  6:158911/122‑1 (MQ=255)
                               ggTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTg                                                                                                                      <  5:230579/149‑1 (MQ=255)
                                      cAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATc                                                                                                                <  5:92833/148‑1 (MQ=255)
                                         cGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGAt                                                                                                            <  5:235431/149‑1 (MQ=255)
                                              gggATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGcc                                                                                                       >  2:29100/1‑149 (MQ=255)
                                               ggATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCTCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGGCGTCTGGCAGATG‑CCGCAACCGTGCTAATCAGATGCGCCCg                                                                                                     <  6:103319/149‑1 (MQ=255)
                                                    cGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGAtg                                                                                                  <  1:233765/148‑1 (MQ=255)
                                                      aTCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGAtggtg                                                                                               <  3:165196/149‑1 (MQ=255)
                                                       tCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGAtggtgg                                                                                              >  5:392510/1‑149 (MQ=255)
                                                          gCAGGCGATGCCAACCTTCCTCGAAATATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGgtggtg                                                                                            >  2:73762/1‑148 (MQ=255)
                                                               cGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAgtcg                                                                                      <  8:148467/149‑1 (MQ=255)
                                                                 aTGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAgtcgtc                                                                                    >  1:340846/1‑149 (MQ=255)
                                                                  tGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAACTTTATTGACTCACCCGATGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGtcgtcg                                                                                   >  4:198124/1‑149 (MQ=255)
                                                                    ccAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCgg                                                                                  <  1:29100/148‑1 (MQ=255)
                                                                     cAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCgg                                                                                  >  1:58464/1‑147 (MQ=255)
                                                                     cAACCTTCATCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCTCCGGAAATTTTTTTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAAACGTGCTAATCAGTTGCGCCCGATGGTGGTGCAGTCGTCGGaa                                                                                <  4:24986/149‑1 (MQ=255)
                                                                      aaCCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCgg                                                                                  >  6:385562/1‑146 (MQ=255)
                                                                       aCCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGa                                                                              <  8:60173/149‑1 (MQ=255)
                                                                             cTCGAAAAATTACTCGATCAGTCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTa                                                                        <  4:259654/149‑1 (MQ=255)
                                                                             cTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCt                                                                         <  2:340846/148‑1 (MQ=255)
                                                                                gAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCt                                                                     <  8:118849/149‑1 (MQ=255)
                                                                                gAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCt                                                                     <  5:385562/149‑1 (MQ=255)
                                                                                 aaaaaTTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTg                                                                    <  4:85522/149‑1 (MQ=255)
                                                                                 aaaaaTTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTg                                                                    >  4:35962/1‑149 (MQ=255)
                                                                                   aaaTTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGCTGGTGCAGTCGTCGGAAGATTTCTATCTGGc                                                                  >  7:201785/1‑149 (MQ=255)
                                                                                   aaaTTACTCGATCAGGCGCGTTACTCCTCGCGCACGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGc                                                                  >  7:207262/1‑149 (MQ=255)
                                                                                        aCTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAg                                                             >  7:16415/1‑149 (MQ=255)
                                                                                        aCTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAg                                                             >  5:200383/1‑149 (MQ=255)
                                                                                          tCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGc                                                                  >  5:373030/1‑142 (MQ=255)
                                                                                          tCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGc                                                                  <  6:373030/142‑1 (MQ=255)
                                                                                          tCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAgcg                                                           >  8:12833/1‑149 (MQ=255)
                                                                                            gATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCCCCCGTAGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGGGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAgcgcg                                                         >  3:322070/1‑149 (MQ=255)
                                                                                            gATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAgcgcg                                                         >  3:20473/1‑149 (MQ=255)
                                                                                            gATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAgcgc                                                          <  6:392510/148‑1 (MQ=255)
                                                                                             aTCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAgcgcgc                                                        <  3:35962/149‑1 (MQ=255)
                                                                                                aGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCAGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGcaca                                                     <  7:255789/149‑1 (MQ=255)
                                                                                                aGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCAGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGcaca                                                     <  7:249597/149‑1 (MQ=255)
                                                                                                 ggCGCGTTACTCCTCGCGCCCGCCGGAAACTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCacac                                                    <  3:198124/149‑1 (MQ=255)
                                                                                                     cgTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTggg                                                <  6:347261/149‑1 (MQ=255)
                                                                                                         aCTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGAt                                                                             >  4:269852/1‑116 (MQ=255)
                                                                                                         aCTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGAt                                             <  6:163333/148‑1 (MQ=255)
                                                                                                         aCTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATg                                            >  5:203770/1‑149 (MQ=255)
                                                                                                          ctcctcGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCGATCTGGCGAACGCGCGCACACTGGGGCTg                                            >  5:241139/1‑148 (MQ=255)
                                                                                                          ctcctcGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGGCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGt                                           <  1:260153/149‑1 (MQ=255)
                                                                                                           tcctcGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGta                                          >  5:316111/1‑149 (MQ=255)
                                                                                                              tcGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACggg                                                 >  6:23853/1‑139 (MQ=255)
                                                                                                                   cccgccGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCg                                  <  2:95320/149‑1 (MQ=255)
                                                                                                                    cccccgGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCAGATGGTGCTGCAGTCGTCTGAAGATTTCTATCTGTCCAAAGCGCGCACAATGGGGATCTATAATTCCg                                  >  4:37342/4‑148 (MQ=255)
                                                                                                                     cgccgGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGAtggt                                                                                                >  7:48678/1‑85 (MQ=255)
                                                                                                                     cgccgGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGAtggt                                                                                                <  8:48678/85‑1 (MQ=255)
                                                                                                                      gccgGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGAc                               >  3:251575/1‑149 (MQ=255)
                                                                                                                          gAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTaa                           >  6:110132/1‑149 (MQ=255)
                                                                                                                          gAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTaa                           <  8:201785/149‑1 (MQ=255)
                                                                                                                           aaaTTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCTGACGTaa                           <  8:207262/148‑1 (MQ=255)
                                                                                                                           aaaTTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAAc                          <  2:58464/149‑1 (MQ=255)
                                                                                                                                   ttGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATGAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACCAGCTCAc                  <  6:316111/149‑1 (MQ=255)
                                                                                                                                   ttGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACCAGCTCa                   <  6:200383/148‑1 (MQ=255)
                                                                                                                                   ttGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACCAGCTCa                   <  8:16415/148‑1 (MQ=255)
                                                                                                                                               ttGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACCAGCTCACAAGTGATTTGCt      >  1:262881/1‑149 (MQ=255)
                                                                                                                                                   cGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACAAGCTCACCAGTGATTTGCTGGAt  >  5:67007/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGAAGATTTCTATCTGGCGAAAGCGCGCACACTGGGGATGTATAATTCCGGACGTAACCAGCTCACCAGTGATTTGCTGGAT  >  NC_000913/2616692‑2616987

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