Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,376,390 A→T 10.8% L392Q (CTG→CAG)  dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,376,3900AT10.8% 175.3 / 10.7 65L392Q (CTG→CAG) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base T (4/3);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGG  >  NC_002947/5376243‑5376533
                                                                                                                                                   |                                                                                                                                               
agcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                                                                                                                                                >  1:148163/1‑149 (MQ=255)
agcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                                                                                                                                                >  3:178986/1‑149 (MQ=255)
 gcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                                                                                                                                                >  8:147584/1‑148 (MQ=255)
  caCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTTGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                                                                                                                                                >  5:100156/1‑147 (MQ=255)
  caCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                                                                                                                                                >  8:5617/1‑147 (MQ=255)
     gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcag                                                                                                                                             >  5:297028/1‑147 (MQ=255)
      tGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                                                                                                                                          <  1:177490/149‑1 (MQ=255)
            gTCACGGCGCTCTGGTTGTCGTCGGCAGTTGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTc                                                                                                                                    <  6:100156/149‑1 (MQ=255)
             tCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCt                                                                                                                                   <  6:297028/149‑1 (MQ=255)
                  gcgcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCAc                                                                                                                              <  5:160287/149‑1 (MQ=255)
                              tcgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCa                                                                                                                   <  4:178986/148‑1 (MQ=255)
                               cgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc                                                                                                                     >  6:51571/1‑145 (MQ=255)
                               cgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc                                                                                                                     <  5:51571/145‑1 (MQ=255)
                                  cgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAcc                                                                                                              >  7:39640/1‑149 (MQ=255)
                                   ggCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcagc                                                                                                                                 <  6:164436/129‑4 (MQ=255)
                                   ggCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcagc                                                                                                                                 >  5:164436/1‑126 (MQ=255)
                                   ggCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCg                                                                                                             <  1:104740/149‑1 (MQ=255)
                                    gCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCTATACACAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTAACGTCACCGGCCAGTCCCGc                                                                                                            >  7:124680/1‑149 (MQ=255)
                                           agaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTgg                                                                                                     <  3:91576/149‑1 (MQ=255)
                                           agaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTg                                                                                                      >  7:144420/1‑148 (MQ=255)
                                            gaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc                                                                                                                     <  8:31647/132‑1 (MQ=255)
                                            gaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc                                                                                                                     >  7:31647/1‑132 (MQ=255)
                                                 ccTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                                                               <  7:147584/149‑1 (MQ=255)
                                                  cTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATg                                                                                                  <  6:224409/145‑1 (MQ=255)
                                                  cTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATg                                                                                                  >  5:224409/1‑145 (MQ=255)
                                                  cTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcag                                                                                              >  5:192978/1‑149 (MQ=255)
                                                           cttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCAt                                                                                       <  6:192978/147‑1 (MQ=255)
                                                              ggTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATgg                                                                                     <  2:199799/146‑1 (MQ=255)
                                                              ggTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCaa                                                                                  >  5:89584/1‑149 (MQ=255)
                                                               gTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCaa                                                                                  <  8:254221/148‑1 (MQ=255)
                                                                  gggATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCaa                                                                                  <  6:228404/145‑1 (MQ=255)
                                                                  gggATGGTGGTGTTATTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCaa                                                                                  >  5:228404/1‑145 (MQ=255)
                                                                   ggATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCt                                                                             <  8:39640/149‑1 (MQ=255)
                                                                          ggtgTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCg                                                                      >  1:255759/1‑149 (MQ=255)
                                                                           gtgtTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCg                                                                      <  7:5617/148‑1 (MQ=255)
                                                                                cttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGTAACGGCTTCGTCCGGGTTg                                                                 <  8:144420/148‑1 (MQ=255)
                                                                                cttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGa                                                                >  7:174205/1‑149 (MQ=255)
                                                                                  tctcGATCAGCGCGGTCATCACGCCGCACATGGTTTCGATACCCAGGGTCAGCGTGCTGACGTCCAGCAGCACGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGTATGCCAGAACCCAAAGAAACGGCATCGTCCGGGTTGAc                                                               >  8:270196/1‑148 (MQ=255)
                                                                                   ctcGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCAcc                                                                                                         >  8:257454/1‑105 (MQ=255)
                                                                                      gATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTcc                                                           <  2:148163/148‑1 (MQ=255)
                                                                                           gcgcGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCggg                                                                    >  1:258968/1‑134 (MQ=255)
                                                                                           gcgcGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCggg                                                                    <  2:258968/134‑1 (MQ=255)
                                                                                               ggTCATCACGCCGCCCATGGTTTCGTTACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGTCCAGTACCGCACGCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACTTCCTTGCGCGctt                                                 <  6:279972/149‑1 (MQ=255)
                                                                                                  catcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGcttctt                                              >  2:257811/1‑149 (MQ=255)
                                                                                                   atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGcttctt                                              >  1:64078/1‑148 (MQ=255)
                                                                                                   atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCttt                                             <  8:174205/149‑1 (MQ=255)
                                                                                                     caCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAc                                           >  2:4723/1‑149 (MQ=255)
                                                                                                         ccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       <  2:274370/149‑1 (MQ=255)
                                                                                                         ccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       <  4:152014/149‑1 (MQ=255)
                                                                                                          cgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaag                                      <  2:64078/149‑1 (MQ=255)
                                                                                                            cccATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaagaa                                    <  1:233653/149‑1 (MQ=255)
                                                                                                                        aTACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtc                        >  2:233503/1‑149 (MQ=255)
                                                                                                                           cccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtctt                      <  2:255759/148‑1 (MQ=255)
                                                                                                                            ccAGGGTCAGCGGGCTGACGTCCAGCGGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttctttc                    >  4:189850/1‑149 (MQ=255)
                                                                                                                            ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTcc                                                           <  4:237740/110‑1 (MQ=255)
                                                                                                                            ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTcc                                                           >  3:237740/1‑110 (MQ=255)
                                                                                                                            ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttctttc                    >  2:2025/1‑149 (MQ=255)
                                                                                                                                ggTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCa                <  1:257811/149‑1 (MQ=255)
                                                                                                                                 gTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCAc               <  6:182867/149‑1 (MQ=255)
                                                                                                                                 gTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCAc               <  1:233503/149‑1 (MQ=255)
                                                                                                                                     gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGc           >  7:196520/1‑149 (MQ=255)
                                                                                                                                        ggCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGc        <  6:89584/149‑1 (MQ=255)
                                                                                                                                        ggCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGc        <  1:4723/149‑1 (MQ=255)
                                                                                                                                        gaaagACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCgg         >  8:225653/5‑148 (MQ=255)
                                                                                                                                           tGACGTCCAGCAGCAGCACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGCAAGCACACATGGCAACGGCTTCGTCCGGGTTGAAGTCCTTGCGCGCTTCTATACATAAAAAGTCGGCCACTACTTTCTGCACCAGCGCCATg     >  6:125021/1‑149 (MQ=255)
                                                                                                                                            gACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       >  8:139697/1‑114 (MQ=255)
                                                                                                                                            gACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       >  8:143791/1‑114 (MQ=255)
                                                                                                                                            gACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       <  7:143791/114‑1 (MQ=255)
                                                                                                                                            gACGGCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       <  7:139697/114‑1 (MQ=255)
                                                                                                                                             aCGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCACTTCTTTACCGAAGAAGt                                  >  1:190925/1‑118 (MQ=255)
                                                                                                                                             aCGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCACTTCTTTACCGAAGAAGt                                  <  2:190925/118‑1 (MQ=255)
                                                                                                                                              agTCCAGCCGCATGAAGTCTGTCAAGTCACCTGACAGTAACGCACCCTGTATGGCATCACCCATTGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTCCCGAAGAATGCGGCCACTTCTTTCTGCACCAGCGGCATGCgg  <  8:124680/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGG  >  NC_002947/5376243‑5376533

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: