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

ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604195‑4604490                                                                                                                                                      |                                                                                                                                                     aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                     <  2:801833/151‑1 (MQ=255)    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                                    <  2:154995/149‑1 (MQ=255)       aGAAAAATGCTTCATGTGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTTCGGCGGTCATTttttctttttct                                                                                                                                               <  2:492133/151‑1 (MQ=255)        gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                                 >  1:321863/1‑148 (MQ=255)           aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc                                                                                                                                            >  1:534578/1‑150 (MQ=255)             aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCTAt                                                                                                                                          >  2:611311/1‑150 (MQ=255)             aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                                                                                         <  1:601464/151‑1 (MQ=255)               gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATgtg                                                                                                                                       >  1:15336/1‑149 (MQ=255)                  tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATTCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                                     <  2:792551/150‑1 (MQ=255)                   cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGAGGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcctcg                                                                                                                                         >  2:119860/1‑142 (MQ=255)                    aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATTCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                                  <  2:777005/151‑1 (MQ=255)                    aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgagga                                                                                                                                     >  2:265246/1‑146 (MQ=255)                          agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAg                                                                                                                                >  1:586837/1‑147 (MQ=255)                           gcaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                           <  1:459/149‑1 (MQ=255)                             aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                         <  1:352018/151‑1 (MQ=255)                             aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                          <  1:391709/150‑1 (MQ=255)                             aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                          <  1:8983/150‑1 (MQ=255)                               cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagag                                                                                                                         >  1:330755/1‑149 (MQ=255)                               cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                                       <  2:683930/151‑1 (MQ=255)                                ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                      >  1:608147/1‑151 (MQ=255)                                ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                      <  2:156961/151‑1 (MQ=255)                                 gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAggg                                                                                                                     >  2:420983/1‑151 (MQ=255)                                   aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGGGAAGGGAGGGCAGAGAggg                                                                                                                     >  1:851468/1‑149 (MQ=255)                                        gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGGGCAGAGAGGGGTTATCt                                                                                                              >  1:442498/1‑151 (MQ=255)                                         agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCtt                                                                                                             >  1:532226/1‑151 (MQ=255)                                         agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                             <  1:323100/151‑1 (MQ=255)                                              aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                         <  2:150270/150‑1 (MQ=255)                                               ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGGAGTGGGCGGCGAGGGGGGTTATCTTTCct                                                                                                         >  1:645538/1‑149 (MQ=255)                                                 cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                      <  2:861992/150‑1 (MQ=255)                                                 cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCAGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                     >  2:417906/1‑151 (MQ=255)                                                   gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGAGAGCGCAGAGGGGGGTTATCTTTCctgcat                                                                                                     >  1:272179/1‑146 (MQ=255)                                                   gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                     >  1:15878/1‑149 (MQ=255)                                                   gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGAGGGATTATCTTGCctgactg                                                                                                    >  1:377010/1‑150 (MQ=255)                                                    aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGaa                                                                                                  <  2:560961/151‑1 (MQ=255)                                                         ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTCTGCCCAGAGAAGTGAGAGCTGAGAGGGGTTATCTTGCCTGTCTGAACttg                                                                                              >  2:687042/1‑148 (MQ=255)                                                         ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTTCCATGAGAAGGTAGCGCAGGGAGGGGTTATTTTGCctgtctga                                                                                                   >  2:367192/1‑145 (MQ=255)                                                         ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGAGCAGGGAGGGAGTATCTTGCCTGGCTTAACt                                                                                                >  2:313420/1‑148 (MQ=255)                                                               gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTaa                                                                                     <  1:402525/151‑1 (MQ=255)                                                                    ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGATTTTCTTGCCTGACTGAACT‑GGC‑GATTAAGttt                                                                                 >  1:317776/1‑150 (MQ=255)                                                                     cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTGCCTGAGGAGAGAGCGCGCAGAAGGAGTATCTTTTCTGATTGAACT‑CGC‑GCTTAAGtttt                                                                                >  2:204599/1‑149 (MQ=255)                                                                      ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCATGTGGAGTGAGGGCAGAGAGGGGTTATTTTTCCTGACTGAAAT‑GGC‑GAtt                                                                                       >  2:135581/1‑142 (MQ=255)                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGTGAAGTGAGCGCGGAGAGGGGTTATCTTGCCTGACTGTACT‑GGC‑GACTAAGttttcacggaac                                                                        >  2:231828/1‑140 (MQ=255)                                                                                cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGGAATTTTTTTTTTTCTTCCCTAGGGAGGGAGGGCAGGGGGGGCGTATCTTGGCTTGCTTAAATGGGG‑GATTAAGTTTCCCGGAACCCg                                                                      >  2:722762/1‑150 (MQ=255)                                                                                   aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTTCCTGACTGAAAT‑GGC‑GATTAAGTTTTCGCGTACCCGGCGtt                                                                 >  1:690490/1‑151 (MQ=255)                                                                                     cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGAGAGCGCAGAGGGGGGTTATTTTTCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGGACCCCGCGttt                                                                >  1:100590/1‑150 (MQ=255)                                                                                        taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATTTTTCCCGACTGAACT‑GGC‑GATTAAGTTTTCGCGTACCCGGCGTTTAGtt                                                            >  1:388733/1‑151 (MQ=255)                                                                                                aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                                    <  1:628320/151‑1 (MQ=255)                                                                                                aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCACGAGAAGCGCGCGCAGCGGGGGGGTTTCTTGTCTCACTGGACT‑GGC‑GGTTAAGCTTCCGCGGAACCGGCGCTTTGTTATTGTggg                                                    >  2:126040/1‑150 (MQ=255)                                                                                                 ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCTTGGTGAGATGGAGCAGGGGGGGGTTTTTTTGTCTGCCCTGACT‑GGC‑GATTTAGGTTCTGCGTACCCCGCGTTTTattatt                                                         >  2:663182/1‑145 (MQ=255)                                                                                                  gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAATTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                  <  2:786114/151‑1 (MQ=255)                                                                                                   aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                                 <  2:564650/151‑1 (MQ=255)                                                                                                          tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCCGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATAtt                                          >  1:835577/1‑151 (MQ=255)                                                                                                           cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                         <  1:12095/151‑1 (MQ=255)                                                                                                            aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCATGAGGAGTGGGCGCAGGGGGGGGTTTTCTTTCCTGAATTAAAT‑GGC‑GGTTAAATTTCCGCGGGTCCGGGGGTTTATTTTTGGGGATTATGATAtt                                          >  2:80472/1‑149 (MQ=255)                                                                                                              gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                      <  1:836489/151‑1 (MQ=255)                                                                                                               aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                     <  1:646424/151‑1 (MQ=255)                                                                                                                tGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGCGAGCGGCGGGAGGGGTTATCTTTCCTTACTGAACT‑GGG‑GAATAATTTTTCCCGTGCCCGGGGGTTTGGTTTTGCGggtggaggttat                                          >  2:344418/1‑137 (MQ=255)                                                                                                                        cgcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCCGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGGTGATGATATTTAAGCTCGGCTg                              >  1:579008/1‑149 (MQ=255)                                                                                                                           ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGGAGTGAGCGCAGGGAGGGATTTTCTTGTCTGTCTGAACT‑GGCGGTTTAAGTTT‑CGCGGACCCGGGGTTTAGGTTTTGCGGGTGGTGATATTGAAGCTCGGCTGATTgg                         >  2:271704/1‑151 (MQ=255)                                                                                                                                  gTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GGTTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTTAAGCTCGGCTGATTGGTTCAtt                   >  1:82986/1‑150 (MQ=255)                                                                                                                                     tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGGGCCGCGGGGGGTTATCTTTTCTTACTGAACT‑TGC‑GATTAAGGTTTCCCGTACCCGGCGGTTTGTTTTTGCGGATGGTGATATTTTAGCTCGGCTTAATTGTTCATTGt                 >  2:774382/1‑149 (MQ=255)                                                                                                                                      gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg               <  2:645383/150‑1 (MQ=255)                                                                                                                                       cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc             <  1:265321/151‑1 (MQ=255)                                                                                                                                           ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa         <  1:628709/151‑1 (MQ=255)                                                                                                                                                 tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg   >  2:99635/1‑151 (MQ=255)                                                                                                                                                  ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  >  1:213511/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑