Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,350	0	C	T	43.5%	57.3 / 35.1	62	S64F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (8/27);  new base T (27/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.53e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGT  >  NC_000913/4604201‑4604497                                                                                                                                                       |                                                                                                                                                     aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑tttttctttttc                                                                                                                                                      <  1:246252/150‑1 (MQ=255) aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑Tttttctttttct                                                                                                                                                     <  1:1149593/151‑1 (MQ=255)    aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGcc                                                                                                                                                  <  2:924915/151‑1 (MQ=255)     aaaggcggtaTGAGTATTAACGGTATTGAGATTTTATTAAGAAAGTTTGCAGGTAACCTTAATAGGCTCAATGACTAAACGGGTTTGGTTTCTTTTCACGATACCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTGTTTCTGcgg                                                                                                                                                 <  2:1097496/142‑3 (MQ=255)      aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCAt                                                                                                                                                <  2:1019184/151‑1 (MQ=255)      aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCAt                                                                                                                                                <  2:704204/151‑1 (MQ=255)       aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCT‑TTTTTCTTTTTTTGCCAt                                                                                                                                                >  1:190193/1‑150 (MQ=255)              aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                         <  2:1140220/150‑1 (MQ=255)                gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                                      <  2:1014303/151‑1 (MQ=255)                    agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTATTTTTTTTCCATGAGAAGTGAgcg                                                                                                                                    >  1:1173802/1‑149 (MQ=255)                      taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                                 <  2:786096/150‑1 (MQ=255)                       aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                               <  1:437737/151‑1 (MQ=255)                       aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                                <  2:924801/150‑1 (MQ=255)                         cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTCTTTTTTTTCCGTGGAAAAGGAGCGCggag                                                                                                                               >  1:1134287/1‑149 (MQ=255)                         cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTTCCTGAGAAGTGAGCGCag                                                                                                                                 >  2:58292/1‑147 (MQ=255)                          ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGAT‑TTTTTTTTTTTTTTCCATGAGGAGGGTGCggagggag                                                                                                                             >  2:23823/1‑150 (MQ=255)                          ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTTTTCCATGAGAAGTGAGCGCAGGGAgg                                                                                                                            >  1:1190526/1‑151 (MQ=255)                                gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                      <  1:583177/151‑1 (MQ=255)                                       tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAT‑TTTTTTTTTTTTTTTCATTAGGAGTGGAGGCGGAGAgggggtatttt                                                                                                                   >  2:513173/1‑139 (MQ=255)                                        aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCct                                                                                                              <  2:806826/150‑1 (MQ=255)                                        aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCct                                                                                                              <  1:875695/150‑1 (MQ=255)                                            aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTATTTTTAATTTTTTGCCTTCAG‑AGGGAGCGCGGAGGGGGATTGTC‑TTGGctgcacg                                                                                                         >  1:1277435/1‑147 (MQ=255)                                             gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAT‑TTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGATTATC‑ttttcttact                                                                                                          >  1:397854/1‑141 (MQ=255)                                              aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCctgactga                                                                                                        <  2:1257967/150‑1 (MQ=255)                                               aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATCAGAAGAGAGCGGCAAGAGGGGCTATC‑TTTCcttactg                                                                                                         >  2:194603/1‑148 (MQ=255)                                                aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCAGGGAGGGGTTATC‑TTTCCTGACTGGACt                                                                                                     >  2:1022422/1‑151 (MQ=255)                                                 gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTATT‑TTGCctgact                                                                                                          >  1:213572/1‑145 (MQ=255)                                                 gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATT‑TTGCCTGACTGaa                                                                                                       >  1:770292/1‑148 (MQ=255)                                                       gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGAtt                                                                                              <  1:709921/151‑1 (MQ=255)                                                         gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCAT‑TTTTTTTTTTTTTGCCCTTGGAAGGGAGGGGCCGGGGGGATTATCATTCCCTG‑CTGAACTGGCgggt                                                                                              >  2:677557/1‑146 (MQ=255)                                                             cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGttt                                                                                        <  2:941064/151‑1 (MQ=255)                                                              ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCG‑TTTTTTTTTTTCTGGCCAGAGAAGGGAGGGCACAGGGGGATTATG‑TTGTCTTTATTTACTGGCGATTAAgct                                                                                         >  2:2227/1‑147 (MQ=255)                                                                 tataAGCTCCATACCTATACCAGTCAGGAGTGATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                      <  1:233742/149‑1 (MQ=255)                                                                 tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATC‑TTGCCTGACTGAAATGGCGGTTAAGttttc                                                                                      >  1:376714/1‑147 (MQ=255)                                                                   taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCCGAGAGGGGTTTTT‑TTTCCTGACTGAACTGGCGATTAAGTTTCcg                                                                                     >  1:639572/1‑148 (MQ=255)                                                                      gCTCCATACCTGTACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAGGGGTTATT‑TTTCCTGACTGAACTGGCGATTAAGTTTTCGCGGAcc                                                                               >  1:26184/1‑151 (MQ=255)                                                                      gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                               <  1:829233/151‑1 (MQ=255)                                                                          cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGGGGGGTTATT‑TTTCCTGACTGAACTGGGGGTTAAGTTTTcgcggaacc                                                                              >  2:930972/1‑143 (MQ=255)                                                                              ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                       <  1:268028/151‑1 (MQ=255)                                                                                   aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                                                  <  2:822629/151‑1 (MQ=255)                                                                                    ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                                 <  1:488651/151‑1 (MQ=255)                                                                                      aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                               <  1:639250/151‑1 (MQ=255)                                                                                       gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                              <  1:408404/151‑1 (MQ=255)                                                                                        tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCATGCGGAGTGAGCGCCGAGAGGGGTTATT‑TTTCCTGACTGAACCTGCGGTTAAGTTTCCGCCTACCCCGCGTTTTGTTTTTGCg                                                             >  1:663196/1‑151 (MQ=255)                                                                                            gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATC‑TTGCCTGACTGAAATGGCGATTAAGTTTTCGCGTACCCCGCGTTTTGTTATTGCggg                                                           >  1:620451/1‑148 (MQ=255)                                                                                              gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATGTGAAGTGAGGGCGGAGGGGGGTTATC‑TTGCCTGACTTAACTGGCGAGTAAGTTTTCGCGTTCCCCGCGCTTTGTTTTTGCgggt                                                          >  1:55531/1‑146 (MQ=255)                                                                                                 aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCCGCGGTTAGTTATTTCGggtgatg                                                      >  1:683905/1‑149 (MQ=255)                                                                                                    tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCAGAGAGGGATTATC‑TTGTCTGACTGAACTGGGGATTAAGTTTTCCCGTACCCGGGGCTTTGTTATTTCGGGTGATGATAtt                                                 >  1:1159192/1‑151 (MQ=255)                                                                                                    acAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAT‑TTTTTTTTTTTTTGCCCTTAGGAGGGGGGGCCTGGAGGGATTTTC‑TTTGCTGACCTAAATGGCGGTTAAGTTTTCGCGTCACCGGCGGTTTGTTATTTCGGATGGTGAAAtt                                                 >  1:865573/2‑151 (MQ=255)                                                                                                             cATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTTTTTTTCTGCCATGGGAAGGGTGTGCTGGGAGGGATTTTC‑TTGCCTGACTGAATTGGCGGTTAAGTTTCCGCGTACCAGGGGTTTCGGTATTTTGGATGTTGATTTTGTAGTTc                                          >  1:135433/1‑149 (MQ=255)                                                                                                              aTGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                       <  1:338614/151‑1 (MQ=255)                                                                                                                gTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                                     <  1:567635/151‑1 (MQ=255)                                                                                                                   gcgGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAGGGGTTTTC‑TTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGTTTTGTTTTTTCGGGTGATGATATTTAAGCTCGGCTggt                                   >  2:923749/1‑148 (MQ=255)                                                                                                                         ctcGTTTGCGGCGGTCAT‑TTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCCGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                              >  1:1107106/1‑149 (MQ=255)                                                                                                                           cGTTTGCGGCGGTCAT‑TTTTTCTTTTTTTGCCATGAGGAGAGAGCGGGCGGAGGGGTTATC‑TTTGCTTAATGTACTGGCGGTTTAGTTTCCGCGGAACCCGCGTCTTGTTATTGCGGGTGATGTTATTGAAGATCgcgcggttggtgtctt                          >  2:437912/1‑136 (MQ=255)                                                                                                                           cGTTTGCGGCCGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCGGAGAGGGATTATC‑TTGCCTGAGTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTTATTGGTTCa                            <  2:1151825/149‑1 (MQ=255)                                                                                                                            gTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg                         <  1:244919/151‑1 (MQ=255)                                                                                                                               tGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGCGAAGTGAGCGCAGAGGGGGGTTTTC‑TTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTGCGGGTGATGATATTTAAGCTCGGCTGGTTGGGTTATTGt                        >  2:804299/1‑149 (MQ=255)                                                                                                                                gcggcgGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg                      <  2:971657/150‑1 (MQ=255)                                                                                                                                     ggTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                 >  1:296951/1‑150 (MQ=255)                                                                                                                                        taT‑TTTTTCTTTTTCTGCGATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCTTGTCTGAACTGGCGGTTAAGTTTCCGCGTCCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg             >  2:292691/2‑151 (MQ=255)                                                                                                                                        cAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg             >  1:607343/1‑151 (MQ=255)                                                                                                                                        cAT‑TTTTTCTTTTTCTGCCATGAGAAGAGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc              >  2:755754/1‑150 (MQ=255)                                                                                                                                          t‑tttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCTTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCTTTGTCGCCATCAc               >  2:493020/1‑147 (MQ=255)                                                                                                                                               ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCGTCACCGCTgg         >  2:528737/1‑149 (MQ=255)                                                                                                                                                     ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt  <  1:2227/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGT  >  NC_000913/4604201‑4604497

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