Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,841 A→T 6.0% E22V (GAA→GTA) ‡ mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,8410AT 6.0% 204.9 / 2.4 67E22V (GAA→GTA) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base T (2/2);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCC  >  NC_002947/635696‑635986
                                                                                                                                                 |                                                                                                                                                 
aCCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAc                                                                                                                                                 <  6:262511/148‑1 (MQ=255)
aCCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAc                                                                                                                                                 <  2:171376/148‑1 (MQ=255)
aCCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACt                                                                                                                                                >  1:271907/1‑149 (MQ=255)
  cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTg                                                                                                                                              <  3:154193/149‑1 (MQ=255)
  cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTg                                                                                                                                              <  1:46783/149‑1 (MQ=255)
   aGGGGCGACTGAGGGAAACAGCGTGTGCGATTTGTCGCGCCGTTGGGTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCAGATTCTCGGTATTTGCGGGACTTTCATGGGTTCGCTGGGGGTGCTGGCCAAGTAACTTgg                                                                                                                                             <  3:16107/149‑1 (MQ=255)
       gCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCAc                                                                                                                                         <  2:312897/149‑1 (MQ=255)
        cGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCagc                                                                                                                                        <  5:41441/149‑3 (MQ=255)
        cGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCAcc                                                                                                                                        <  7:121339/149‑1 (MQ=255)
           cTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcgcg                                                                                                                                     >  7:344498/1‑149 (MQ=255)
              aGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGt                                                                                                                                    <  1:66855/147‑1 (MQ=255)
                 gAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                          >  2:60302/1‑138 (MQ=255)
                 gAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCa                                                                                                                                          <  1:60302/138‑1 (MQ=255)
                      cgCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGGGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                          <  1:311370/148‑1 (MQ=255)
                      aGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                          >  6:271969/1‑149 (MQ=255)
                        cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGa                                                                                                                        <  3:346556/149‑1 (MQ=255)
                         gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGAAGTTTTGCGAGGTCGACCATGCATATTCACATTATAGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACAGGCTCTGAc                                                                                                                       >  1:233254/1‑149 (MQ=255)
                           gTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                       >  7:168051/1‑147 (MQ=255)
                           gTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                       <  8:168051/147‑1 (MQ=255)
                           gTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGc                                                                                                                     >  3:296493/1‑149 (MQ=255)
                                 tttGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGt                                                                                                               <  6:324009/149‑1 (MQ=255)
                                    gTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                          >  4:76857/1‑135 (MQ=255)
                                    gTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGctct                                                                                                                          <  3:76857/135‑1 (MQ=255)
                                       gcgcCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATcc                                                                                                         <  1:167762/149‑1 (MQ=255)
                                           ccTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGATCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCg                                                                                                     <  1:262570/149‑1 (MQ=255)
                                                gCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAg                                                                                                >  8:69245/1‑149 (MQ=255)
                                                 cTACCCTGCCTCGATGTTTTGCGAGTTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGc                                                                                               <  6:303231/149‑1 (MQ=255)
                                                     ccTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                  >  1:85218/1‑94 (MQ=255)
                                                     ccTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                  <  2:85218/94‑1 (MQ=255)
                                                      cTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCa                                                                                          >  6:68998/1‑149 (MQ=255)
                                                         ccTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCt                                                                                       >  6:367062/1‑149 (MQ=255)
                                                          cTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTc                                                                                      >  8:205002/1‑149 (MQ=255)
                                                          cTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTc                                                                                      >  5:42991/1‑149 (MQ=255)
                                                           tCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGaccaccaggtctg                                                                                     >  8:132343/1‑142 (MQ=255)
                                                             gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                   <  2:271907/149‑1 (MQ=255)
                                                                 ttttGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGccc                                                                               <  1:79687/149‑1 (MQ=255)
                                                                   ttGCGAGGTCGACCATGCATATTCACATTCTCGGTAATTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                              <  8:329612/148‑1 (MQ=255)
                                                                    tGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAg                                                                             <  1:364338/148‑1 (MQ=255)
                                                                     gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                            <  1:134095/148‑1 (MQ=255)
                                                                     gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                            <  3:165794/148‑1 (MQ=255)
                                                                     gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                            <  3:171497/148‑1 (MQ=255)
                                                                      cGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAggg                                                                           >  2:209178/1‑148 (MQ=255)
                                                                         ggTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                  <  6:241401/74‑1 (MQ=255)
                                                                         ggTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                  >  5:241401/1‑74 (MQ=255)
                                                                              aCCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACt                                                                  <  5:271969/149‑1 (MQ=255)
                                                                                  tGCATATTCACATTCTCGGTATTTGCTGCTCTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCGACCGCGTCTCCGGCTCTGACGCCAACGTCTATCCCCCGATGAGGAGCCAGCTCGAAGCCCTGGGTATCGAAGTGAcc                                                              <  2:233254/149‑1 (MQ=255)
                                                                                     atatTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGAccc                                                             >  3:141854/1‑147 (MQ=255)
                                                                                      tatTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAg                                                          <  4:141854/149‑1 (MQ=255)
                                                                                           acaTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                     >  5:120864/1‑149 (MQ=255)
                                                                                           acaTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAACCCCAGGGTATCGAACTGACCCAAGGCTa                                                      <  1:209178/148‑1 (MQ=255)
                                                                                                     tATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCTCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGcc                                            >  2:207494/1‑148 (MQ=255)
                                                                                                       tttGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                              >  6:132671/1‑112 (MQ=255)
                                                                                                       tttGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                              <  5:132671/112‑1 (MQ=255)
                                                                                                              cACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                  >  1:327300/1‑149 (MQ=255)
                                                                                                               aCTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATcc                                 >  3:134407/1‑149 (MQ=255)
                                                                                                                    cATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGaa                                                                    >  6:335568/1‑109 (MQ=255)
                                                                                                                    cATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGaa                                                                    <  5:335568/109‑1 (MQ=255)
                                                                                                                     aTGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                            <  6:292045/148‑1 (MQ=255)
                                                                                                                               tGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgc                              >  5:173922/1‑136 (MQ=255)
                                                                                                                               tGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgc                              <  6:173922/136‑1 (MQ=255)
                                                                                                                                 gCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggcg                 <  6:120864/147‑3 (MQ=255)
                                                                                                                                     tGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATc           <  5:68998/149‑1 (MQ=255)
                                                                                                                                      ggtGGCCAAGTTCCTTGGCCAGCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCg          >  6:142140/3‑149 (MQ=255)
                                                                                                                                       cTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg         <  7:69245/149‑1 (MQ=255)
                                                                                                                                        tGGCCAAGGAACTTGGCCACCGCGGCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc        <  6:42991/149‑1 (MQ=255)
                                                                                                                                           ccAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAAc     >  7:258303/1‑149 (MQ=255)
                                                                                                                                            cAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCg                                <  4:295979/121‑1 (MQ=255)
                                                                                                                                            cAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCg                                >  3:295979/1‑121 (MQ=255)
                                                                                                                                            cAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACg    >  4:270494/1‑149 (MQ=255)
                                                                                                                                             aaGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCCACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCCGCCCAGCTCGATCCGGAGCCAGACCTCGAGGTCATCGGCAACGc   >  5:304098/1‑149 (MQ=255)
                                                                                                                                              aGGAAGTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCGCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGcc  <  8:170/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCC  >  NC_002947/635696‑635986

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: