Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,344	0	C	T	40.0%	66.3 / 30.9	55	S62F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (3/30);  new base T (22/0);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.77e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604197‑4604490                                                                                                                                                    |                                                                                                                                                     tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAtttttttttt                                                                                                                                                   >  1:854429/1‑151 (MQ=255) tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttt                                                                                                                                                   >  2:769280/1‑151 (MQ=255)   aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttt                                                                                                                                                 >  1:1149885/1‑151 (MQ=255)   aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                                 <  1:1217273/151‑1 (MQ=255)      gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                                                                              <  2:970197/151‑1 (MQ=255)       aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCtttttttttttt                                                                                                                                                 >  1:641237/1‑147 (MQ=255)       aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                                                                              <  1:596074/150‑1 (MQ=255)        aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                                                                            <  2:1052578/151‑1 (MQ=255)          aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCTttttttttttttt                                                                                                                                                >  1:203150/1‑144 (MQ=255)             gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTTCCa                                                                                                                                           >  2:797369/1‑147 (MQ=255)              cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGcc                                                                                                                                            >  2:509559/1‑145 (MQ=255)              cGTCATGAGTAGTACCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                      <  1:168777/151‑1 (MQ=255)                 cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                   <  1:718963/151‑1 (MQ=255)                  aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCCTGAGAAgg                                                                                                                                   >  1:33445/1‑149 (MQ=255)                   tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgttaa                                                                                                                                     >  2:638269/1‑143 (MQ=255)                           aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCag                                                                                                                           >  2:144999/1‑149 (MQ=255)                                 aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                                                                                    <  2:591105/150‑1 (MQ=255)                                   tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAtt                                                                                                                  <  1:540995/150‑1 (MQ=255)                                    gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTAt                                                                                                                >  2:1197878/1‑151 (MQ=255)                                      gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                                                                                              <  1:772557/151‑1 (MQ=255)                                       agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTCTTTTTTTGCCATGAGAAGGGGGGGCAGAGAGGGAGTAtttt                                                                                                             >  2:1214512/1‑151 (MQ=255)                                        gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                            <  2:82882/151‑1 (MQ=255)                                          ttATTCAGAAAGTTTGCGGGTAACCCTTAGAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                           <  2:690571/150‑1 (MQ=255)                                              tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctga                                                                                                       <  1:435681/150‑1 (MQ=255)                                                 gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTACATGAGAAGGGAGGGCAGAGAGGGGTTTTCTTTCctgaat                                                                                                     >  2:540833/1‑147 (MQ=255)                                                       ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                              <  1:1018691/150‑1 (MQ=255)                                                       ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                              <  1:1114582/150‑1 (MQ=255)                                                           gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGTGAGGGCGGAGAGGGGTTATTTTGCCTGACTGAACTGGCGG‑tt                                                                                        >  1:656850/1‑151 (MQ=255)                                                               aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGAGGGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCGA‑TTaa                                                                                      >  2:974305/1‑149 (MQ=255)                                                                aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGGTATCTTTCCTGACTGAACTGGGGA‑TTAAGtt                                                                                   >  1:529461/1‑151 (MQ=255)                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCcg‑cg                                                                            <  1:13655/151‑1 (MQ=255)                                                                        aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGGGGGGGCAGAGAGGGGTTATTTTTCCCGACTGAAATGGCGG‑TTAAGttttc                                                                                >  2:503651/1‑144 (MQ=255)                                                                         aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTAc                                                                         <  1:359370/151‑1 (MQ=255)                                                                             ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCg                                                                      <  2:776073/150‑1 (MQ=255)                                                                               aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCg                                                                   <  1:211424/151‑1 (MQ=255)                                                                                tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGGAGGGAGGGAATTTTTTTCCTGACTGAAATGGCGG‑TTAAAGTTTCG‑CGTAccccg                                                                     >  2:1160208/1‑146 (MQ=255)                                                                                       aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGGAGGGGGGGGCGAGAGGGATTATCTTGCCTGACTGCAACGGCGA‑TTAAGTTTCCGCCGAACCCgggttt                                                                 >  2:716645/1‑142 (MQ=255)                                                                                        ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGtt                                                            <  1:509544/149‑1 (MQ=255)                                                                                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTGTTTTTTTTTTTTCTTGCTTGCGAAGTGAGGCCCGGGAGGGGTTTTCTTTTCTTACTGAACTGGCGACTTAATTTTCC‑‑CGCTCCCCGCGTTTAgctttt                                                         >  1:938650/1‑146 (MQ=255)                                                                                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCg                                                      <  2:1205055/150‑1 (MQ=255)                                                                                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCg                                                      <  2:180479/150‑1 (MQ=255)                                                                                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGTTTTTTTGCCTGACTGAACTGGCGA‑TTAAGTTTTCG‑CGTACCCCGCGGTTTGTTTTTTCgg                                                     >  2:684584/1‑149 (MQ=255)                                                                                                    cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTCCCCTGAGAAATGAGCGCAGAGAGGGGTTATTTTGTCTGACTGAACTGGCGG‑ATAAGTTTTCG‑CGGACCCCGCGTTTTGTTTTTTCg                                                      >  1:745674/1‑143 (MQ=255)                                                                                                     aTTTCAGGATGTCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAGCTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTGTAGTTATTGCGgatgatg                                               <  1:197067/149‑1 (MQ=255)                                                                                                     aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTGCCTTACTGAAATGGCGA‑TTAAGTTTTCG‑CGTACCCGGCGTTTTGTTTTTGCGgatgatg                                               >  2:988427/1‑149 (MQ=255)                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGGGGAGGGAGGGATTATTTTGCCTGACTGAAATGGCGG‑TTTAGTTTTCG‑CGTACCCGGCGGTTTGTTATTGCGGGTGATGGTAtt                                          >  2:85617/1‑150 (MQ=255)                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                          <  2:634172/150‑1 (MQ=255)                                                                                                           ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTTCGGATGATGATATTGaa                                       <  1:872431/151‑1 (MQ=255)                                                                                                                 cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGGGGGCAGAGAGGGGGTATCTTTCCCGAGTTAAATGGGGA‑TTTAGTTTTCC‑CGGACCCCGCGTTTTGTTTTTTCGGGTGATGATATTTAAGGTCgg                                 >  2:265871/1‑151 (MQ=255)                                                                                                                    gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                               <  1:39106/150‑1 (MQ=255)                                                                                                                    gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                              <  1:1082102/151‑1 (MQ=255)                                                                                                                       gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGCAGAGGAGGATTATTTTTCCTGAATTTACTGGCGG‑TTAATTTTTCC‑CGTACCCCGCGTTTTGTTTTTTCGGGTGTTGTTATTTAAGCTCGGCTGGtt                           >  2:1018502/1‑151 (MQ=255)                                                                                                                          gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGAGAGCGCAGAGAGGGATTATATTGTCTGATTGAACTGGTGG‑TTAAGTTTCCG‑CGGACCCGGCGGTTTGTTTTTTCGGGTGATGATATTTAAGCTCGGCTGAtt                           >  1:71796/1‑148 (MQ=255)                                                                                                                              tcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGCGCGCCGGGGGGGGTTATTTTTCCTGACTGAAATGGCGA‑GTAAGTTTTCC‑TGTAACCGGCGTTTTGTTATTTTGGGTGGTGGTATTGAAGCCCGGCTGATTTGtt                       >  1:712210/1‑148 (MQ=255)                                                                                                                               cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                    <  1:692802/150‑1 (MQ=255)                                                                                                                                gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg                  <  1:431556/151‑1 (MQ=255)                                                                                                                                   tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                <  1:38731/150‑1 (MQ=255)                                                                                                                                    gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc              <  2:1113149/151‑1 (MQ=255)                                                                                                                                            cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg      <  1:265011/151‑1 (MQ=255)                                                                                                                                             aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc     >  1:925190/1‑151 (MQ=255)                                                                                                                                              ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt    >  1:690636/1‑151 (MQ=255)                                                                                                                                                ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  >  2:221220/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCG‑CGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604197‑4604490

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