Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,692	0	C	A	43.1%	79.6 / 42.2	58	intergenic (+318/+1055)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (5/28);  new base A (25/0);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 8.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAC  >  NC_000913/2903542‑2903840                                                                                                                                                        |                                                                                                                                                        gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAAc                                                                                                                                                         <  2:153192/151‑1 (MQ=255)     ttACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTAAAAACa                                                                                                                                                       >  1:958518/1‑149 (MQ=255)        cAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTAAAAacacac                                                                                                                                                    >  1:14815/1‑147 (MQ=255)          aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTa                                                                                                                                               <  1:400040/151‑1 (MQ=255)           ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTTAAAACACTCA‑tt                                                                                                                                                >  2:162746/1‑149 (MQ=255)              atatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCt                                                                                                                                          <  1:460556/151‑1 (MQ=255)                   aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCaaa                                                                                                                                      <  1:378445/150‑1 (MQ=255)                    aCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTAAAAACACTCA‑TTA‑ACCCCaaaa                                                                                                                                     >  1:81456/1‑150 (MQ=255)                                     tAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TAATGTTAAACCCACCC‑TTC‑GCCCCGAAACAAAAAAAAAATTcctg                                                                                                                    >  1:840795/1‑148 (MQ=255)                                           aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTTTTAAAAAAAC‑CA‑TTT‑TCCTCCAAAAAAAAAAAAAAATCCAG‑ATttatt                                                                                                            >  2:176044/1‑151 (MQ=255)                                           aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAAAAACCCA‑TTA‑ACCTCAAAAAAAAAACAAAATTCCAG‑Atttttt                                                                                                            >  1:460797/1‑151 (MQ=255)                                               gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACtt                                                                                                        <  2:687945/151‑1 (MQ=255)                                                   aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTAAAAACAAACA‑TTA‑TCCCCAAAACAAAAAAAAAAATCCAC‑ATTTATTAAtt                                                                                                        >  2:779049/1‑147 (MQ=255)                                                   aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTTAAAACAACCA‑TTA‑AGCTCAAAAAAAAAAAAAAATTTCAA‑AGttatt                                                                                                            >  1:104567/1‑143 (MQ=255)                                                     aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCa                                                                                                   <  1:176066/150‑1 (MQ=255)                                                        aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATATAAAAACACCCA‑TTA‑TCCTCCAAAAAAAAAAAAAATTTCAG‑ATTTATTACTTAATCaa                                                                                                  >  1:422987/1‑148 (MQ=255)                                                           cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTAAAAACACCCA‑TTT‑ACCTCAAAAAAAAAAAAAAATTCCCG‑ATTTTTTACTTAAACaaaaaa                                                                                              >  1:277951/1‑145 (MQ=255)                                                                 ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTTAAAACAATCC‑TTA‑ACCCCCAAAAAAAAAAAAAAATTCAG‑ATTTATTAATTTAACAAAGAATAaaaa                                                                                        >  1:254586/1‑149 (MQ=255)                                                                   aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCtt                                                                                    <  2:498691/151‑1 (MQ=255)                                                                    tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTAAAAACACCCATTTA‑ACCCCAAAAAAAAAAAAAATTCCCAA‑TTTTTTTATTAAACCattaataaaaaacat                                                                                     >  1:593098/1‑136 (MQ=255)                                                                    tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTAAAAACACCCA‑TTT‑ACCTCAAAAAAAAAAAAAAAATCCAG‑ATTTTTTTATTAAACAATAAATAaaaaaatt                                                                                    >  1:958592/1‑146 (MQ=255)                                                                      aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTaaa                                                                                 <  2:73040/151‑1 (MQ=255)                                                                      aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTaaa                                                                                 <  2:16474/151‑1 (MQ=255)                                                                         tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTaaaaa                                                                               <  1:605958/150‑1 (MQ=255)                                                                            cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTAAAAACACAAC‑TTT‑TCCTCCAAAAAAAAAAAAAACTCCCG‑AATTTTTTCTTTAACAATGAAAAAAAAAGATaaaaaa                                                                              >  2:979824/1‑148 (MQ=255)                                                                                tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATAc                                                                         <  1:631729/149‑1 (MQ=255)                                                                                  aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                                                                     <  2:603584/151‑1 (MQ=255)                                                                                   aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGt                                                                      <  2:461821/149‑1 (MQ=255)                                                                                   aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAAACACTCA‑TTA‑ACCTCAAAACAAAAACAAAATTCCAG‑ATTTTTTATTTAAACAAAGAAAAAAAAACTTTAAAAAAATACTGt                                                                      >  1:461872/1‑149 (MQ=255)                                                                                    gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTTAAAACACACA‑TTA‑ACCTCCAAACAAAAAAAAACTTCCAG‑ATTTTTTTCTTAAACAAAGAATAACAAAACTAAAAAAAATAATGtttt                                                                   >  2:933402/1‑149 (MQ=255)                                                                                    gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATTTAAAAACAAACA‑TTA‑ACCCCCAAAAAAAAACAAAAATCCAG‑ATTTTTTTTTTAAACAAAGAAAAAAAAACTTTAAAAAAAAAATGtt                                                                     >  1:813730/1‑149 (MQ=255)                                                                                     tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTTAAAAAACCCA‑TTA‑ACCCCCAAAAAAAAACAAAATTCCCA‑ATTTTTTTTTTTAACAAATAATAAAAAACTTAAAAAAaataattttttt                                                                  >  2:848842/1‑143 (MQ=255)                                                                                       tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑GTGTGTTAAAACCCACA‑ATT‑TACTCAAAAAAAAAAAAAACTTTCCC‑AGTTTTTTAATAAACAAAGAAAAAAAAACTTTAAAAAAATACTAttt                                                                    >  1:117986/1‑146 (MQ=255)                                                                                       tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑CTTTATAAAAAAACAAA‑ATT‑TCGTTAAAAAAAAAAAAAAATTCCCC‑ATTTATTATTTTAACAAAGAAAAAAAAAATAAAAAAAAAAAATattatttt                                                                >  2:594714/1‑149 (MQ=255)                                                                                        gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCt                                                                <  2:87250/150‑1 (MQ=255)                                                                                          aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                             <  1:34678/151‑1 (MQ=255)                                                                                            tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTGTTAAAAAAAACA‑ATA‑ACCTCCAAAAAAAAACAAACCTTCAA‑ATTTTTTAAATTAACAAAAAAAAAAAAAATTAAAAAAAAAAATTTTATTTTaaaa                                                            >  2:152075/1‑149 (MQ=255)                                                                                                      tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTTTTTAAAAACACACA‑TTATCCCCCAAAAAAAAAAAAAAATTCCGA‑ATTTTTTCTTTAACAAAAAAAAAAAAAATTAAAAAAAAAAATTTTTTTTTTAAAAAAAAAAta                                                    >  1:572124/1‑148 (MQ=255)                                                                                                       aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACa                                                <  1:179788/151‑1 (MQ=255)                                                                                                           aTAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCa                                            <  2:81445/151‑1 (MQ=255)                                                                                                           aTAAGAACTATTCATCAACAGCGTATAGAGGCGG‑ATCAGTCAAAACACTGA‑CTA‑ATCTCCAAAACAAGAAAAACGTCCAGAATTTTTTAATTAA‑CAACGAAAAACAAATCAAAAAAAAATACAGTTTTTTTAAATAAAAAAATaaaaaac                                             >  1:81011/1‑145 (MQ=255)                                                                                                             aaGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAtt                                          <  2:742826/151‑1 (MQ=255)                                                                                                             aaGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAt                                           <  2:239434/150‑1 (MQ=255)                                                                                                                    aTTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTc                                     <  1:894898/149‑1 (MQ=255)                                                                                                                      tcatcaACAGCGTATAGAGGCGG‑TTATTTAAAAACACACA‑TTA‑ACCTCAAAAAAAAAACAAAAATCCAG‑ATTTTTTTCTTAAACAAAGAATAAAAAACTTTAAAAAAAAAATTTTTTTTTTAAAAAAAAAATAAAAAGCAATaa                                        >  1:206105/1‑144 (MQ=255)                                                                                                                       catcaACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                <  2:70785/151‑1 (MQ=255)                                                                                                                       catcaACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                <  1:611417/151‑1 (MQ=255)                                                                                                                       catcaACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                                 <  2:788080/150‑1 (MQ=255)                                                                                                                         tcaACAGCGTATAGAGGCGG‑TTTTGTAAAAACACTCA‑TTA‑TTCCCCAAAAAAAAAAAAACATCCAG‑ATTTATTTTTTAAACAAAAAATAAAAAACTTAAAAAAAATAATTTTTTTTTTAAAAAAAAAATAAAAAACATTAATTaaa                                   >  2:195105/1‑146 (MQ=255)                                                                                                                            aCAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCGATATTaaaa                           <  1:673452/151‑1 (MQ=255)                                                                                                                                   ataGAGGCGG‑TTATGTAAAAACAACCA‑CTT‑GCCTCCAAAAAAAAAAAAAAATCCCC‑ATTTTTTTTTTTAACAAAAAAAAAAAAACTTAAAAAAAATAATGTTATTTTTTAAAAAAAAATAAAAAGCATTAATTAAATATTaaaaaaaaaa                     >  2:332257/1‑150 (MQ=255)                                                                                                                                     agagGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAataata                   <  2:141078/150‑1 (MQ=255)                                                                                                                                      gagGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa                 <  2:161814/151‑1 (MQ=255)                                                                                                                                          cGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa             <  1:787591/151‑1 (MQ=255)                                                                                                                                              ttATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaa                          <  1:944930/135‑1 (MQ=255)                                                                                                                                              ttATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaa                          >  2:944846/1‑135 (MQ=255)                                                                                                                                                aTGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt        <  2:890323/151‑1 (MQ=255)                                                                                                                                                  gTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt      >  2:703319/1‑151 (MQ=255)                                                                                                                                                     aaaCCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa    >  2:575520/1‑150 (MQ=255)                                                                                                                                                      aaCCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAc  >  2:793782/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGG‑TTATGTAAAACCACTCA‑TTA‑GCCTCAAAACAAAGACAAACTTCCAG‑ATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAC  >  NC_000913/2903542‑2903840

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