Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,309 A→G intergenic (+299/‑180) sokC → / → nhaA small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,3090AG97.5% 119.4 / ‑6.1 40intergenic (+299/‑180)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/17);  minor base C (1/0);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTT  >  NC_000913/17076‑17542
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
tAACCCATTCTGTTTTATCGATTCTAATCCTGAAAACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCg                                                                                                                                                                                                                             <  2:484162/250‑1 (MQ=255)
                         aaTCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGcc                                                                                                                                                                                                              <  2:426846/240‑1 (MQ=255)
                         aaTCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGcc                                                                                                                                                                                                              >  1:426846/1‑240 (MQ=255)
                          aTCCTGAAGACGGCTCGCATTTTTGTGGCGTCATTTTTTATTGATTTAATTATTTAACTTTAATTTATCTATCTATAGCAATTATTGCCGACAAGCTGGATTATTTTTGAACTATTGCCCTAACAAATAACGCAGACTGACAAGACATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATATGAGCTTCGGCTCTCTCGTGAGCGCTTACAGCCGTCa‑aatacctc                                                                                                                                                                                                 >  2:329045/1‑245 (MQ=255)
                                aaGACGCCTCGCATTTTTGTGGAGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCa‑aaaaactcatctc                                                                                                                                                                                            >  2:428952/1‑241 (MQ=255)
                                      ccTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTGTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCGCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTg                                                                                                                                                                                     <  1:269619/251‑1 (MQ=255)
                                                 tgtgGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTCTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAAt                                                                                                                                                                          >  2:613947/1‑251 (MQ=255)
                                                  gtgGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAATTATTGTCCTAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTAACCTTAAGCTATGATTTTCTGTGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCACTGATGGCGCAAAtt                                                                                                                                                                         >  2:380938/1‑251 (MQ=255)
                                                          aTTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAAtg                                                                                                                                                                 >  1:119527/1‑250 (MQ=255)
                                                          aTTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAAtg                                                                                                                                                                 <  2:209526/251‑2 (MQ=255)
                                                                tAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGAAATTATTGAAGACAATCTGGTTTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAATAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAATCTATGATTATCTGTGATTAGGGTCACTCGTGCACGCTTACAGCTGTCA‑AAAAACTCTTCTCACCGCTGATGTCGCCATTTCTTCAAtggatcgt                                                                                                                                                           >  2:356044/1‑244 (MQ=255)
                                                                       ttAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTTAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCATAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTAATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaaa                                                                                                                                                     <  2:388717/250‑1 (MQ=255)
                                                                                     ttAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCt                                                                                                                                                                                      <  2:219522/203‑1 (MQ=255)
                                                                                     ttAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCt                                                                                                                                                                                      >  1:219522/1‑203 (MQ=255)
                                                                                     ttAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTac                                                                                                                                       <  2:254263/250‑1 (MQ=255)
                                                                                                   taTCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAttt                                                                                                                        >  1:491399/2‑251 (MQ=255)
                                                                                                    aTCGCAATTAATGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAttt                                                                                                                        <  2:216223/250‑1 (MQ=255)
                                                                                                         aaTTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACgt                                                                                                                  >  2:347926/1‑250 (MQ=255)
                                                                                                            tattGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAAAGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACgatg                                                                                                                >  1:78231/1‑250 (MQ=255)
                                                                                                             attGACGACAAGCTGGATTATTTTTGAAATATTGTCCTATCAAATAACGCCGACTGACAAGAAATTTATTATTACTTTTCCTAATGAATCACTCAGAAATCCTCACCTTAAGCTATGATTATCTGGACTTAGGTTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGACTCAAATTCTTCAATTTATCGTAAAAAACGAGTTATTCCTCCACTATAATCTGATTTTAACgatgat                                                                                                              <  2:342178/251‑1 (MQ=255)
                                                                                                                      acgCTGGCTTATTTTTGAAATATAGACCTACCACTTAAAGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCATCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGTCACAAAC‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAACAAACGCGTTATTCCTACACTATAATCTGATTTTAACgatgctttgtgaggg                                                                                                     >  2:570014/3‑240 (MQ=255)
                                                                                                                         cTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAg                                                                                                                                                                                                                <  1:284/142‑1 (MQ=255)
                                                                                                                         cTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAg                                                                                                                                                                                                                >  2:284/1‑142 (MQ=255)
                                                                                                                                 tttttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGt                                                                                           >  2:420499/1‑250 (MQ=255)
                                                                                                                                      gAAATATTGTCCTAACAAATAACGCCGAATGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTaaaat                                                                                      <  2:205211/250‑2 (MQ=255)
                                                                                                                                         atatTGTCCTAACAAATAACGCCGACTGACAAGCAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGAt                                                                                  >  1:497876/1‑251 (MQ=255)
                                                                                                                                                        ataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGaa                                                                     >  1:39926/5‑249 (MQ=255)
                                                                                                                                                        ataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCACTCAGAAATCCTCACCTTAAGCTCTGTTTATCTCGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGCTTTGCGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGaa                                                                     >  1:414704/5‑249 (MQ=255)
                                                                                                                                                                gACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATaa                                                           <  1:380938/251‑1 (MQ=255)
                                                                                                                                                                       aagaaaTTAATTATTTCTTTTCCTAATTAATCCCTCAGAAAGCCGCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGACA‑AAAAACTCATCTCTCCGCTGATTGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACATTATAAGCTGATTTTAACGATGATTTGGGCGGTGTAAAATAGTAAAAACGATCTAGTCGCCTGAAAGAGAAAGAAAAAGTGa                                                    <  2:564695/248‑1 (MQ=255)
                                                                                                                                                                         gaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGaa                                                   >  1:408485/2‑250 (MQ=255)
                                                                                                                                                                               aattaTTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCt                                             <  1:624850/250‑1 (MQ=255)
                                                                                                                                                                                     taCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCg                                       >  1:386136/1‑250 (MQ=255)
                                                                                                                                                                                            ccTAATGAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCATAGTGTCAATCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGGGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGTTTCTTTa                               <  2:51532/251‑1 (MQ=255)
                                                                                                                                                                                             cTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCttt                                >  2:160993/1‑249 (MQ=255)
                                                                                                                                                                                                  taatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTg                         >  2:6187/1‑251 (MQ=255)
                                                                                                                                                                                                         tCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                  >  2:374203/1‑251 (MQ=255)
                                                                                                                                                                                                          cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGTCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACTCTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATATTAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                  <  2:554746/246‑1 (MQ=255)
                                                                                                                                                                                                          cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCg                 >  1:141432/5‑251 (MQ=255)
                                                                                                                                                                                                                         gaaGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCtt  <  1:595755/250‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
TAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTT  >  NC_000913/17076‑17542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑