Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,692	0	C	A	46.8%	53.6 / 57.0	62	intergenic (+318/+1055)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (8/25);  new base A (29/0);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

ACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAC  >  NC_000913/2903548‑2903840                                                                                                                                                 |                                                                                                                                                       aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTc                                                                                                                                                   <  1:109619/150‑1 (MQ=255)       aaTATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAAAAATCATTatac                                                                                                                                            >  1:776594/1‑148 (MQ=255)       aaTATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGcc                                                                                                                                            <  2:848684/151‑1 (MQ=255)         tatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAGcc                                                                                                                                            >  2:784582/1‑149 (MQ=255)            tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACa                                                                                                                                                      >  1:990773/1‑136 (MQ=255)                 acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACTCCtttaactccaaaaaa                                                                                                                                   >  1:234423/1‑137 (MQ=255)                 acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAACCATTAAACTC‑Caaaaaa                                                                                                                                  >  2:578756/1‑150 (MQ=255)                   acTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCCTTAAACTC‑AAAACaaaa                                                                                                                                >  1:990118/1‑149 (MQ=255)                       gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAAACATTAACCTC‑aaaaaaaaaaaaa                                                                                                                            >  1:397987/1‑145 (MQ=255)                           tCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAACCAATCATTATCCTC‑AAaacaaaaacaaa                                                                                                                           >  1:92817/1‑147 (MQ=255)                           tCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCATTTACCTC‑AAAACAaaaaaaaa                                                                                                                           >  1:214104/1‑147 (MQ=255)                            cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTAACCCC‑aaaaaaaaaaaaaa                                                                                                                           >  2:166724/1‑146 (MQ=255)                              ttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACCCATTATACTC‑CAaaaaaaaaaaaaaatcc                                                                                                                      >  2:689047/1‑144 (MQ=255)                               tAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAg                                                                                                                   <  1:1000184/150‑1 (MQ=255)                                aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACCCATTAACCTC‑CAAAAAAAAAAAAAA‑TTTCCGAt                                                                                                                 >  2:52560/1‑151 (MQ=255)                                   aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATtta                                                                                                              <  1:772523/151‑1 (MQ=255)                                   aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATtta                                                                                                              <  1:31495/151‑1 (MQ=255)                                   aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAGCCTC‑AAAAAAAAAACAAAC‑TTTCAGAttt                                                                                                               >  2:505176/1‑150 (MQ=255)                                     aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATttatt                                                                                                            <  2:289584/151‑1 (MQ=255)                                        aGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCTTCAACAGCGTATAGAGGCGGCTATTTAAAAACACTCATTATCCTC‑AAAAAAAAAACAAAA‑ATCCAGAtttttt                                                                                                            >  1:1110815/1‑148 (MQ=255)                                             aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAAATCATTAACCTC‑AAAAAAAAAAAAAACATTCCaatttttt                                                                                                             >  1:563490/1‑136 (MQ=255)                                               aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAAcc                                                                                                    <  2:805198/149‑1 (MQ=255)                                                cGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAAt                                                                                                 <  2:633126/151‑1 (MQ=255)                                                 gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACCATCCCTAGCCTC‑TAAAAAAAAACAAAC‑ATTCCCAGTTAttaattaatc                                                                                                    >  2:874165/1‑145 (MQ=255)                                                  aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCACTAAACTC‑AAAACAAAAACAAAA‑TTTCAGATTTTTTTCtt                                                                                                        >  2:688834/1‑142 (MQ=255)                                                  aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAACCTCAAAAAAAAAAAAAAAC‑TTCC‑GATTTTTTACTTTAACaa                                                                                                  >  2:177325/1‑148 (MQ=255)                                                   aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTCTGTAAAAACACATACTAAACTC‑AAAAAAAAAAAAAAA‑TAACCGCTTTATTAttttta                                                                                                     >  2:848529/1‑141 (MQ=255)                                                      gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTTAACTCAAAAAAAAAAAAAAAC‑TCCCAATTTTTTTATTTAACaaaaaaaa                                                                                             >  1:81975/1‑145 (MQ=255)                                                        aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACa                                                                                         <  1:908782/151‑1 (MQ=255)                                                          tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCCTTAACCTC‑AAAAAAAAAAAAAAA‑TTTCAGATTTTTTATTTTACCAAAGAATAaaaa                                                                                        >  2:418232/1‑150 (MQ=255)                                                           ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                        <  1:755323/149‑1 (MQ=255)                                                             aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTAACCCC‑AAAAAAAAAAAAAAA‑TTCCCGATTTTTTTCTTAAACaaaaaaaaaaaaa                                                                                       >  2:264801/1‑137 (MQ=255)                                                                     tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAat                                                                             <  2:420787/150‑1 (MQ=255)                                                                     tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCATTAGCCTC‑AAAAAAAAAAAAAAA‑CTCCCGATTTTTTACTTAAACAAAGAAAAACAAAATTaaaaaa                                                                              >  2:323971/1‑149 (MQ=255)                                                                      cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAACCATTAGCCCC‑CAAACAAAAACAAAA‑TTTCAGATTTTTTACTTTAACAATTAAAAAAAAACTTaaaaaaaa                                                                            >  2:1116629/1‑148 (MQ=255)                                                                       aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACACATTAGTCTC‑AAAAAAAAAAAAAAC‑ATTCAGATTTATTACTTTAACAAAGAAAAAAAAACTTaaaaaa                                                                              >  1:285537/1‑147 (MQ=255)                                                                            aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                                                                     <  1:521895/151‑1 (MQ=255)                                                                              gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAt                                                                   <  2:291299/151‑1 (MQ=255)                                                                               tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCCTTAAACTC‑AAAACAAAAACAAAC‑TTCCAGATTTTTTTATTTAACAAAGAATAAAAAACTTAAAAAAAATAATTTTAtt                                                                  >  1:350436/1‑151 (MQ=255)                                                                               tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAAAACCCATTAACCTC‑AAAAAAAAAAAAAAA‑ATCCAGATTTTTTACTTTACCAATGAAAAAAAAACTTAAAAAAAATAATGTTAtt                                                                  >  2:495265/1‑151 (MQ=255)                                                                                aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAACCTC‑AAAAAAAAAAAAAAC‑TTCCAGGTTTTTTTCTTAAACAATGAATAAAAAACTTAAAAAAAAAAATGTTAtt                                                                  >  2:266401/1‑150 (MQ=255)                                                                                aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAtc                                                                  <  2:580336/150‑2 (MQ=255)                                                                                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCCTTAACCTC‑CAAAAAAAAAAAAAA‑ATCCCGATTTTTTTCTTTACCAAAGAATAAAAAACTTTAAAAAAATAATGTtttt                                                                  >  2:662492/1‑149 (MQ=255)                                                                                  gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTa                                                               <  2:347875/151‑1 (MQ=255)                                                                                     cTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTTTTTAAAACAACACATTTGACCC‑AAAAAAAAAAAAAACATTCCGGATTT‑TTACATTAACAAAAAATAAAAAAATAAAAAAAAAAACTGTTCtttttaa                                                             >  1:959876/1‑145 (MQ=255)                                                                                           atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTAGCCTC‑CAAAAAAAAAAAAAA‑ATCCAGATTTTTTACTTTAACAAAGAAAAAAAAACTTAAAAAAAATACTGTTTTTTTTaaaaaaaa                                                       >  2:847491/1‑150 (MQ=255)                                                                                                 aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACa                                                <  1:906092/151‑1 (MQ=255)                                                                                                 aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACa                                                <  1:252935/151‑1 (MQ=255)                                                                                                  ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCCTTAGCCTC‑CAAAAAAAAAAAAAA‑TTCCCGATTTTTTACTTAAACAAAGAATAAAAAACTTAAAAAAAATACTGTTATTTTAAAAAAAAAAATaaaa                                                >  2:334184/1‑148 (MQ=255)                                                                                                        aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTa                                         <  2:718569/151‑1 (MQ=255)                                                                                                           aCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAAtt                                      <  2:368654/151‑1 (MQ=255)                                                                                                            cTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAAAACTCCATACCCTC‑CAAACAAAAACAAAC‑ATTCCGATTTATTTTTTTATCAAAAAAAAACAAACTTAAAAAAAATACTGCTTTTTTAAAAAAAAAAATAACAAGCATTaaat                                      >  2:426128/1‑148 (MQ=255)                                                                                                              aTTCATCAACAGCATATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCa                                    <  1:617291/150‑1 (MQ=255)                                                                                                                tcatcaACAGCGTATAGAGGCGGTTATTTAAAAACACTCCTTAACCTC‑AAAAAAAAAACAAAA‑TTCCCGATTTTTTATTTTAACAAAGAAAAAAAAAAATAAAAAAAAAAATTTTTTTTTTAAAAAAAAAATAAAAAGCAATAATTTaa                                   >  2:493911/1‑149 (MQ=255)                                                                                                                      aCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaa                           <  1:96889/151‑1 (MQ=255)                                                                                                                          cGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaat                       <  2:945081/151‑1 (MQ=255)                                                                                                                               agagGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAg                  <  1:639915/151‑1 (MQ=255)                                                                                                                                  ggCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAAAATTAAAAAAATAATAGAtt               <  2:1016321/151‑1 (MQ=255)                                                                                                                                   gCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt               >  1:837052/1‑150 (MQ=255)                                                                                                                                        tATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAtt         >  2:856579/1‑151 (MQ=255)                                                                                                                                             aaaaCCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa    <  2:214108/151‑1 (MQ=255)                                                                                                                                               aaCCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAc  >  2:994051/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                       ACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTC‑AAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAC  >  NC_000913/2903548‑2903840

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