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

AGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAT  >  NC_000913/2903551‑2903832                                                                                                                                              |                                                                                                                                                  aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACa                                                                                                                                                 >  1:1240838/1‑144 (MQ=255) aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCAtt                                                                                                                                           >  2:408418/1‑150 (MQ=255)      tatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTTTTAAAAACACCCATT‑aac                                                                                                                                       >  2:982731/1‑146 (MQ=255)        ttAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACCCAtttac                                                                                                                                        >  1:343616/1‑143 (MQ=255)          aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑aaaa                                                                                                                              <  1:57930/151‑1 (MQ=255)                  tCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAAACCCTT‑ATCCT‑C‑aaaaaaaaaa                                                                                                                        >  1:980590/1‑147 (MQ=255)                   cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑a                                                                                                                     <  2:1115986/150‑1 (MQ=255)                        tCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACACATT‑AACCT‑C‑aaaaaaaaaaa‑aaaaat                                                                                                                >  2:1231317/1‑147 (MQ=255)                          tttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGGTAAAACACCCCCT‑AACCT‑C‑CAaaaaaaaaa‑aaa                                                                                                                   >  2:130336/1‑145 (MQ=255)                           ttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCa                                                                                                             <  2:163835/150‑1 (MQ=255)                              gAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTATATAAAAAAACACCTT‑AATCT‑C‑CAAACAAAAAA‑AAAAATCCa                                                                                                             >  2:130025/1‑147 (MQ=255)                                  aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACAACTCATT‑AATCT‑C‑AAAACAaaaaa‑aaaaattcaactttttt                                                                                                     >  1:466856/1‑137 (MQ=255)                                      gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACtt                                                                                                <  1:821997/151‑1 (MQ=255)                                       aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAATAAAGAC‑AAACTTCCAGA‑TTTATTACTTa                                                                                               <  1:463372/151‑1 (MQ=255)                                             cGAATCGCACTTTATCAATGTTCACTTAAAGTATGAACGGGATATCTTTATTCATAAGAATGATTCATCAACAGCGTATAGAGGCGTTTATGTAAAACCACACATT‑AGACT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCaa                                                                                          <  1:218154/150‑1 (MQ=255)                                              gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATg                                                                                        <  1:1158250/151‑1 (MQ=255)                                                aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGaa                                                                                      <  1:248033/151‑1 (MQ=255)                                                   gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATaa                                                                                   <  2:357423/151‑1 (MQ=255)                                                     aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACa                                                                                 <  2:492818/151‑1 (MQ=255)                                                      cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCCTT‑ATCCT‑C‑AAAAAAAAAAA‑AAACATCCCCATTTTTTTATTTAACAAAt‑aataaaa                                                                                 >  2:386270/1‑148 (MQ=255)                                                        ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTTTAAAAACACTCCTT‑TTTCC‑C‑AAAACAAAAAA‑AAAAATCCAGA‑GTTTTTTCTTAACCAAAGAATAAc                                                                                  >  2:994801/1‑147 (MQ=255)                                                           tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACCATCACT‑TACCT‑C‑CAAAAAAAAAA‑AAAAATCCCCA‑TTTTTTTATTTAACAAAGAATaaaa                                                                                 >  2:794851/1‑143 (MQ=255)                                                              aTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATT‑AGCCT‑C‑CAAAAAAAAAA‑AAACCTTCAGA‑TTTATTAATTAATCaaaaaaaaaaaaac                                                                              >  2:974964/1‑133 (MQ=255)                                                                   cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACAAACCTT‑AACCT‑C‑AAAAAAAAAAA‑AAAATTCCAGA‑TTTTTTATTTAAACAAAAAATAAAaaaaataaaaaata                                                                    >  1:427753/1‑150 (MQ=255)                                                                   cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAACACTCATT‑ATCCT‑C‑CAAACAAAAACAAAACTTCC‑GA‑TTTTTTTTTTTAACAAAGAAAAAAAAACTTtaaaaaaa                                                                    >  2:1131893/1‑148 (MQ=255)                                                                   cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                                                   <  2:714103/151‑1 (MQ=255)                                                                      ttAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACCCATT‑AGCCT‑C‑AAAACAAAAAA‑AAACTTCCAGA‑TTTTTTTTTTTACCAAAAAATAAAAAACTTaaaaaaaa                                                                    >  2:982163/1‑145 (MQ=255)                                                                       tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCCTTAACCCT‑C‑AAAAAAAAAAA‑AAAATTCAGAA‑TTTTTTTTTTAACC‑ATTAAAAAAAAAATTaaaaaa                                                                      >  2:778758/1‑144 (MQ=255)                                                                        aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTg                                                               <  2:173051/150‑1 (MQ=255)                                                                        aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTg                                                               <  1:751041/150‑1 (MQ=255)                                                                        aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTg                                                               <  2:958728/150‑1 (MQ=255)                                                                         aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAACACCCATT‑TACCT‑C‑CAAACAAAAAA‑AAAAATCCCCA‑TTTTTTAATTTTACAAAAAAAAAAAAACCTAAAAAAaataatttt                                                             >  2:1037502/1‑146 (MQ=255)                                                                            tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACCACTCCTT‑AACCT‑C‑CAAAAAAAAAC‑AAAATTCCCCA‑ATTTTTTTTTTACCCAAAAATAAAAAACTTTAAAAAAATAATGTtttt                                                          >  2:544744/1‑151 (MQ=255)                                                                                aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaa                                                      <  2:1166444/151‑1 (MQ=255)                                                                                 aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                     <  2:471120/151‑1 (MQ=255)                                                                                        atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaa                                               <  2:1205444/150‑1 (MQ=255)                                                                                            ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAAAAAAT‑AAACT‑CAAAAAAAAAAAA‑AAATTTCCAAT‑TTTTTAATTTAACAAAAAAAAAAAAAAATTAAAAAAAAAATTTTTTTTTTAaaaaaaaa                                               >  1:13733/1‑147 (MQ=255)                                                                                               ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACTCATT‑AACCTCC‑AAAAAAAAAAA‑AAAATTCC‑GA‑TTTTTTTTTTAACCAAAAAATAAAAAACATAAAAAAAAAAATTTTTTTTTTAAAAAAAAAATaaaa                                        >  2:574603/1‑148 (MQ=255)                                                                                               ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAAATCAAT‑TTCCC‑C‑CAAAAAAAAAA‑AAAAATCCAGA‑TTTATTACTTAAACAAATAAAAAAAAACTTAAAAAAAAAAATGTTTTTTTTAAAAAAAAAATaaaa                                        >  1:318548/1‑148 (MQ=255)                                                                                                 cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGc                                     <  2:701007/151‑1 (MQ=255)                                                                                                    aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCCTT‑TACCT‑C‑AAAAAAAAAAA‑AAAAATTCAGA‑ATTTTTTTTTTAACAATTAAAAAAAAAATTAAAAAAAATAATGTTTTTTTAAAAAAAAAAATaaaaagaatt                                  >  2:1067449/1‑143 (MQ=255)                                                                                                       aaCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAATAAAAAATTAACAGGCATTAAt                               <  1:783309/151‑1 (MQ=255)                                                                                                       aaCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAAt                               <  2:130328/151‑1 (MQ=255)                                                                                                           aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCaa                           <  2:274539/151‑1 (MQ=255)                                                                                                             tcatcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAat                          <  2:126506/150‑1 (MQ=255)                                                                                                              catcaACATCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                         <  2:283893/150‑1 (MQ=255)                                                                                                              catcaACAGCGTATAGAGGCGGTTATGTTAAAACACCCCCT‑AACCC‑C‑AAAAAAAAAAA‑AAAATTTCCCA‑TTTTTTATTTTAACAATAAAAAAAAAACCTTAAAAAAAAAATGTTATTCTTAAGAAAAAAAAAAAAAGCAATAATTTTatat                        >  2:1031665/1‑151 (MQ=255)                                                                                                              catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                         <  2:450643/150‑1 (MQ=255)                                                                                                               atcaACAGCGTATAGAGGCGGTTATATTAAAAAACACCTT‑TGCCT‑C‑CAAAAAAAAAA‑AAAATTCCACA‑TTTATTTCTTTAACAATAAAAAAAAAAAATAAAAAAAATAATTTTATTCTTAAAAAAAAAATAAAAAGCATTAATTTAATAtt                       >  1:20236/1‑151 (MQ=255)                                                                                                                 caACAGCGTATAGAGGCGGTTTTGTAAAAACAATCCTT‑AACCT‑C‑AAAAAAAAAAC‑AAACCTCCCGA‑ATTATTTCTTAAACAAAAAAAAAAAAACTTAAAAAAAATACTTTTTTTTTAAAAAAAAAAATAAAAAGCATTAATTTAatat                        >  2:144502/1‑148 (MQ=255)                                                                                                                 caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaa                     <  1:151563/151‑1 (MQ=255)                                                                                                                  aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaa                    <  2:74796/151‑1 (MQ=255)                                                                                                                  aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaa                    <  1:1235746/151‑1 (MQ=255)                                                                                                                            agagGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAataata           <  1:624550/150‑1 (MQ=255)                                                                                                                            agagGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAataata           <  1:721234/150‑1 (MQ=255)                                                                                                                             gagGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa         <  1:254030/151‑1 (MQ=255)                                                                                                                                gCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt       >  2:825145/1‑150 (MQ=255)                                                                                                                                gCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa      <  2:466883/151‑1 (MQ=255)                                                                                                                                 cGGGTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACATAGAC‑AAACTTCCAGC‑TTTATTACTTAACCAATGACTAACAAGCTAAAAAAATATACTGTTATTCTAGAGAAAAAATTAACACGCATTAATTAAATATGAAAAAAATAATAGATTa      >  1:58831/1‑150 (MQ=255)                                                                                                                                     tATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt  <  2:1086245/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  AGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATT‑AGCCT‑C‑AAAACAAAGAC‑AAACTTCCAGA‑TTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAT  >  NC_000913/2903551‑2903832

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