Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,344	0	C	T	50.0%	39.6 / 29.3	54	S62F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (6/20);  new base T (26/0);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 3.65e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604195‑4604490                                                                                                                                                      |                                                                                                                                                   aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                   <  2:101999/151‑1 (MQ=255)  cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt                                                                                                                                                  >  1:155218/1‑151 (MQ=255)        gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                                                                            <  1:453858/151‑1 (MQ=255)         aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                             >  2:633970/1‑147 (MQ=255)         aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttt                                                                                                                                               >  2:155949/1‑147 (MQ=255)           aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCtttttttttttttttcc                                                                                                                                          >  1:578232/1‑145 (MQ=255)           aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc                                                                                                                                          >  1:463543/1‑150 (MQ=255)                   cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGaa                                                                                                                                   >  2:901097/1‑149 (MQ=255)                   cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGaa                                                                                                                                   >  1:736033/1‑149 (MQ=255)                      gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAg                                                                                                                                  >  1:927936/1‑147 (MQ=255)                         tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGGgcg                                                                                                                            >  1:552275/1‑150 (MQ=255)                             aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCag                                                                                                                         >  2:172293/1‑149 (MQ=255)                                       agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                             <  2:892125/151‑1 (MQ=255)                                          gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGGCGAGGGGGGTTAtttt                                                                                                           >  2:745261/1‑150 (MQ=255)                                             tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGCGAAGGGAGCGCAGAGGGGGGATATCTTTcc                                                                                                        >  2:841172/1‑150 (MQ=255)                                                tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTGCct                                                                                                       >  1:836584/1‑148 (MQ=255)                                                  aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCATGAGAAGTGGGCGCAGAGAGGGGTTTTTTTGcc                                                                                                        >  1:482184/1‑145 (MQ=255)                                                   gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                 <  1:330976/151‑1 (MQ=255)                                                    aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                   >  1:406362/1‑148 (MQ=255)                                                     aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTTCCTGACTGaa                                                                                                >  2:515807/1‑150 (MQ=255)                                                       gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                                                                             <  1:414988/151‑1 (MQ=255)                                                         ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                            <  1:323335/150‑1 (MQ=255)                                                                 aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATTAGAAGTGAGCGCAGAGAGGGATTTTCTTGCCTGGCTGAAATGGCGa                                                                                         >  1:113662/1‑145 (MQ=255)                                                                     cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCTGCTTTTAAATTTCTGCCATGAGAAGCGTGGGCCGAGAGGGATTATTGTTCCTGACTGAACTGGCGGTTAAGTTTcc                                                                               >  1:961146/1‑151 (MQ=255)                                                                       tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTATGCCGCGGTCTGTTTTTAATTTTTTAACATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                               >  1:732083/1‑149 (MQ=255)                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGATTAAGTTTcc                                                                               >  1:960556/1‑146 (MQ=255)                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTACGCGt                                                                           <  1:640815/150‑1 (MQ=255)                                                                               ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg                                                                      <  2:207895/150‑1 (MQ=255)                                                                                     cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGAGGGCCGAGGGGGATTATTTTGCCTGACTGAACTGGCGATTTAGTTTCCGCGTACCCGGCGttt                                                                >  2:733900/1‑150 (MQ=255)                                                                                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                               <  2:955138/149‑1 (MQ=255)                                                                                           cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                         <  2:960516/151‑1 (MQ=255)                                                                                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTTTTTTTTTGCTCTGAAAAGTGAACGAAGAGAGGGGTTATCTTGTCTGACTTAACACGGGGTTTAGTTTTCGGGTACCCGGccgtttgtt                                                            >  2:743280/1‑139 (MQ=255)                                                                                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTTGCGGTTAAGTTTTCGCGTACCCGGCGTTTTGTtttt                                                         >  1:1006443/1‑149 (MQ=255)                                                                                               cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                      <  1:75714/150‑1 (MQ=255)                                                                                                      cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGGAGTGAGCGCCGAGAGGGATTATTTTTCCTGACTGGACTTGGGGTTTAGTTTTCGCGTACCCGGCGTTTTGTTTTTTCGgatgatg                                               >  2:229982/1‑150 (MQ=255)                                                                                                      cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                               <  1:585269/150‑1 (MQ=255)                                                                                                       aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGGGAAGTGAGCGGAGAGAGGGGGTATCTTGCCTGAGTGAACTGGCGATTAAGTTTTCGCGGACCCCGCGGTTTGTTATTTCGggtg                                                  >  2:81185/1‑146 (MQ=255)                                                                                                       aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGTGGGATTATCTTACCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                             <  1:314276/151‑1 (MQ=255)                                                                                                       aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                             <  1:626015/151‑1 (MQ=255)                                                                                                            aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCCTGGGAAGTGAGCGCAAAGAGGGATTATTTTGCCTGACTGAACTGGCGATTAAATTTCCGCGTACCCGGCGTTTAGTTTTTGCGGATGATGATAtt                                          >  1:584766/1‑149 (MQ=255)                                                                                                               aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                     <  2:553545/151‑1 (MQ=255)                                                                                                               aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                     <  1:387085/151‑1 (MQ=255)                                                                                                               aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                     <  1:912285/151‑1 (MQ=255)                                                                                                                tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATCAGAAGTGAGCGCGGAGAGGGAATATCTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGGGGTTAGTTATTGCGGGTGATGATATTTaa                                       >  1:321694/1‑148 (MQ=255)                                                                                                                   cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGGGAAGGGGGCGCAGAGGGGGATTATCTTTCCTGACTGAAATGGCGGTTAAGTTTTCGCGTACCCGGCGTTTTGGTTTTGCGGGTGGTGATATTGAAGCTCgg                                 >  1:224622/1‑151 (MQ=255)                                                                                                                     tGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCTTAAGGAGGGAGCGCAGAGAGGGATTATTTTTTTTTGCTGAACTTGGGATTTAGCTTCCGCTTACTCCGCGTCTAGTTTTTGAGGGTGATGTTATTGCATCTCTGc                                >  2:280390/1‑150 (MQ=255)                                                                                                                      gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                               <  1:466081/150‑1 (MQ=255)                                                                                                                       tCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCCTGAGAAGGGAGCGCCGGGAGGGGTTATCTTGCCTTACTTAACTGGCGGTTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGGGTGATGATATTTAAGCTCGGCt                               >  2:614741/1‑149 (MQ=255)                                                                                                                         gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAAATGGCGGTTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGGtt                           >  1:877807/1‑151 (MQ=255)                                                                                                                           ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATATTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTGACCGGCGTTTAGTTTTTGCGGGTGATGATATTGAAACTCGGCTGGtt                           >  1:701192/1‑149 (MQ=255)                                                                                                                                      gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg               <  1:133293/150‑1 (MQ=255)                                                                                                                                       cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc             <  2:460193/151‑1 (MQ=255)                                                                                                                                        ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc             >  2:459636/1‑150 (MQ=255)                                                                                                                                         gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt           >  2:503167/1‑151 (MQ=255)                                                                                                                                         gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt           >  2:1003960/1‑151 (MQ=255)                                                                                                                                             tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc       >  1:641276/1‑151 (MQ=255)                                                                                                                                               aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc     <  2:747366/151‑1 (MQ=255)                                                                                                                                                ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt    >  2:1032717/1‑151 (MQ=255)                                                                                                                                                ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt    >  1:865202/1‑151 (MQ=255)                                                                                                                                                  ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  >  1:583494/1‑151 (MQ=255)                                                                                                                                                  ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  >  2:379168/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604195‑4604490

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