Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	85,481	T→C	R11G (AGA→GGA)	yabR ←	protein YabR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	85,481	0	T	C	90.7%	109.6 / 3.8	43	R11G (AGA→GGA)	yabR	protein YabR
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (22/17);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-02

GGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCG  >  NC_000913/85239‑85714                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          ggCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATgg                                                                                                                                                                                                                                    >  1:133735/1‑250 (MQ=255)  gCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATggag                                                                                                                                                                                                                                  <  2:91854/251‑1 (MQ=255)            gCCGGGCGCTATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTAt                                                                                                                                                                                                                          <  2:542404/249‑1 (MQ=255)                              cATCAGTGGATGGAAGAGCAAATAGTCTAACTTTGCAAACGCTAACTGCTTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATATAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTACTCCTCCGTCATGGAGGATTTTATCGTTTCTTTTCACCTTtcct                                                                                                                                                                                                      >  2:684364/1‑251 (MQ=255)                                           aaGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGACTGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGAGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCtt                                                                                                                                                                                         >  1:464036/1‑251 (MQ=255)                                            agagCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCtta                                                                                                                                                                                        <  1:167569/251‑1 (MQ=255)                                              agCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCttatt                                                                                                                                                                                      >  1:496372/1‑251 (MQ=255)                                               gCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCttatt                                                                                                                                                                                      >  1:664516/1‑250 (MQ=255)                                                 aaTTAGTCTCAATTTGCAAACGCTAACTGATTGCCGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAACAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTCATGGAGGATTTTCTTGTTTCTTTTCACCTTTCCTCATGTTTATTCTTATTAcc                                                                                                                                                                                   >  2:390071/1‑251 (MQ=255)                                                     aGTCTCAATTTGTAAACGCTAACTGATTGCAGAATAGGTTAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGTATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGATGATTTTATCGTTTCTTTTCTCCTTTCCTCCTGTTTATTCTTATTAccctg                                                                                                                                                                                <  2:24935/250‑3 (MQ=255)                                                                    aCGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGc                                                                                                                                                                >  1:578404/1‑251 (MQ=255)                                                                     cGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCt                                                                                                                                                               <  1:497536/251‑1 (MQ=255)                                                                       cTAACTGATTTCATAATACGTCATACATGAATGTCTGTTTTATTCTGCATTTTTTATTGAATTTAGAATTGTATTCTTAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGTGATTCAGCCTATTTATTATCAATTTAATCCTCCGTAATGGAGTATTTTATCGTTTCTTTTCCCCTTTCCTCCTGTTTATTCTTATTCCCCCGTGTTTATGTCTCTgtct                                                                                                                                                               <  1:497519/249‑4 (MQ=255)                                                                             gATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTACTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTg                                                                                                                                                       >  2:113486/1‑251 (MQ=255)                                                                                tGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATATTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTACGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGTTTTTCTCGTTTCTTTTCAATTTTCATCCTGTTTATTCTTTTTACCCCGTGTTTATTTCTCTGGCTGCCATTTgttt                                                                                                                                                    >  2:587286/1‑248 (MQ=255)                                                                                         ggTCAGACATGAATGTCTGGTTTATTCTTCATTGTTTAGTGAATGTAGAATTTTATTCTGAATGTGGGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATGGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATGTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTAATACCCCTCGTATATGGCTCTGGCTGCCAATTGCTTAAGCAAGa                                                                                                                                            <  2:185435/250‑1 (MQ=255)                                                                                          gTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTACTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATc                                                                                                                                          <  1:113486/251‑1 (MQ=255)                                                                                            cAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCg                                                                                                                                         <  1:281571/250‑1 (MQ=255)                                                                                                    aattgcTGGTTTATTATGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATGAAAAAAATAGAGAAATTGCTGTAAGTTGTGTGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTAGCGTTTCTTTTCACCGTTCCTCCTGTTTATTCGTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTaa                                                                                                                                <  2:316861/246‑1 (MQ=255)                                                                                                    aaTGTCTTGTTTTTTCTGCATTTTGTATTGAATGTAGAATTTTAATGTGAATGTTTGGGCTCCCTATATTAGGATTAATTTAAAAAATAGAGAAAGTGCTGTAAGTTGTGGGACTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGTAGGATTTTATCGTTTCGTTTCACCTTTCCTCCTGTTTATTTGTATTACGCCGAGTTTATTTCTCTGGCAGTCAATGGCTTAAGCAAGATCGGACGGTTaa                                                                                                                                <  1:464266/251‑1 (MQ=255)                                                                                                                     gCATTTGTTATTGAATGTAGAATTGTATTCTGAATGGGTGGGCTCTCGATTTTAGGATTAAATAAAAAAATAGAGAAAGTGCTGTAAGGTGTGGGATTCAGCCGATTTATTAGCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGGTTATTATTATTAACCCGTGTTTATGTCTCTGGCTGCCAATTGAGTAAGCAAGATCGGACGGTTAATGTGTTTTACACAtttt                                                                                                               <  2:155227/251‑1 (MQ=255)                                                                                                                                         aaTTTTATTATGAGTGTGTGGGCACTCTAATTTAGGCTCAATTAATAAAATCGAGAAATTGCTGTAAGTTGTGGGATTCAGACGCTTTATTATCAATTTAAAACTCCGTAATGGAGGCTTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACTCCGTGTTTATGTATATGGCTGCCATTTGCTTAAGCTGGATCGGACGTTTACTGTGTTTTAACCATTTTTTcattcaccaagtgaggc                                                                                            >  2:304124/1‑234 (MQ=255)                                                                                                                                             ttATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTATGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACGGTGCGGCAGGcc                                                                                       >  2:501443/1‑251 (MQ=255)                                                                                                                                                 tCTGAATGTGTGGGCTCTCTATTTTAGGCTTAATTAAAAAACTAGAGGAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGTCGGTTAATGTGTTTTACACATTTTTTCCGTCAAACCGTGAGGCAGGCCATgg                                                                                   >  2:203430/1‑251 (MQ=255)                                                                                                                                                        gtgtgGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGAtgtt                                                                            >  2:186952/1‑251 (MQ=255)                                                                                                                                                                     aTTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCg                                                                 >  1:446111/1‑249 (MQ=255)                                                                                                                                                                       tttAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCg                                                                 >  1:273620/1‑247 (MQ=255)                                                                                                                                                                           ggATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGgtc                                                         <  1:594203/251‑1 (MQ=255)                                                                                                                                                                                aattaaAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCg                                                    >  1:601267/1‑251 (MQ=255)                                                                                                                                                                                      aaaaaTAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGCGATGGTCGTCCGATCGCt                                              >  1:561247/1‑251 (MQ=255)                                                                                                                                                                                         aaTAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTAt                                           >  1:18894/1‑251 (MQ=255)                                                                                                                                                                                          aTAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCGCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGGCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATc                                          <  2:255695/251‑1 (MQ=255)                                                                                                                                                                                                  aTTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAggg                                  >  1:91223/1‑251 (MQ=255)                                                                                                                                                                                                    tGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGc                                 <  1:143139/250‑1 (MQ=255)                                                                                                                                                                                                     gCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAAtgt                                                                                                                             <  2:497162/157‑1 (MQ=255)                                                                                                                                                                                                     gCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAAtgt                                                                                                                             >  1:497162/1‑157 (MQ=255)                                                                                                                                                                                                        gTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTa                             >  1:497380/1‑250 (MQ=255)                                                                                                                                                                                                               ggtggATTCAGTCTAGTTATTATCAATTTAATCCTCCGTAATGGAGGAGTTTATCGTTTCTTTTCCCCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTTGCTGCCAATTTCCTAAGCAAGATCTGACGGTTAATGTGTTTTACACATTTTTTCCGGCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGt                      <  2:82251/247‑1 (MQ=255)                                                                                                                                                                                                                 gggATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGTCGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGCTGTTGTCTGGTGCCGAGGTGGTCGTCCGATCGCTTATCGCTCAGGGCGTTACACAAGTAtt                   >  2:176731/1‑251 (MQ=255)                                                                                                                                                                                                                       cAGCCTATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTGTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACGCATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTa             <  1:494091/251‑1 (MQ=255)                                                                                                                                                                                                                            gATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCg        >  2:680914/1‑251 (MQ=255)                                                                                                                                                                                                                                tattatCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAg     <  2:677527/250‑1 (MQ=255)                                                                                                                                                                                                                                  ttatCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTCTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGgcg  >  1:95196/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          GGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCG  >  NC_000913/85239‑85714

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