Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,359,064 T→A 5.0% E402D (GAA→GAT nusA ← transcription termination/antitermination factor L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,359,0640TA 5.0% 248.9 / 2.3 79E402D (GAA→GATnusAtranscription termination/antitermination factor L
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (2/2);  total (37/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 = 7.08e-01

CAGTTCAGCCGCCAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATGTCACCGGTTT  >  NC_002947/5358920‑5359212
                                                                                                                                                 |                                                                                                                                                      
cAGTTCAGCCGCCAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGt                                                                                                                                                    <  5:277065/149‑1 (MQ=255)
 aGTTCAGCCGCCAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGGGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTc                                                                                                                                                   <  6:117165/149‑1 (MQ=255)
  gTTCAGCCGCCAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGAAGACATCGTc                                                                                                                                                   >  7:134357/1‑148 (MQ=255)
     cAGCCGCCAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGAAGACATCGTCaat                                                                                                                                                <  8:134357/148‑2 (MQ=255)
        ccgccAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCa                                                                                                                                            <  7:198165/149‑1 (MQ=255)
          gccAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTGTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCAt                                                                                                                                           >  1:99236/1‑149 (MQ=255)
           ccAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATc                                                                                                                                          >  4:186552/1‑148 (MQ=255)
               gTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGt                                                                                                                                     >  5:252181/1‑149 (MQ=255)
                tCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGtt                                                                                                                                    <  1:42552/149‑1 (MQ=255)
                tCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGtt                                                                                                                                    <  6:252181/149‑1 (MQ=255)
                    tNTCCATGCNNTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTg                                                                                                                                   <  8:19128/146‑1 (MQ=255)
                    tGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGa                                                                                                                                  <  8:92967/147‑1 (MQ=255)
                        cATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCAttt                                                                                                                            <  7:231471/149‑1 (MQ=255)
                        cATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCAtt                                                                                                                             >  8:231471/1‑148 (MQ=255)
                              cTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATttcttc                                                                                                                        >  2:305939/1‑147 (MQ=255)
                              cTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATttcttc                                                                                                                        <  1:305939/147‑1 (MQ=255)
                              cTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTcc                                                                                                                       >  3:257908/1‑148 (MQ=255)
                              cTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCAATCGATGTTGAGCATTTCTTcc                                                                                                                       >  6:187158/1‑148 (MQ=255)
                              cTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGACATCGATGTTGAGCATTTCTTcc                                                                                                                       >  8:95260/1‑148 (MQ=255)
                                    ggAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCaaa                                                                                                                                                >  2:290069/1‑117 (MQ=255)
                                    ggAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCaaa                                                                                                                                                <  1:290069/117‑1 (MQ=255)
                                    ggAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTa                                                                                                               <  3:186552/149‑1 (MQ=255)
                                       gagCTGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGt                                                                                                            <  7:32422/149‑1 (MQ=255)
                                         gCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGtt                                                                                                           <  5:95473/148‑2 (MQ=255)
                                            ggTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGc                                                                                                        >  4:152899/1‑148 (MQ=255)
                                             gTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGt                                                                                                                <  5:309132/139‑1 (MQ=255)
                                             gTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGt                                                                                                                >  6:309132/1‑139 (MQ=255)
                                             gTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCaa                                                                                                      >  5:185249/1‑149 (MQ=255)
                                                  cGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATttc                                                                                                  <  4:204030/148‑1 (MQ=255)
                                                    gCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGTTGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATttctt                                                                                                >  1:220259/1‑148 (MQ=255)
                                                      cGGGTGGTCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCaa                                                                                             <  4:257908/149‑1 (MQ=255)
                                                      cGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCaa                                                                                             <  4:105293/149‑1 (MQ=255)
                                                      cGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCaa                                                                                             <  6:185249/149‑1 (MQ=255)
                                                       gggTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAg                                                                                            <  3:152899/149‑1 (MQ=255)
                                                       gggTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAg                                                                                            <  6:278527/149‑1 (MQ=255)
                                                       gggTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTAAGTAGGCAATTTCTTCAAg                                                                                            >  8:99525/1‑149 (MQ=255)
                                                               gTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAAAGGTACGTAGGAAATTTCTAAAAGTCTtgtg                                                                                     >  2:221794/1‑148 (MQ=255)
                                                                tCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGaa                                                                                   <  2:99236/149‑1 (MQ=255)
                                                                    ccAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCt                                                                               >  8:20047/1‑149 (MQ=255)
                                                                    ccAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCt                                                                               >  1:284697/1‑149 (MQ=255)
                                                                             tcttcGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGa                                                                        <  2:226887/148‑1 (MQ=255)
                                                                                   gTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGa                                                                        >  6:198874/1‑142 (MQ=255)
                                                                                   gTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGa                                                                        <  5:198874/142‑1 (MQ=255)
                                                                                    tAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGCGCTCATTGACGATGTCTTCGTCAAAGCCCTCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACt                                                               >  6:46468/1‑149 (MQ=255)
                                                                                    tAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACt                                                               >  4:27609/1‑149 (MQ=255)
                                                                                       cGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTgg                                                            <  2:220259/149‑1 (MQ=255)
                                                                                         aTGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGATGGTGAAGCCTTCGTCGACCAGCACTTGGGc                                                          >  7:213306/1‑149 (MQ=255)
                                                                                            gCCTTTGTCAACAAACGGTCCTTTGCTCGAGCGCGGATCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTTAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAg                                                       <  1:221794/149‑1 (MQ=255)
                                                                                             ccTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGAc                                                                       >  2:134710/1‑133 (MQ=255)
                                                                                               ttGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGctc                                                    >  1:254018/1‑149 (MQ=255)
                                                                                                  gTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGctcctc                                                 <  2:254018/149‑1 (MQ=255)
                                                                                                  gTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGctcctc                                                 <  3:27609/149‑1 (MQ=255)
                                                                                                  gTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGctcctc                                                 >  2:218740/1‑149 (MQ=255)
                                                                                                   tCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCtcg                                                >  8:103204/1‑149 (MQ=255)
                                                                                                     aacaaACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCtcgtc                                              >  5:224006/1‑149 (MQ=255)
                                                                                                     aacaaACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCtcgtc                                              >  2:95547/1‑149 (MQ=255)
                                                                                                          aCGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACtt                                         >  4:113258/1‑149 (MQ=255)
                                                                                                             gTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTc                                        >  5:148242/1‑147 (MQ=255)
                                                                                                             gTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTc                                        <  6:148242/147‑1 (MQ=255)
                                                                                                                    gCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCATAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCt                                <  3:52531/148‑1 (MQ=255)
                                                                                                                    gCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCt                                <  1:218740/148‑1 (MQ=255)
                                                                                                                        gAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGCTGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTtcgt                                                                           <  4:253045/102‑1 (MQ=255)
                                                                                                                        gAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTtcgt                                                                           >  3:253045/1‑102 (MQ=255)
                                                                                                                        gAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGa                             <  7:20047/147‑1 (MQ=255)
                                                                                                                        gAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGAt                            >  6:297663/1‑148 (MQ=255)
                                                                                                                          gcgaGGATCTCATTGACGATGTCTTCGTCAAAGCCATCGATTTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATg                           <  8:213306/147‑1 (MQ=255)
                                                                                                                                 tggctttGACGAAGACATCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTgcg                   <  2:198695/143‑1 (MQ=255)
                                                                                                                                 tggctttGACGAAGACATCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTgcg                   >  1:198695/6‑148 (MQ=255)
                                                                                                                                     aTTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCt                                >  6:315811/1‑131 (MQ=255)
                                                                                                                                     aTTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGa              <  2:284697/149‑1 (MQ=255)
                                                                                                                                      ttGAAGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCATGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGAt             <  3:37563/149‑1 (MQ=255)
                                                                                                                                         aCGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATGTc          <  7:184766/149‑1 (MQ=255)
                                                                                                                                          cGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAg                                     >  1:73506/1‑122 (MQ=255)
                                                                                                                                          cGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAg                                     <  2:73506/121‑1 (MQ=255)
                                                                                                                                           gATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAg                >  6:170710/1‑141 (MQ=255)
                                                                                                                                           gATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAg                <  5:170710/141‑1 (MQ=255)
                                                                                                                                           gATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGACAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATTTCAc        >  7:257528/1‑149 (MQ=255)
                                                                                                                                            aTGTCTTCGTCAAAGCCATCTATGTTTAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCTTGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATGTCAc        <  1:95547/149‑1 (MQ=255)
                                                                                                                                            aTGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATGTCAcc       <  7:95260/149‑1 (MQ=255)
                                                                                                                                                 ttCGTCAAAGACATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATGTCACCGGttt  >  6:135157/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CAGTTCAGCCGCCAGGTCTTTGTCCATGCCCTCAAGGGAGAGCAGGTCTTCGGCCGGGTGGGCGTCTGCCAGTTTTTCTTCGGTAGCGATGGCCTTGGTCAACAAACGGTCCTTGGCTCGAGCGCGGAGCTCATTGACGATGTCTTCGTCAAAGCCATCGATGTTGAGCATTTCTTCCAACGGTACGTAGGCAATTTCTTCAAGGCTGGTGAAGCCTTCGTCGACCAGCACTTGGGCCAGCTCCTCGTCGACTTCCAGTTCCTCGATGAAATTGCGCAGGATGTCACCGGTTT  >  NC_002947/5358920‑5359212

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: