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%	68.8 / 27.5	55	S62F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base T (22/0);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.15e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604196‑4604494                                                                                                                                                     |                                                                                                                                                        cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt                                                                                                                                                       >  1:828905/1‑151 (MQ=255) cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATtttttttt                                                                                                                                                       >  1:1542/1‑151 (MQ=255)   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                                     <  2:653168/151‑1 (MQ=255)    aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttt                                                                                                                                                    >  1:338348/1‑151 (MQ=255)       gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                                                                                   >  1:1242275/1‑148 (MQ=255)         aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttttttc                                                                                                                                               >  2:1126262/1‑146 (MQ=255)               cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATGCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                         <  1:129798/151‑1 (MQ=255)                  cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                      <  2:1065293/151‑1 (MQ=255)                     gagaAGTATCGGTATTGAGCGTTTATTAAGACAGTTTGCGGGTAACCCTTATAAGCCCCATACCTATACCAGTCAGCAATCATTTGAGGATCCCATGTCGCGGATCTCGTTTGCGTCGGGGATTTTGTTTTTTTTTgtaatgcgga                                                                                                                                        >  2:720729/1‑137 (MQ=255)                       gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                                 <  1:858599/151‑1 (MQ=255)                            aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                            <  1:627755/151‑1 (MQ=255)                                   ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCCGGGAGGGGtt                                                                                                                     >  1:1164349/1‑151 (MQ=255)                                    tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAtt                                                                                                                     <  2:446651/150‑1 (MQ=255)                                      agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGGGGGGGGTTAt                                                                                                                   >  1:1128596/1‑150 (MQ=255)                                      agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                  <  2:297007/151‑1 (MQ=255)                                            tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCCGGGGGGGGTTATCtt                                                                                                                >  1:789648/1‑147 (MQ=255)                                             aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                            <  2:990582/150‑1 (MQ=255)                                               tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGAAGGGAGCGCCGAGGGGGGTTATTTTTCctgc                                                                                                          >  1:116469/1‑149 (MQ=255)                                               tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGcc                                                                                                             >  1:313803/1‑147 (MQ=255)                                                 aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGTCATGAGAAGTGAGGGCAGAGAGGGGTTAttct                                                                                                                >  1:877562/1‑139 (MQ=255)                                                  gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGGCAGGAGTCATTTCAGGATGCCATGTCGCGGATCGCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGAGATTATCTTGCctgccgga                                                                                                      <  2:175749/151‑6 (MQ=255)                                                   aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCAAGGGGAGGGAGAGCAGGGAGGGGGTTTCTTGGccgaccg                                                                                                       >  1:193685/1‑147 (MQ=255)                                                    aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTCGCCATCTGATGAGTGGGGAGAGAGGGGTTATGttttctcttt                                                                                                        >  2:882759/1‑139 (MQ=255)                                                     aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                                                                                   <  2:432091/151‑1 (MQ=255)                                                            gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGTGAAGTGAGCGCGGAGAGGGGTTATATTTTCTGACTGGACTGTGGA‑tt                                                                                           >  2:681134/1‑151 (MQ=255)                                                               tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTTTTTTTCTTCCATGAGAAATGAGGGGCGAGAGGGATTATCTTGCCTGACTGAACTGGGGA‑ttta                                                                                         >  1:701274/1‑148 (MQ=255)                                                                aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAACTGGCGG‑TTaa                                                                                         >  1:376303/1‑149 (MQ=255)                                                                  cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTTCCTGAGAAATGAGCGCCGGGAGGGGGTATTTTTCCTGTCTGAACTGGGGA‑TTAAGttt                                                                                     >  2:999363/1‑151 (MQ=255)                                                                      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGGAGAGGGGGGTTATCTTGCCCGACTGAACTGGGGG‑TTAAGttt                                                                                     >  1:113051/1‑147 (MQ=255)                                                                          aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATTGGGAGAGAGCCCGGAGGGGGGGTTTTTTTCCTGGCTGAACTTGGGG‑TTTAGTTTTcgcggtc                                                                             >  2:256552/1‑148 (MQ=255)                                                                          aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATCAGAAGAGAGGGGGGAGAGGGGTTATCTTGGCCGACTTAACTGGCGAGTTAAGTTTcgcg                                                                                 >  1:636196/1‑145 (MQ=255)                                                                              ccATACCTATACCAGTCAGGAGTCATTCCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCTATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCg                                                                          <  1:1085258/150‑1 (MQ=255)                                                                               cATCCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCgg                                                                         <  1:551151/150‑1 (MQ=255)                                                                                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAATGGCGA‑TTAAGTTTCCGCGTACCCGGc                                                                        >  1:386675/1‑150 (MQ=255)                                                                                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCg                                                                       <  1:1143425/151‑1 (MQ=255)                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGtt                                                                     <  1:1262454/151‑1 (MQ=255)                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGtt                                                                     <  1:150712/151‑1 (MQ=255)                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGtt                                                                     <  2:828875/151‑1 (MQ=255)                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGt                                                                      <  2:306870/150‑1 (MQ=255)                                                                                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGGGGGCGCCGAGGGGGGTTATCTTTCCTTACTGAACTGGCGA‑TTAAGTTTTCCCGGTCCCGGCGGTTTGt                                                                 >  2:1249819/1‑151 (MQ=255)                                                                                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                 <  2:1163779/151‑1 (MQ=255)                                                                                       taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGGGAGGGATTTTTTTGCCTGACTGAACTGGGGA‑TTAAGGTTTCGCGGACCCGGCGGTTAg                                                                  >  2:565395/1‑149 (MQ=255)                                                                                        aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                <  1:714147/150‑1 (MQ=255)                                                                                         ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                              <  1:228638/151‑1 (MQ=255)                                                                                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                                   <  2:671507/150‑1 (MQ=255)                                                                                                      aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCCTTAGAAGTGGGGGACGAGAGGGGTTATCTTGCCTGCCTGAAATGGGGAGTTAGGTTTCC‑CGGCACCGGCGGTTTGTTATTGGGggtgatg                                                   >  2:78823/1‑149 (MQ=255)                                                                                                      aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGTGAAGTGAGCGCCGGGAGGGGTTTTCTTGCCTTTCTGAACTTGCGG‑TTAAGGTTTCTCGTACCCCGGGGTTTGTTATTGCGggttatggt                                                 >  2:485023/1‑149 (MQ=255)                                                                                                         tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                               <  1:732291/150‑1 (MQ=255)                                                                                                           aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGGGGGGTTTTTTTTTCTGGCTGAACTGGCGG‑TTTAGGTTTCCCGGACCCCGCGTTTTGGTTTTGCGGGTGGTGATAtt                                              >  2:429128/1‑149 (MQ=255)                                                                                                              aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                         <  1:1003968/151‑1 (MQ=255)                                                                                                                gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                        <  2:989778/150‑1 (MQ=255)                                                                                                                   aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACTGAACTGGGGA‑TTTAATTTTCGCGTACCCCGCGTTTTGTTATTGGGGGTGGTGATATTTAAGCTCgg                                     >  1:194157/1‑150 (MQ=255)                                                                                                                     gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGGGGGCAGGGAGGGGTTATCTTTTCTGACTGAAATTGCGG‑TTTAGTTTTCCCGGACCCCGCGTTTTGGTATTTTGGGGGATGATATTGAAGGTccg                                     >  2:458635/1‑146 (MQ=255)                                                                                                                     gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                                  <  1:632989/151‑1 (MQ=255)                                                                                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                               <  1:438204/151‑1 (MQ=255)                                                                                                                          ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  2:769241/151‑1 (MQ=255)                                                                                                                                    tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                    <  1:809182/150‑1 (MQ=255)                                                                                                                                        gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt               <  1:510552/151‑1 (MQ=255)                                                                                                                                          ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              >  1:1013921/1‑150 (MQ=255)                                                                                                                                            tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc           >  1:28416/1‑151 (MQ=255)                                                                                                                                             cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc           <  2:1036463/150‑1 (MQ=255)                                                                                                                                                tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg       >  1:1134815/1‑151 (MQ=255)                                                                                                                                                  tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg      >  1:807544/1‑150 (MQ=255)                                                                                                                                                  tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg      >  2:908371/1‑150 (MQ=255)                                                                                                                                                  tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg      <  1:1237640/150‑1 (MQ=255)                                                                                                                                                   ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc    >  1:588986/1‑151 (MQ=255)                                                                                                                                                    tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc    <  1:831941/150‑1 (MQ=255)                                                                                                                                                    tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGGTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGAGGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc    <  2:1005607/150‑1 (MQ=255)                                                                                                                                                     ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGGTTCCGCGTACCCGACGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg  >  1:1065544/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                        CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604196‑4604494

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: ‑