Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,364,891 G→C 5.9% V83V (GTC→GTG)  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,8910GC 5.9% 242.5 / 3.7 80V83V (GTC→GTG) PP_3841MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base C (2/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATCGAAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCGATAG  >  NC_002947/4364743‑4365038
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                               <  2:246789/147‑1 (MQ=255)
                          cAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGcaa                                                                                                                               <  6:35688/147‑1 (MQ=255)
                              gCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCaatcaatc                                                                                                                         <  6:207950/149‑1 (MQ=255)
                                aTTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCg‑t                                                                                                                      >  2:95622/1‑149 (MQ=255)
                                        cccATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGAcc                                                                                                                                                    <  7:26387/112‑1 (MQ=255)
                                        cccATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGAGTGGTGATCTTCGTTTGAATCTCGGCCAACATGGCGGCCAcc                                                                                                                                                    >  8:26387/1‑112 (MQ=255)
                                              ggCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTg                                                                                                         <  1:119788/148‑1 (MQ=255)
                                              ggCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTg                                                                                                         <  1:28250/148‑1 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATTTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTgg                                                                                                        <  1:95622/148‑1 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTgg                                                                                                        <  4:32804/148‑1 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGc                                                                                                       >  7:22193/1‑149 (MQ=255)
                                               gCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGc                                                                                                       >  3:120788/1‑149 (MQ=255)
                                                cAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgc                                                                                                                                   <  8:119214/121‑1 (MQ=255)
                                                cAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATATTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgcgc                                                                                                                                   >  7:119214/1‑121 (MQ=255)
                                                 aaTAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgc                                                                                                                                     <  2:114740/118‑1 (MQ=255)
                                                 aaTAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTgc                                                                                                                                     >  1:114740/1‑118 (MQ=255)
                                                  aTAGCACTGGCTTCAAGCTTGGCACGCTCCGCGACCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGCCGACCAAGACGGCGCCTTGCGCCCAATCAATCG‑TGTTCCCCCGCCTGGCCAc                                                                                                    >  8:101803/1‑149 (MQ=255)
                                                      cACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAg                                                                                                 >  1:216768/1‑148 (MQ=255)
                                                        cTGGCTTCAAGCTTGGCACGCTCCGCGTACAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCTACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAt                                                                                              <  1:285773/149‑1 (MQ=255)
                                                           gCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGTGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGaa                                                                                               <  6:277130/145‑1 (MQ=255)
                                                           gCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGTGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGaa                                                                                               >  5:277130/1‑145 (MQ=255)
                                                           gCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTc                                                                                          >  3:318832/1‑149 (MQ=255)
                                                                   cTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAAtct                                                                                  <  8:285049/149‑1 (MQ=255)
                                                                   cTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCACAAtct                                                                                  >  8:155641/1‑149 (MQ=255)
                                                                      ggCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCa                                                                                <  2:307694/148‑1 (MQ=255)
                                                                         aCGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGc                                                                              >  2:127454/1‑147 (MQ=255)
                                                                         aCGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGc                                                                              >  5:337675/1‑147 (MQ=255)
                                                                            cTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTg                                                                         <  4:318832/149‑1 (MQ=255)
                                                                                gcgTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCGCAATCTCAGCAAGTGaaaa                                                                     <  2:216768/149‑1 (MQ=255)
                                                                                      aGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGTCTcc                                                               >  2:336475/1‑149 (MQ=255)
                                                                                        aaaaaCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTccc                                                                                                                >  4:302086/1‑99 (MQ=255)
                                                                                        aaaaaCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTccc                                                                                                                <  3:302086/99‑1 (MQ=255)
                                                                                        aaaaaCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGc                                                             >  1:304384/1‑149 (MQ=255)
                                                                                            aCATAAGCGTCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGaa                                                                       <  2:254622/135‑1 (MQ=255)
                                                                                            aCATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCAACGCGGCGACCAAGACGGCGCCATGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGaa                                                                       >  1:254622/1‑135 (MQ=255)
                                                                                              aTAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTggg                                                       >  6:92630/1‑149 (MQ=255)
                                                                                                 aGCGGCCTGTAGCGTGGTGATCTTCGTTTGACTCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCa                                                     <  8:91930/148‑1 (MQ=255)
                                                                                                  gCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCAct                                                   >  2:249876/1‑149 (MQ=255)
                                                                                                     gCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTg                                                >  4:263056/1‑149 (MQ=255)
                                                                                                      ccTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTCAAAAGCCTCCGCCCTGGGCACTCTGt                                               <  1:183916/149‑1 (MQ=255)
                                                                                                       cTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAATCGTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGt                                               <  1:249876/149‑1 (MQ=255)
                                                                                                          tAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTAt                                             <  6:337675/147‑1 (MQ=255)
                                                                                                          tAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATg                                            >  6:9890/1‑148 (MQ=255)
                                                                                                             cGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTTGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGaaa                                         <  7:155641/148‑1 (MQ=255)
                                                                                                                 gtgATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGt                                               >  8:287487/1‑138 (MQ=255)
                                                                                                                 gtgATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGt                                               <  7:287487/138‑1 (MQ=255)
                                                                                                                       ttCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAg                               >  2:75571/1‑148 (MQ=255)
                                                                                                                          gTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCaa                            <  4:120788/148‑1 (MQ=255)
                                                                                                                           tttGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACg                          >  5:344632/1‑149 (MQ=255)
                                                                                                                           tttGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACg                          <  2:304384/149‑1 (MQ=255)
                                                                                                                           tttGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACg                          <  5:9890/149‑1 (MQ=255)
                                                                                                                             tGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACgg                         <  8:22193/148‑1 (MQ=255)
                                                                                                                              gAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCGCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGt                        <  3:263056/148‑1 (MQ=255)
                                                                                                                                aTCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTg                     >  5:168398/1‑149 (MQ=255)
                                                                                                                                 tctcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGt                    >  3:286596/1‑149 (MQ=255)
                                                                                                                                   tcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCAct                                                   <  5:124378/116‑1 (MQ=255)
                                                                                                                                   tcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCAct                                                   >  6:124378/1‑116 (MQ=255)
                                                                                                                                   tcGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGtt                   >  3:251487/1‑148 (MQ=255)
                                                                                                                                     ggCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCa                 <  6:344632/148‑1 (MQ=255)
                                                                                                                                     ggCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCa                 <  5:201548/148‑1 (MQ=255)
                                                                                                                                     ggCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCGGTAGGAAAGCCAAGTCAGCAACGGTCTGTTCAc                >  5:279736/1‑149 (MQ=255)
                                                                                                                                     cgCCGTCTTGGTCGCCGCCAAGATGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGt                        <  5:199711/140‑1 (MQ=255)
                                                                                                                                     cgCCGTCTTGGTCGCCGCCAAGATGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGt                        >  6:199711/2‑141 (MQ=255)
                                                                                                                                       ccATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTg              >  3:160128/1‑149 (MQ=255)
                                                                                                                                       ccATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTg              <  7:101803/149‑1 (MQ=255)
                                                                                                                                           cTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCGTTGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTa            >  1:65649/1‑148 (MQ=255)
                                                                                                                                              ggcggcGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTg                     >  5:130694/1‑135 (MQ=255)
                                                                                                                                              ggcggcGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTg                     <  6:130694/135‑1 (MQ=255)
                                                                                                                                               gcggcgACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCg      <  1:127454/149‑1 (MQ=255)
                                                                                                                                                cggcgACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCg      >  8:160534/1‑148 (MQ=255)
                                                                                                                                                   cgACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCGATAg  <  5:92630/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCGAAAATAGTCACTGACAGGACTTCAAGGCATTATCTTCCCATTGGCAATAGCACTGGCTTCAAGCTTGGCACGCTCCGCGTCCAGAAAAACATAAGCGGCCTGTAGCGTGGTGATCTTCGTTTGAATCTCGGCCATCTTGGCGGCGACCAAGACGGCGCCTTGCGCGCAATCAATCG‑TGTTCCCCCGCCTGGCCACTAGAAT‑GTCCCCAATCTCAGCAAGTGAAAAGCCTCCGCCCTGGGCACTCTGTATGAAAGCCAAGTCAGCAACGGTCTGTTCACTGTAACTGCGATAG  >  NC_002947/4364743‑4365038

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