Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,364,887 C→T 7.7% A85T (GCC→ACC)  PP_3841 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,364,8870CT 7.7% 227.4 / 6.2 78A85T (GCC→ACC) PP_3841MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base T (2/4);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

ATCGAAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCG  >  NC_002947/4364743‑4365034
                                                                                                                                                |                                                                                                                                                     
atcgaAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTggcggc                                                                                                                                                    >  2:119788/1‑148 (MQ=255)
atcgaAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGgcggcg                                                                                                                                                   <  1:273489/149‑1 (MQ=255)
 tcgaAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGa                                                                                                                                                  >  5:207950/1‑149 (MQ=255)
   gaAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGAc                                                                                                                                                 <  5:347539/148‑1 (MQ=255)
           tCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCAGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGTCGCCGCCAAGATGgc                                                                                                                                        <  8:246867/149‑1 (MQ=255)
            cACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGTCACGCTCCGCGTCCAGAAAAACATAAGCAGCCTGTAGCGTGGTTATCTTCGTTTGAATCTCGGCCATCTTGGTCGCCGCCAAGATGgc                                                                                                                                        <  8:246877/148‑1 (MQ=255)
             aCTGACAGGACTTAAAGTAATTATCTTCCAATTGGCAATAGCACTGTCTTCAAGATTGGCAAGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATATTCGTTTGAATCTCGTCCATCTTGGCGGCGACCAAGAAGgcgc                                                                                                                                      >  2:285773/1‑149 (MQ=255)
              cTGACAGGACTTCAAGGCATTATCTTCCCATTCGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGcc                                                                                                                                     <  7:114868/149‑1 (MQ=255)
               tGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCt                                                                                                                                    >  6:201548/1‑149 (MQ=255)
                    ggACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgc                                                                                                                               <  1:296130/149‑1 (MQ=255)
                    ggACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGATCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTTATCTTCGTTTGAATCTCGGCCATCATGCCGGCGACCAAGACGGCGCCTTgcgc                                                                                                                               >  1:203221/1‑149 (MQ=255)
                     gACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgcg                                                                                                                              >  2:28250/1‑149 (MQ=255)
                     gACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgcg                                                                                                                              >  1:307694/1‑149 (MQ=255)
                          cAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGcaat                                                                                                                          >  7:285049/1‑148 (MQ=255)
                          cAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGcaa                                                                                                                           <  6:35688/147‑1 (MQ=255)
                          cAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGcaa                                                                                                                           <  2:246789/147‑1 (MQ=255)
                              gCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCaatcaatc                                                                                                                     <  6:207950/149‑1 (MQ=255)
                                aTTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCgt                                                                                                                  >  2:95622/1‑149 (MQ=255)
                                        cccATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGAcc                                                                                                                                                <  7:26387/112‑1 (MQ=255)
                                        cccATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGAGTGGTGATCTTCGTTTGAATCTCGGCCAACATGGCGGCCAcc                                                                                                                                                >  8:26387/1‑112 (MQ=255)
                                              ggCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTg                                                                                                     <  1:119788/148‑1 (MQ=255)
                                              ggCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTg                                                                                                     <  1:28250/148‑1 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATTTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTgg                                                                                                    <  1:95622/148‑1 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTgg                                                                                                    <  4:32804/148‑1 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGc                                                                                                   >  7:22193/1‑149 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGc                                                                                                   >  3:120788/1‑149 (MQ=255)
                                                cAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgc                                                                                                                               <  8:119214/121‑1 (MQ=255)
                                                cAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATATTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgc                                                                                                                               >  7:119214/1‑121 (MQ=255)
                                                 aaTAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgc                                                                                                                                 >  1:114740/1‑118 (MQ=255)
                                                 aaTAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgc                                                                                                                                 <  2:114740/118‑1 (MQ=255)
                                                  aTAGCACTGGCTTCAAGCTTGGCACGCTCCGCGACCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGCCGACCAAGACGGCGCCTTGCGCCCAATCAATCGTGTTCCCCCGCCTGGCCAc                                                                                                >  8:101803/1‑149 (MQ=255)
                                                      cACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAg                                                                                             >  1:216768/1‑148 (MQ=255)
                                                        cTGGCTTCAAGCTTGGCACGCTCCGCGTACAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCTACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAAt                                                                                          <  1:285773/149‑1 (MQ=255)
                                                           gCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGTGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGaa                                                                                           <  6:277130/145‑1 (MQ=255)
                                                           gCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGTGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGaa                                                                                           >  5:277130/1‑145 (MQ=255)
                                                           gCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTc                                                                                      >  3:318832/1‑149 (MQ=255)
                                                                   cTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAAtct                                                                              <  8:285049/149‑1 (MQ=255)
                                                                   cTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCACAAtct                                                                              >  8:155641/1‑149 (MQ=255)
                                                                      ggCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCa                                                                            <  2:307694/148‑1 (MQ=255)
                                                                         aCGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGc                                                                          >  2:127454/1‑147 (MQ=255)
                                                                         aCGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGc                                                                          >  5:337675/1‑147 (MQ=255)
                                                                            cTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTg                                                                     <  4:318832/149‑1 (MQ=255)
                                                                                gcgTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCGCAATCTCAGCAAGTGaaaa                                                                 <  2:216768/149‑1 (MQ=255)
                                                                                      aGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGTCTcc                                                           >  2:336475/1‑149 (MQ=255)
                                                                                        aaaaaCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTccc                                                                                                            <  3:302086/99‑1 (MQ=255)
                                                                                        aaaaaCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTccc                                                                                                            >  4:302086/1‑99 (MQ=255)
                                                                                        aaaaaCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGc                                                         >  1:304384/1‑149 (MQ=255)
                                                                                            aCATAAGCGTCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGaa                                                                   <  2:254622/135‑1 (MQ=255)
                                                                                            aCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCAACGCGGCGACCAAGACGGCGCCATGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGaa                                                                   >  1:254622/1‑135 (MQ=255)
                                                                                              aTAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTggg                                                   >  6:92630/1‑149 (MQ=255)
                                                                                                 aGCGGCCTGTAGCGTGGTGATCTTCGTTTGACTCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCa                                                 <  8:91930/148‑1 (MQ=255)
                                                                                                  gCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCAct                                               >  2:249876/1‑149 (MQ=255)
                                                                                                     gCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTg                                            >  4:263056/1‑149 (MQ=255)
                                                                                                      ccTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTCAAAAGCCTCCGCCCTGGGCACTCTGt                                           <  1:183916/149‑1 (MQ=255)
                                                                                                       cTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATCGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGt                                           <  1:249876/149‑1 (MQ=255)
                                                                                                          tAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTAt                                         <  6:337675/147‑1 (MQ=255)
                                                                                                          tAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATg                                        >  6:9890/1‑148 (MQ=255)
                                                                                                             cGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTTGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGaaa                                     <  7:155641/148‑1 (MQ=255)
                                                                                                                 gtgATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGt                                           <  7:287487/138‑1 (MQ=255)
                                                                                                                 gtgATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGt                                           >  8:287487/1‑138 (MQ=255)
                                                                                                                       ttCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAg                           >  2:75571/1‑148 (MQ=255)
                                                                                                                          gTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCaa                        <  4:120788/148‑1 (MQ=255)
                                                                                                                           tttGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACg                      >  5:344632/1‑149 (MQ=255)
                                                                                                                           tttGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACg                      <  2:304384/149‑1 (MQ=255)
                                                                                                                           tttGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACg                      <  5:9890/149‑1 (MQ=255)
                                                                                                                             tGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACgg                     <  8:22193/148‑1 (MQ=255)
                                                                                                                              gAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCGCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGt                    <  3:263056/148‑1 (MQ=255)
                                                                                                                                aTCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTg                 >  5:168398/1‑149 (MQ=255)
                                                                                                                                 tctcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGt                >  3:286596/1‑149 (MQ=255)
                                                                                                                                   tcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCAct                                               <  5:124378/116‑1 (MQ=255)
                                                                                                                                   tcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCAct                                               >  6:124378/1‑116 (MQ=255)
                                                                                                                                   tcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGtt               >  3:251487/1‑148 (MQ=255)
                                                                                                                                     ggCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCa             <  6:344632/148‑1 (MQ=255)
                                                                                                                                     ggCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCa             <  5:201548/148‑1 (MQ=255)
                                                                                                                                     ggCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCGGTAGGAAAGCCAAGTCAGCAACGGTCTGTTCAc            >  5:279736/1‑149 (MQ=255)
                                                                                                                                     cgCCGTCTTGGTCGCCGCCAAGATGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGt                    <  5:199711/140‑1 (MQ=255)
                                                                                                                                     cgCCGTCTTGGTCGCCGCCAAGATGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGt                    >  6:199711/2‑141 (MQ=255)
                                                                                                                                       ccATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTg          >  3:160128/1‑149 (MQ=255)
                                                                                                                                       ccATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTg          <  7:101803/149‑1 (MQ=255)
                                                                                                                                           cTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTa        >  1:65649/1‑148 (MQ=255)
                                                                                                                                              ggcggcGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTg                 >  5:130694/1‑135 (MQ=255)
                                                                                                                                              ggcggcGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTg                 <  6:130694/135‑1 (MQ=255)
                                                                                                                                               gcggcgACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCg  <  1:127454/149‑1 (MQ=255)
                                                                                                                                                cggcgACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCg  >  8:160534/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ATCGAAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTGTTCCCCCGCCTGGCCACTAGAATGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCG  >  NC_002947/4364743‑4365034

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: