Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,555,119	0	A	G	43.9%	77.5 / 30.9	66	I182V (ATC→GTC)	uxuR	fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (10/27);  new base G (29/0);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCAACGTCTGC  >  NC_000913/4554974‑4555269                                                                                                                                                  |                                                                                                                                                        gCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCa                                                                                                                                                     <  1:373480/150‑1 (MQ=255)    tGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAACAAACCAAAGTGGGTACAGtt                                                                                                                                                  >  1:865788/1‑150 (MQ=255)    tGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGGTCCAGtt                                                                                                                                                  >  1:21398/1‑150 (MQ=255)      aaCTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGca                                                                                                                                               <  2:1095010/151‑1 (MQ=255)       aCTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAATCCAATGTGGGTcccggttcac                                                                                                                                              >  1:147473/1‑143 (MQ=255)          ggCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGGTCCAGTTTCACAg                                                                                                                                            >  1:476186/1‑150 (MQ=255)           gCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGcc                                                                                                                                          <  2:4175/151‑1 (MQ=255)              tCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCAGTTGCACAGCCATc                                                                                                                                       >  1:651368/1‑151 (MQ=255)                    gcgcCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATg                                                                                                                                 <  2:262291/151‑1 (MQ=255)                        cGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAATCCAATGTGGGTCCCGGTGCACAGCCCCCTGGAttga                                                                                                                               >  1:783430/1‑146 (MQ=255)                                agcGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATc                                                                                                                     <  2:472674/151‑1 (MQ=255)                                   gAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑Ca                                                                                                                 <  2:811261/151‑1 (MQ=255)                                       gCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAAAGTGGGTCCAGTTGCCCAGCCCTCTGGATGACAAGCTGTATCG‑Caaa                                                                                                               >  2:52885/1‑149 (MQ=255)                                           tGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAAACCAATGTGGAACCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTgg                                                                                                         >  2:895427/1‑151 (MQ=255)                                             aCATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGtt                                                                                                       <  2:1099392/151‑1 (MQ=255)                                             aCATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAATCCAATGTGGATCCAGTTGCACAGCCATCTGGGTGACAGCCTGTATCG‑CAAAGAgggg                                                                                                         >  2:298488/1‑149 (MQ=255)                                               aTGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCCATGTGGGTTCAGTTTCCCAGCCATCTGGGTGACAGCCTTTATTG‑CAAAGAGTgggt                                                                                                       >  1:881484/1‑147 (MQ=255)                                               aTGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAAAGTGGGTCCCGTTGGCCAGCCACCTGGGTGAACGCCTGTATCG‑CAAAGAGGGGGTgg                                                                                                     >  1:790327/1‑151 (MQ=255)                                                 gCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGc                                                                                                   <  2:652352/151‑1 (MQ=255)                                                    gTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATATGGAGCCGGTTGCACAGCCATCTGGATGAGAACCGGTTTCGCCAAAAAGTGTTggggggt                                                                                                 >  1:363429/1‑148 (MQ=255)                                                      tCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAACAATTCAATATGGAGACGGTGGACCCCACATCTCGGTGAAACCCTCTATAG‑CAAAGAGGGGGTGGGCgcccc                                                                                              >  2:373587/1‑147 (MQ=255)                                                       ccATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATca                                                                                              <  2:749350/150‑1 (MQ=255)                                                       ccATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATca                                                                                              <  2:166924/150‑1 (MQ=255)                                                       ccATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAAACAATGTGGGTTCAGTTGCACAACCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGGTGGGGGAGca                                                                                              >  2:1179089/1‑150 (MQ=255)                                                         aTCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAACCAATGTGGGTCCAGTTGCACAGCCATCTGGATGAAAGCCTGTATCG‑CAAAGAGGGGGTGGGCGAAcaca                                                                                            >  2:411134/1‑150 (MQ=255)                                                         aTCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCAGTTGCACAGCCATCTGGGTGACAGCCCGTATCG‑CAAAGAGGGGTTGGGGGGTcaca                                                                                            >  1:425451/1‑150 (MQ=255)                                                           ctcGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGGTTCAACGCCATCTGGATGACAGCCTGTATCG‑CGAAGAGGGGTTGGGGGATCACaaa                                                                                          >  2:696552/1‑150 (MQ=255)                                                             cGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAg                                                                                       <  2:132466/151‑1 (MQ=255)                                                              gCTAGTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGa                                                                                      <  1:977017/151‑1 (MQ=255)                                                                 aTTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAATACAATGGGGATCCAGTTTCACAGCCCTCTGGGTGACAGCCTGTATCG‑CAAAGAGTGGTTGGGGGGTCACAAACAGAAcc                                                                                   >  2:898637/1‑151 (MQ=255)                                                                   tGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAAACCAATTTGGGTCCCGGTGCACAGCCCTCTGGGTGACAGCCTGTATCG‑CAAAGAGGGGTTGGGCGGTCACAAACAGAACCTc                                                                                 >  2:439800/1‑151 (MQ=255)                                                                    gCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCCGTTGCACAGCCATCTGGGTGGCAGCCCGTATCG‑CAAAGAGGGGGTGGGGGATCACAAACAGAAcccc                                                                                 >  2:661012/1‑148 (MQ=255)                                                                       gAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAATCCAATGTGGGTCCCGTTGCACCGCCCTCTGGGTGACAGCCTGTATCG‑CCAAGAGTGGGTGGGCGGTCACAAACAGATCCTcgcc                                                                              >  1:1169059/1‑150 (MQ=255)                                                                          gCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCCGTTGCACAGCCCTCTGGATGACAGCCCGTATCG‑CAAAAAGGGGGTGGGGGATCACAAACAGATCCCCGCcgcg                                                                           >  2:110956/1‑150 (MQ=255)                                                                            aaCGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTa                                                                        <  2:881784/151‑1 (MQ=255)                                                                            aaCGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTa                                                                        <  1:840275/151‑1 (MQ=255)                                                                             aCGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGGTCCAGGTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGGGGTTGGGCGGTCACAAACAAATCCTCGCcgcggt                                                                         >  2:341190/1‑147 (MQ=255)                                                                                cATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCAGGTGCACAGCCCTCTGGGTGACAGCCTGTATTG‑CAAAGAGGGGGTGGGGGATCACAAAAAGATCCTCGCCGCGTTAATc                                                                     >  1:557024/1‑150 (MQ=255)                                                                                   aaCAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAAAACAATGTGGGTCCAGGTGCACAGCCATCTGGGTGACAGCCTGTATTG‑CCAAAAGGGGTTGGGGGATCACAAAAAGATCCTCGCCGCGTTAATCaaa                                                                  >  1:269591/1‑150 (MQ=255)                                                                                   aaCAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAATTCATGTGGATCCACTTGGACCGGCCTCTGGTTTGCAACCTGTATCG‑CAAAGAGGTGGTGGGGGGTGACAAAAAAATCCTCGCCCCGCTaaataaa                                                                  >  2:850523/1‑145 (MQ=255)                                                                                    aCAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCaaaaa                                                                <  1:87765/151‑1 (MQ=255)                                                                                    aCAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCCGGTTCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGGGGTTGGGCGGTCACAAAAAAATCCCCGCCCCGGTTATCaaaa                                                                 >  2:207129/1‑150 (MQ=255)                                                                                     cAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGCAAAAAAAACCATGTGGGTCCAGTTGCACAGCCATCTTGGTGAAAAACTGTATCG‑CCAAGAGGGGGTGGGCGGTCACCAACAAAACCTCCCCGCGGTAATCcaaaa                                                                >  1:467922/1‑150 (MQ=255)                                                                                     cAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGGGGGTCCAGTTGCACAGCCCTCTGGATGACAGCCCGTATTG‑CAAAGAGGGGGTGGGCGATCACAAAAAGATCCTCGCCGCGCTAATCaaaaa                                                                >  2:188600/1‑150 (MQ=255)                                                                                     cAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCCATGTGGGTCCCGGTGCACAGACCTCTGGGTGGCAGCCTGTATCG‑CAAAGAGGGGGTGGGGGATCACAAACAAATCCTCCCCCCGTTAATCaaaaa                                                                >  1:868329/1‑150 (MQ=255)                                                                                      aGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAAAGCAATGTGGGACCAGGTGGACAACCCACTGGAGGACAGCCCGTATTG‑CAAAGAGGGGGTGGGGGATTACAAACAGAACCTCGCCGCGGTAATAaaaaaa                                                               >  2:196961/1‑150 (MQ=255)                                                                                         aTGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATg                                                           <  1:300231/151‑1 (MQ=255)                                                                                           gCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGc                                                          <  1:1119415/150‑1 (MQ=255)                                                                                              ggtggAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGGGGGTCCCGGTGCACAGCCCTCTGGGGGACAGCCCGTATCG‑CAAAGAGGGGGTGGGCGGTCCCAAACAAATCCCCCCCGCGTTAATCaaaaaaaaaggcc                                                        >  1:162504/1‑142 (MQ=255)                                                                                               gtggAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAg                                                     <  1:46072/151‑1 (MQ=255)                                                                                                    gCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATGTGGGTCCCGTTGCACAGCCGTCCGGAGGGCAGCCTGTTTCG‑CAAAAAGAGGTGGGGCGGTCCAAAACAAATCCTCGCCGCGTTAATCAAAAAAAAGGCgcgggcgggg                                                >  1:678727/1‑149 (MQ=255)                                                                                                      tGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCCAAGTGGGTCCCGGTGCACAGCCCTCTGGGTGACAGCCTGTATCG‑CAAAGAGGGGGTGGGCGATCACAAAAAGATCCTCGCCGCGGTAATCAAAAAAAATGCCCGAgcggcg                                                >  1:418325/1‑149 (MQ=255)                                                                                                       gTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGaa                                              <  2:1111561/150‑1 (MQ=255)                                                                                                       gTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCCATGTGGGTACAGTTGCCCAGCCCCCTTGGGGAAAGCCCGTATCC‑CAAAGAGGGGGTGGGCGGTCACAAAAAAAAACTCCCCCCGGTTATTAAAAAAAAAGTCCGCgcgg                                                  >  1:821740/1‑146 (MQ=255)                                                                                                           cGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGAGGCGAAGCTg                                          <  2:796755/150‑1 (MQ=255)                                                                                                             tCAGTCCTGGCAGTGGCGGGAAAAAAAAACAAAGTGGGTCCCGTTGCACAGCCATCTGGGTGGCAGCCCGTATTG‑CAAAGAGGGGGTGGGGGGTCACCAAAAAAACCTCCCCCCGGTAATCAAAAAAAATGCCCGAGCGGCGAAGCTgg                                         >  1:72134/1‑149 (MQ=255)                                                                                                              cAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCa                                       <  1:415325/150‑1 (MQ=255)                                                                                                                gTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAAtg                                    <  1:818353/151‑1 (MQ=255)                                                                                                                               ggAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGaa                     <  1:210982/151‑1 (MQ=255)                                                                                                                                     cAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGtt               <  1:391573/151‑1 (MQ=255)                                                                                                                                         ccAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGc           >  1:707257/1‑151 (MQ=255)                                                                                                                                          cAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCa          <  2:983586/151‑1 (MQ=255)                                                                                                                                           aaTGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCa          >  1:1045316/1‑150 (MQ=255)                                                                                                                                              gtgGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCAACg       <  1:871056/150‑1 (MQ=255)                                                                                                                                                ggATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCAACGTct    >  1:1110039/1‑151 (MQ=255)                                                                                                                                                ggATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCAACGTc     <  2:776088/150‑1 (MQ=255)                                                                                                                                                  aTCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCAACGTctgc  >  2:597002/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        GCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCG‑CAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAAACGTTAAGCAACGTCTGC  >  NC_000913/4554974‑4555269

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