Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	70,434	0	C	T	48.6%	‑6.5 / 115.9	72	N16N (AAC→AAT)	araC	ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (22/15);  new base T (20/15);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

CGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCT  >  NC_000913/70284‑70582                                                                                                                                                       |                                                                                                                                                      cGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAAc                                                                                                                                                       >  2:541180/1‑151 (MQ=255)  ggTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACg                                                                                                                                                      >  2:319931/1‑151 (MQ=255)    ttGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGcc                                                                                                                                                    >  1:259040/1‑151 (MQ=255)     ttggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGcc                                                                                                                                                    <  1:152242/148‑1 (MQ=255)       ggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCAt                                                                                                                                                 >  1:133041/1‑151 (MQ=255)          tAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                                                                                                               <  1:33279/150‑1 (MQ=255)             cGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtg                                                                                                                                           <  1:276429/151‑1 (MQ=255)                   gagCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGt                                                                                                                                     <  1:347702/151‑1 (MQ=255)                    agCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGtt                                                                                                                                    <  1:509245/151‑1 (MQ=255)                      ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                                                                                                   >  2:159976/1‑150 (MQ=255)                       cAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa                                                                                                                                 <  1:579231/151‑1 (MQ=255)                        aGTAAAAGACGCAGTGACGGCAGTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAAc                                                                                                                                >  1:481000/1‑151 (MQ=255)                          tAAAAGACGCAGTGACGGCAATGTCTGATGCCATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGGTTAACGc                                                                                                                              >  1:582986/1‑151 (MQ=255)                            aaaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGcc                                                                                                                             >  2:336377/1‑150 (MQ=255)                               gACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                                                                                                         >  1:98105/1‑151 (MQ=255)                                  gCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                                                                                                        >  1:540049/1‑149 (MQ=255)                                  gCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                                                                                                        >  1:283707/1‑149 (MQ=255)                                       gACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCa                                                                                                                  <  2:458796/150‑1 (MQ=255)                                        aCGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                                                                                                 >  2:281497/1‑150 (MQ=255)                                         cGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                                                                                                               <  2:283700/151‑1 (MQ=255)                                              aTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAt                                                                                                          <  2:621656/151‑1 (MQ=255)                                               tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtc                                                                                                         >  2:257508/1‑151 (MQ=255)                                               tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACCGGTAt                                                                                                          >  2:123803/1‑150 (MQ=255)                                                gTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtc                                                                                                         <  1:558898/150‑1 (MQ=255)                                                      tGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTa                                                                                                           <  1:182667/142‑1 (MQ=255)                                                      tGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGAGGCCAACGGtt                                                                                                            >  2:182662/1‑141 (MQ=255)                                                       gCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑tttt                                                                                                <  2:600014/151‑1 (MQ=255)                                                             tGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATc                                                                                          >  2:493129/1‑151 (MQ=255)                                                             tGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACGGTTATCTCGA‑ttttttttt                                                                                           >  2:126348/1‑148 (MQ=255)                                                                 cAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATcgacc                                                                                      <  1:571093/151‑1 (MQ=255)                                                                   aTTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATcgaacg                                                                                     >  1:626491/1‑150 (MQ=255)                                                                      ggTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCgaccgacc                                                                                  <  2:605873/150‑1 (MQ=255)                                                                         ttcttcTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTg                                                                              >  1:265168/1‑151 (MQ=255)                                                                         ttcttcTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTg                                                                              >  1:117918/1‑151 (MQ=255)                                                                             tctcTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                                                                                                         <  2:248982/105‑1 (MQ=255)                                                                             tctcTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                                                                                                         >  1:248988/1‑105 (MQ=255)                                                                              ctctGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAAt                                                                         >  2:253195/1‑151 (MQ=255)                                                                              ctctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTgggga                                                                          >  2:520088/1‑148 (MQ=255)                                                                                   aaTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAg                                                                    <  2:361418/151‑1 (MQ=255)                                                                                   aaTGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAg                                                                    <  2:34338/151‑1 (MQ=255)                                                                                    aTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTTGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTTATCGACCGACCGCTGGGAATGaaa                                                                     <  1:257515/150‑1 (MQ=255)                                                                                    aTGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAgg                                                                   <  1:33907/151‑1 (MQ=255)                                                                                       gtggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTta                                                                <  2:243306/151‑1 (MQ=255)                                                                                       gtggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTta                                                                <  1:545147/151‑1 (MQ=255)                                                                                          ggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTatat                                                             <  2:126754/151‑1 (MQ=255)                                                                                            aGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATtc                                                           <  1:336384/151‑1 (MQ=255)                                                                                             gTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATtct                                                          <  2:485572/151‑1 (MQ=255)                                                                                              tATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATtctc                                                         >  1:144471/1‑151 (MQ=255)                                                                                               aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCa                                                        <  1:613224/151‑1 (MQ=255)                                                                                                 gAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCaa                                                       <  1:439473/150‑1 (MQ=255)                                                                                                  agaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAtc                                                     >  1:411666/3‑151 (MQ=255)                                                                                                   aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAtct                                                    >  1:107391/1‑151 (MQ=255)                                                                                                    aaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAtctc                                                   >  1:216907/1‑151 (MQ=255)                                                                                                      gTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCAc                                                 >  2:30509/1‑151 (MQ=255)                                                                                                      gTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCAc                                                 >  1:585654/1‑151 (MQ=255)                                                                                                        aTGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCa                                               <  1:440896/151‑1 (MQ=255)                                                                                                          ggCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCAt                                              <  2:171550/150‑1 (MQ=255)                                                                                                            cTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTcg                                           >  1:242784/1‑151 (MQ=255)                                                                                                             tGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTcgc                                          >  1:154504/1‑151 (MQ=255)                                                                                                                         tGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGggtg                              >  2:295267/1‑151 (MQ=255)                                                                                                                             ctccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGgtggtg                           >  2:26513/3‑150 (MQ=255)                                                                                                                             ccccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGTCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGgtggtg                           >  1:296399/1‑150 (MQ=255)                                                                                                                             ccccTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGgtggtg                           >  2:69479/1‑150 (MQ=255)                                                                                                                                  gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAAt                     >  2:151653/1‑151 (MQ=255)                                                                                                                                    tgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATc                    >  2:623154/1‑150 (MQ=255)                                                                                                                                     gcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCa                   <  2:48070/150‑1 (MQ=255)                                                                                                                                        gggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAggg                >  1:448798/1‑150 (MQ=255)                                                                                                                                        gggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAggg                >  1:139437/1‑150 (MQ=255)                                                                                                                                           aTACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACga            <  1:30509/151‑1 (MQ=255)                                                                                                                                             aCTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACg             >  1:458455/1‑148 (MQ=255)                                                                                                                                                 gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAttt      >  2:17737/1‑151 (MQ=255)                                                                                                                                                 gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAttt      >  1:192373/1‑151 (MQ=255)                                                                                                                                                     aaCGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCt  >  2:274286/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGA‑TTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCT  >  NC_000913/70284‑70582

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