Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,376,396 A→T 10.8% L390Q (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,3960AT10.8% 175.1 / 10.9 65L390Q (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 = 5.49e-01

TGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGG  >  NC_002947/5376249‑5376533
                                                                                                                                                   |                                                                                                                                         
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                                                                                                                     <  5:51571/145‑1 (MQ=255)
                         cgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc                                                                                                                     >  6:51571/1‑145 (MQ=255)
                            cgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAcc                                                                                                              >  7:39640/1‑149 (MQ=255)
                             ggCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcagc                                                                                                                                 >  5:164436/1‑126 (MQ=255)
                             ggCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcagc                                                                                                                                 <  6:164436/129‑4 (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                                                                                                                     >  7:31647/1‑132 (MQ=255)
                                      gaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc                                                                                                                     <  8:31647/132‑1 (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                                                                    <  2:258968/134‑1 (MQ=255)
                                                                                     gcgcGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCggg                                                                    >  1:258968/1‑134 (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                                                           >  3:237740/1‑110 (MQ=255)
                                                                                                                      ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTcc                                                           <  4:237740/110‑1 (MQ=255)
                                                                                                                      ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttctttc                    >  2:2025/1‑149 (MQ=255)
                                                                                                                          ggTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCa                <  1:257811/149‑1 (MQ=255)
                                                                                                                           gTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCAc               <  1:233503/149‑1 (MQ=255)
                                                                                                                           gTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCAc               <  6:182867/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:143791/1‑114 (MQ=255)
                                                                                                                                      gACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       >  8:139697/1‑114 (MQ=255)
                                                                                                                                      gACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       <  7:143791/114‑1 (MQ=255)
                                                                                                                                      gACGGCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                                       <  7:139697/114‑1 (MQ=255)
                                                                                                                                       aCGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCACTTCTTTACCGAAGAAGt                                  <  2:190925/118‑1 (MQ=255)
                                                                                                                                       aCGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCACTTCTTTACCGAAGAAGt                                  >  1:190925/1‑118 (MQ=255)
                                                                                                                                        agTCCAGCCGCATGAAGTCTGTCAAGTCACCTGACAGTAACGCACCCTGTATGGCATCACCCATTGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTCCCGAAGAATGCGGCCACTTCTTTCTGCACCAGCGGCATGCgg  <  8:124680/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGG  >  NC_002947/5376249‑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: ‑