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

TACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903547‑2903841                                                                                                                                                   |                                                                                                                                                         tACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑t                                                                                                                                                     <  2:670891/150‑1 (MQ=255)     aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑AAAAACAC‑TCaat                                                                                                                                                 >  2:958746/1‑148 (MQ=255)     aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑AAAAACAC‑ACAtt                                                                                                                                                 >  2:31659/1‑150 (MQ=255)     aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTa                                                                                                                                                <  1:925733/151‑1 (MQ=255)     aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACCC‑TCAtt                                                                                                                                                 >  2:77955/1‑150 (MQ=255)           atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑a                                                                                                                                         <  2:978291/151‑1 (MQ=255)               aCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑ACATTAGCCTC‑AAAAc                                                                                                                                     <  2:480362/151‑1 (MQ=255)                       cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CCCTTAACCCC‑CAaaaaaaa                                                                                                                                 >  2:398332/1‑146 (MQ=255)                         aTGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑TCATTTTACTC‑aaaaaaaaaaaaaa                                                                                                                            >  1:379957/1‑150 (MQ=255)                                   aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGAttt                                                                                                                 <  1:203072/151‑1 (MQ=255)                                     tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑TAAAACAC‑TCATTAACCTC‑AAAAAAAAAAAAAAAATTCAGAttttt                                                                                                               >  2:266619/1‑149 (MQ=255)                                      aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CCCTTAGCCCC‑AAAAAAAAAAAAAAATTCCAGAAttttt                                                                                                              >  2:51785/1‑151 (MQ=255)                                      aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑TAAAACAA‑TCATTAGCCCC‑AAAAAAAAAAAAAAATTCCAGAtttttt                                                                                                              >  1:692504/1‑151 (MQ=255)                                      aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATttatt                                                                                                              <  1:721120/151‑1 (MQ=255)                                       aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATttatta                                                                                                             <  2:629003/151‑1 (MQ=255)                                       aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAAAAC‑CCATTTACCTC‑CAAAAAAAAACAAAAATCCCGAtttttt                                                                                                              >  2:614709/1‑150 (MQ=255)                                          gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑AAAAAAAA‑CCATTATTCTC‑AAAAAAAAAAAAAAAATCCAGAATTATTACtt                                                                                                          >  1:313058/1‑151 (MQ=255)                                            cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑TCATTAACCCC‑CAAAAAAAAAAAAAAATCCAAATTTTTTACtt                                                                                                          >  2:153617/1‑149 (MQ=255)                                             aaTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAAAAC‑TCAATTACCCC‑AAAAAAAAAAAAAAAATCCACATTTTTTACtt                                                                                                          >  1:273353/1‑148 (MQ=255)                                              aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAAcc                                                                                                      <  1:16775/151‑1 (MQ=255)                                               tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑CTATTAACCTC‑CAAAAAAAAAAAAACCTCCAGAATTTTTTCtt                                                                                                          >  2:41659/1‑146 (MQ=255)                                                        cACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAAAAC‑ACATTAACCTC‑AAAACAAAAAAAAACATTCAGATTTATTTATTAAACAATaaaaaa                                                                                             >  1:730545/1‑150 (MQ=255)                                                            ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑AAAAACAC‑TCATTAGCCTC‑CAAAAAAAAAAAAAAATCCAGATTTTTTTCTTAAACAATAAATAaa‑aa                                                                                         >  2:720778/1‑149 (MQ=255)                                                            ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑AAAAACAC‑CCATTAGCCTC‑AAAAAAAAAAAAAAATTCCAGAATTATTTCTTAAACAATGAATAaa‑aa                                                                                         >  2:164378/1‑149 (MQ=255)                                                            ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑aa                                                                                         <  1:614500/149‑1 (MQ=255)                                                              aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCt                                                                                      <  2:878915/150‑1 (MQ=255)                                                                    tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTaaaaa                                                                                <  1:33696/150‑1 (MQ=255)                                                                            aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGt                                                                       <  1:340180/151‑1 (MQ=255)                                                                                 aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTc                                                                  <  2:778330/151‑1 (MQ=255)                                                                                   gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCt                                                                 <  2:782806/150‑1 (MQ=255)                                                                                    aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACAC‑CCATTAACCTC‑CAAAAAAAAAAAAAAATCCCGATTTTTTTCTTAAACAAAAAAAAAA‑AAACTTAAAAAAAATAATTTTAtttt                                                                 >  2:675350/1‑147 (MQ=255)                                                                                     aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑AAAAACAC‑TCATTAGCCTC‑CAAAAAAAAACAAAATTCCCGATTTTTTTTTTTAACAAAGAAAAAA‑AAACTTAAAAAAAAAAATGTTTTTTTaa                                                               >  2:729259/1‑150 (MQ=255)                                                                                         gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTT‑TAAAACAA‑ACATTTTACTC‑AAAACAAAAAAAAAACTTCAGGTTTTTTAATTAAACAAATAATAAC‑AAACATAAAAAAAATACTTTTATTTTTAAGaa                                                           >  1:226007/1‑150 (MQ=255)                                                                                            atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGaaaaa                                                        <  2:585374/150‑1 (MQ=255)                                                                                                ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑AAAAACAC‑CCATTAGCCTC‑AAAAAAAAAAAAAACATCCCCAATTTTTACTTTAACAAAAAAAAAA‑AAACTTTAAAAAAATAATGTTTTTTTTAaaaaaaaa                                                       >  2:197214/1‑147 (MQ=255)                                                                                                ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAtt                                                     <  2:746135/149‑1 (MQ=255)                                                                                                  aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAAc                                                  <  2:854608/150‑1 (MQ=255)                                                                                                  aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATAT‑AAAAAAACACCATTAACCCC‑AAAA‑AAAAAAAAAATTTCAGATTTATTACTTTAACAAAAAATAAA‑AAAAATAAAAAAAATAAAGTTATTTTAAAAAAAAAAATaa                                                   >  1:690996/1‑149 (MQ=255)                                                                                                       tAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑TCATTAACCTC‑AAAAAAAAAAAAAAATTCCCAAATTATTTCTTAAACAAAGAAAAAA‑AAACTTAAAAAAAATAATGTTTTTTTAAAAAAAAAAATAAAAAGCa                                             >  2:775115/1‑150 (MQ=255)                                                                                                            aCTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACAC‑ACCTTATCCTC‑AAAACAAAAAAAAACTTTCCGATTTATTAATTTACCAATAAATAAA‑AAAATTAAAAAAAATAATGTTATTTTTAAGAAAAAAATaaaaagaaata                                          >  1:272258/1‑139 (MQ=255)                                                                                                             cTATTCATCAACAGCGTATAGAGGCGGTTATGT‑TAAAACCC‑TCCTTTTCCTC‑AAAAAAAAAACAAACATCCAGATTTTTTTATTAACCAAAAAAAAACAAAACTTAAAAAA‑ATACTATTTTTTTTAAAAAAAAAATaaaacacatt                                           >  1:1010901/1‑138 (MQ=255)                                                                                                              tATTCATCAACAGCGTATAGAGGCGGCTTTGT‑TAAAACAC‑CCATTAGGCTC‑AAAACAAAAAAAAAAATCCAGAATTATTTCTTTAACAAAGAATAAA‑AAACTTTAAAAAAATACTTTTTTTTTTAACAAAAAAATAAAAAGCAATAAttt                                      >  2:286241/1‑149 (MQ=255)                                                                                                               aTTCATCAACAGCGTATAGAGGCGGATATATAAAAAACAT‑CCATTCG‑CTC‑CAAAAAAAAAGAAAAATCCTGAATTTTTAATTAATCAACAAATAAC‑AAAGTTAAAAAAAAACATGTTATTTTTAAGAAAAAAATAAAAAGCATTAAtt                                       >  2:111433/1‑148 (MQ=255)                                                                                                                 tcatcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAat                                   <  1:741764/150‑1 (MQ=255)                                                                                                                  catcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                 <  1:77951/151‑1 (MQ=255)                                                                                                                  catcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                 <  1:248564/151‑1 (MQ=255)                                                                                                                   atcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAtt                                <  1:740507/151‑1 (MQ=255)                                                                                                                   atcaACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGACTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAtt                                <  1:977650/151‑1 (MQ=255)                                                                                                                     caACAGCGTATAGAGGCGGTTTTTT‑TAAAACAC‑TCACTAACCCCAAAAAAAAAAAAAAAATCCCACTTTTTTTCTTTAAACAA‑GAATAAC‑AAACTTAAAAAAAAAACTTTTTTTCTTAAGAAAAAAATAAAAAGCATTAATTTCatat                                 >  2:544458/1‑148 (MQ=255)                                                                                                                     caACAGCGTATAGAGGCGGTTATTT‑AAAAACAC‑TCCTTAGCCTC‑AAAACAAAAAAAAACATTCAGATTTTTTACTTAAACAATGAATAAC‑AAACTTAAAAAAAATAATGTTTTTTTTAAGAAAAAATTAAAAAGCATTAATTTAATAtt                                >  1:187827/1‑149 (MQ=255)                                                                                                                        cAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaa                           <  1:31657/151‑1 (MQ=255)                                                                                                                         aGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                          <  2:692608/151‑1 (MQ=255)                                                                                                                            gTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaata                       <  2:25695/151‑1 (MQ=255)                                                                                                                              ataGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATGCTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataat                     <  1:41657/151‑1 (MQ=255)                                                                                                                                agagGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAataata                    <  1:69211/150‑1 (MQ=255)                                                                                                                                 gagGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAA‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataat                     >  2:819011/1‑148 (MQ=255)                                                                                                                                     cagTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGTCATTAATTCAATATTAAAAAAATAATAGATTa               <  1:95232/148‑1 (MQ=255)                                                                                                                                     cGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa               >  1:21800/1‑150 (MQ=255)                                                                                                                                          aTGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt         <  1:473181/151‑1 (MQ=255)                                                                                                                                             t‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt      >  1:377865/1‑151 (MQ=255)                                                                                                                                             t‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt       <  1:51781/150‑1 (MQ=255)                                                                                                                                                aaaCCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa    <  2:387792/151‑1 (MQ=255)                                                                                                                                                aaaCCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa     >  2:855164/1‑150 (MQ=255)                                                                                                                                                   ccAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACg  >  2:224853/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                         TACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCAC‑TCATTAGCCTC‑AAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903547‑2903841

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