Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,692	0	C	A	50.0%	67.3 / 44.3	64	intergenic (+318/+1055)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (4/28);  new base A (32/0);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903549‑2903841                                                                                                                                                |                                                                                                                                                         cAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACTc                                                                                                                                                     >  1:551400/1‑149 (MQ=255)  aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTTTTAAAAACaaatatt                                                                                                                                                  >  2:350291/1‑145 (MQ=255)   aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTa                                                                                                                                                 <  1:288539/151‑1 (MQ=255)    ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑g                                                                                                                                               <  2:176396/151‑1 (MQ=255)     gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAacacccctt                                                                                                                                                  >  1:950038/1‑143 (MQ=255)         atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTA‑Tcccc                                                                                                                                           >  1:659360/1‑148 (MQ=255)           tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCaa                                                                                                                                         <  1:564263/150‑1 (MQ=255)           tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCaa                                                                                                                                         <  1:740673/150‑1 (MQ=255)               tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTccatatcccccaaaa                                                                                                                                       >  2:754706/1‑135 (MQ=255)               tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACACATTA‑TCCTCaaaaaa                                                                                                                                     >  2:78833/1‑148 (MQ=255)                     cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAc                                                                                                                                <  1:662511/149‑1 (MQ=255)                       aTGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTT‑ACCTCCAaacaaaaac‑aaa                                                                                                                            >  1:299018/1‑150 (MQ=255)                         gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACACATTA‑GACTCAAAAAAAAAAC‑aaa                                                                                                                            >  1:790558/1‑148 (MQ=255)                         gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAAAACTAATTT‑ACCTTaaaaaaaaaaa‑aaaa                                                                                                                           >  2:489044/1‑148 (MQ=255)                            tttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAAACCTTA‑AACTCAAAAAAAAAAA‑AAACTTcc                                                                                                                       >  2:757865/1‑150 (MQ=255)                               aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTA‑AGCTCAAAAAAAAAAC‑aaaaattc                                                                                                                       >  1:12247/1‑142 (MQ=255)                                   tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCATTA‑TTCTTCAAACAAAAAA‑AAAcattcaaaatt                                                                                                                 >  2:297924/1‑139 (MQ=255)                                    aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACCCCTCATTT‑ACCTCAAAACAaaaaa‑aaaaatccccatttttt                                                                                                              >  1:248400/1‑137 (MQ=255)                                        gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACtt                                                                                                         <  1:613415/151‑1 (MQ=255)                                          cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTaa                                                                                                       <  1:592690/151‑1 (MQ=255)                                            aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAAcc                                                                                                     <  1:465817/151‑1 (MQ=255)                                             tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAACACATCATTATTCCCAAAAAAAAAAAA‑AAACTTCC‑GA‑TTTTTTACttta                                                                                                       >  2:42566/1‑146 (MQ=255)                                             tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAAAACTCACTA‑TCCTCCAAAAAAAAAC‑AAAAATCCAGA‑ATTATTATTTaaac                                                                                                     >  1:665242/1‑148 (MQ=255)                                               cGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAAt                                                                                                  <  2:187620/151‑1 (MQ=255)                                                gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACACATTA‑GCCTCAAAAAAAAAAC‑AAAAATCCCAA‑TTTATTTCTTAAc                                                                                                      >  1:609542/1‑146 (MQ=255)                                                      cACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTA‑ACCTCAAAAAAAAAAA‑AAAACTTCAGA‑TTTTTTTCTTTAACaaaaaaaaa                                                                                            >  1:287336/1‑143 (MQ=255)                                                       aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACAATCATTT‑TACTCAAAAAAAAAAC‑AAACTTTCAGA‑TTTATTTTTTTAACAAAGAATaa                                                                                            >  2:52550/1‑149 (MQ=255)                                                       aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACa                                                                                          <  2:432861/151‑1 (MQ=255)                                                        cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACACAATT‑CCCTCACAAAAAAAAA‑AAAATTCCAGA‑TTTTTTATTTTAACaaataaaaaaa                                                                                          >  2:441074/1‑141 (MQ=255)                                                        cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACCCATTT‑TCCTCAAAAAAAAAAA‑AAAAATCCCGA‑TTTTTTTCTTTAACaaaaaaaaaaa                                                                                          >  1:637960/1‑141 (MQ=255)                                                           tATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGc                                                                                       <  2:665075/150‑1 (MQ=255)                                                            aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTT‑GCCTCAAAAAAAAAAA‑AAAAATCCCAA‑TTTTTTTCTTTACCAAAAAATAACaa                                                                                         >  1:1130429/1‑147 (MQ=255)                                                                 ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTT‑TCCCCAAAAAAAAAAC‑AAAATTTCCAA‑TTTTTTTTTTTAACAATGAATAaaaaaaataaaa                                                                                 >  2:946702/1‑142 (MQ=255)                                                                    tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTAGTCATCAACAGCGTAGAGAGGCGGTTATGTAAAACCACTCATTA‑GCGTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAat                                                                              <  2:510375/150‑1 (MQ=255)                                                                     cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACACATTT‑ACCTCAAAAAAAAAAA‑AAAAATCCAGA‑TTTTTTTATTTACCAAAGAATAAAAAAATTaaaaaaaa                                                                             >  2:555877/1‑148 (MQ=255)                                                                      aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTATAAAAACAATCAATATCCCCCAAAAAAAAAAC‑AAA‑TTCTTGA‑TTTTTTTTTTAACCAATGAAAAAAAAAACTaaaaaaaa                                                                             >  1:140680/1‑147 (MQ=255)                                                                      aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTA‑ACCTCAAAAAAAAAAA‑AAAACTTCAGA‑TTTTTTACTTTACCAATAAATAACAAACTTaaaaaaaa                                                                             >  1:918033/1‑147 (MQ=255)                                                                          aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTTTAAAAACACCCCTTA‑TCCTCCAAACAAAAACAAAACTTCAAGT‑TTTTTTTTTTAACCAAGAAAAAAAAACCTAaaaaaaaaaatgt                                                                        >  1:503950/1‑144 (MQ=255)                                                                            aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTGTTTAAAACCACCCATCA‑GTCTCTAAAAAAAAAC‑CAACCTCCAGC‑TTTTTTACTTAACAAATAAAAAAAAAGAATAAAAAATATACttttt                                                                     >  1:54770/1‑150 (MQ=255)                                                                            aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTATGTTAAAAAAAACATTT‑TACCCAAAAAAAAAAA‑AAAATTTCCGA‑TTTATTATTTTAACAAAGAATAAAAAAATTAAAAAAAATAATGtt                                                                      >  2:116137/1‑150 (MQ=255)                                                                             gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAt                                                                    <  1:620306/151‑1 (MQ=255)                                                                              tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACCCCCATTA‑TCCTCAAAAAAAAAAA‑AAACATCCCGA‑ATTTTTATTTTAACAAAAAATAAAAAACTTAAAAAAAATACTGTtttt                                                                   >  1:974812/1‑151 (MQ=255)                                                                                tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCATTA‑ACCTCAAAAAAAAAAA‑AAAAATCCAAA‑TTTTTTAATTTACCAATGAATAAAAAACTTAAAAAAAATAATTTTAtttt                                                                 >  1:3447/1‑149 (MQ=255)                                                                                   aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                              <  2:90732/151‑1 (MQ=255)                                                                                      ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                                                             <  1:908765/149‑1 (MQ=255)                                                                                        atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaa                                                         <  2:745405/151‑1 (MQ=255)                                                                                             tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCATTA‑TTCTCAAAAAAAAAAA‑AAAAATCCAGA‑ATTTTTACTTAAACAAAAAAAAAAAAAATTAAAAAAAAAACTTTTTTTTTTAaaaaaaaa                                                       >  2:74794/1‑148 (MQ=255)                                                                                                    aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCa                                             <  1:245009/151‑1 (MQ=255)                                                                                                    aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCa                                             <  1:158984/151‑1 (MQ=255)                                                                                                     tAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTT‑ACCTCCAAAAAAAAAC‑AAAAATTCAGATTTTTTTTTTTAACAAATAAAAAAAAACTTTAAAAAAAAAATTTTTTTTTTaaaaaaaaaa                                                       >  1:157028/1‑141 (MQ=255)                                                                                                      aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAt                                            <  2:1052371/150‑1 (MQ=255)                                                                                                      aaGAACTATTCATCAACAGCGTATAGAGGCGGGTTTGTAAAAACACTCATTA‑GCCCCAAAAAAAAAAA‑AAAAATTCAGA‑TTTATTTATTTAACAATGAAAAAAAAAAATAAAAAAATAAATTTTTTTTTAAAAAAAAAATAAAAAgcatt                                            >  1:677540/1‑146 (MQ=255)                                                                                                         aaCTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACCCCCCTTA‑TACTCCAAAAAAAAAA‑AAAATTCCCAA‑TTTTTTACTTAAACAAAAAATAAAAAAATTTAAAAAAAAATTGTTTTTTTAGAGAAAAAAATAAAAAGCATTaa                                         >  1:972823/1‑150 (MQ=255)                                                                                                             aTTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAACCATTT‑ACCTCAAAAAAAAAAA‑AAAACTTCAGA‑TTTATTTCTTAACCAAATAAAAAAAAAATTAAAAAAAAAAATGTTTTTTTAAAAAAAAAAATAAAAAGCAATAAtt                                       >  1:624543/1‑148 (MQ=255)                                                                                                                catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                 <  2:949797/151‑1 (MQ=255)                                                                                                                catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                 <  1:757165/151‑1 (MQ=255)                                                                                                                     aCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaa                            <  1:314880/151‑1 (MQ=255)                                                                                                                       aGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                          <  1:1129747/151‑1 (MQ=255)                                                                                                                            ataGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataa                      <  2:30945/150‑1 (MQ=255)                                                                                                                             taGAGGCGGTTATTTAAAAACACTCATTA‑GCCCCAAAAAAAAAAA‑AAAAATCCAGA‑ATTATTACTTAAACAAAGAATAAAAAAAATAAAAAAAATAATGTTATTTTTAAGAAAAAAATAAAAAGAATTAATTCAATATTAAAAAaataa                      >  2:502909/1‑149 (MQ=255)                                                                                                                               gagGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa                  <  2:47406/151‑1 (MQ=255)                                                                                                                                agGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa                  <  2:47398/150‑1 (MQ=255)                                                                                                                                  gCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa               <  1:57990/151‑1 (MQ=255)                                                                                                                                   cGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa               >  1:1109039/1‑150 (MQ=255)                                                                                                                                       tATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt           <  1:576527/150‑1 (MQ=255)                                                                                                                                                ccACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACg  >  2:883674/1‑150 (MQ=255)                                                                                                                                                ccACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACg  >  1:281682/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                         CAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTA‑GCCTCAAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903549‑2903841

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