Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	40.7%	100.9 / 12.3	59	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (7/28);  new base G (24/0);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366384‑366674                                                                                                                                                |                                                                                                                                                     ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGgttttt                                                                                                                                              >  2:1036384/1‑150 (MQ=255) ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttc                                                                                                                                              <  1:397714/151‑1 (MQ=255)                 gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGGGGCAGGGGGGGTTTTCTTTTTACCAGGgag                                                                                                                                >  1:602771/1‑149 (MQ=255)                   gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCTGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                            <  2:571383/151‑1 (MQ=255)                         tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGaa                                                                                                                        >  2:648813/1‑149 (MQ=255)                           acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGC‑TTGGGGGGCCGGGGGGGTTTTCTTTTTACCCGGGAGGCCGGCCAc                                                                                                                       >  1:247317/1‑148 (MQ=255)                             aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCaa                                                                                                                        >  1:960596/1‑145 (MQ=255)                             aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGACTTGGCGGCCGGGGGGTTTTTTCTTTTC‑CCAGGGGGGCGGGGAACAACTg                                                                                                                  >  1:889063/1‑151 (MQ=255)                                     cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                          <  1:888474/151‑1 (MQ=255)                                        tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGG‑TTGGGCGCCCGGGGGGTTTTTCTTTTTACCAGTGAGGCGGGGCACAGCTGGTTGccc                                                                                                           >  1:357041/1‑147 (MQ=255)                                           gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCg                                                                                                    <  1:223069/151‑1 (MQ=255)                                           gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGCGGCCGGGGGGGTTTTCTTTTCACCCGTGAGAAGGGCCACAGCTGGTTGCCCTTCAcc                                                                                                     >  1:705353/1‑150 (MQ=255)                                           gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTA‑TTTGGCGCCCGGGGGGTTTTTTTTTTTCCCAGTGAGACGGGCAACAACTGATTGCCCTTTAcc                                                                                                     >  2:229466/1‑150 (MQ=255)                                               aCTGCCCGCTTTCCAGTCGGGAACCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc                                                                                                  <  1:192307/149‑1 (MQ=255)                                                 tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTT‑TTTGGCGGCAGGGGGGGTTTTTTTTTTACCAAGGGGACGGGCCACAGCTGGTTTTCCTTCTACGcc                                                                                                  >  2:385788/1‑147 (MQ=255)                                                  gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                                                                                              <  1:391280/150‑1 (MQ=255)                                                  gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTGGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCGGGGGACGGGCAACAGCTGGTTGGCCtttcccccc                                                                                                  >  2:370675/1‑139 (MQ=255)                                                          tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGCG‑TTGGGCGGCCGGGGGGGTTTTCTTTTTCCCCGGGAGGCGGGCACCAGCCGGCTTGCCTTCACCGCCTGGCCCtcggag                                                                                      >  1:786794/1‑145 (MQ=255)                                                               tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTTGGGCCCGGGGGGTTTTTTTTTTTTCCCCGGAGACGGGCACCAGCTGCTTGCCCCTCACCCCCCTGccccgcgggggtt                                                                                  >  1:697714/1‑139 (MQ=255)                                                                cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                               <  2:843646/151‑1 (MQ=255)                                                                   gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGGGCCCGGGGGGGTTTTCTTTTCACCAGAGAGACGGGCAACAGCTGATTGCCCCTCACCCCCTGGCCCTGAGAGGGGTgca                                                                               >  1:53293/1‑148 (MQ=255)                                                                     aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCCTT‑TTTGGGGCCCCGGGGGGTTTTCTTTTTACCAAGGAGACCGGCGACAGCTCAATTCCCTTCTCCGCCTGGGCCCGCGAGAGtt                                                                                  >  2:312120/1‑143 (MQ=255)                                                                        cTGTCGGGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGTA‑TTTGGGGCCCGGGGGGTTTTTTTTTTTACCACGGAGAAGGGCAACAGCTGAGTTCCCCTCCCCCCCAGGCCCTcggagaagttccgcaagagg                                                                       >  1:571162/1‑137 (MQ=255)                                                                          gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCTTT‑GTTGGGGCCACGGGGGGGTTTCTTTTTACCAGCGGGAGGGGGCACACATGATTGGCCCTCTCCCCCTTGCCCTGCGAGGGTGTCCGCAAGCg                                                                        >  1:462013/1‑148 (MQ=255)                                                                          gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                                       <  2:919213/149‑1 (MQ=255)                                                                             gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGGCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGGGAGAGGTGCAGCAAGCGGTcc                                                                    >  2:145958/1‑149 (MQ=255)                                                                             gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCAc                                                                  <  1:800480/151‑1 (MQ=255)                                                                               gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACg                                                                 <  2:900417/150‑1 (MQ=255)                                                                                   gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                            <  1:791316/151‑1 (MQ=255)                                                                                      gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtt                                                           <  2:717964/149‑1 (MQ=255)                                                                                          tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTTACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGTCCACGCTGGTTTGgccc                                                     >  2:613700/1‑151 (MQ=255)                                                                                          tAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGTA‑TTGGGCTCCCGGGGGGTTTTTCTTTTTTCCAATGTGACGGGCAACACCCTGTTTCCCTTCACCGCCTTGCCCCCGGAGAGGTGCAAGAAGAGGTCCCCGGTGGTTTGccc                                                      >  2:807469/1‑150 (MQ=255)                                                                                           aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCCCGCTGGGTTTcccc                                                     >  2:303153/1‑150 (MQ=255)                                                                                              gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCgcagca                                                 <  1:635322/151‑1 (MQ=255)                                                                                                  cGGCCAACGCGCGGGGAGAGGCGGTTTGGGGA‑TTTGGGGGCCGGGGGGTTTTTCTTTTTACCAGGGAGACGGGGAACAACCGGTTGCCCCTCCCCCCCTTGCCCCGGGAGAGGTGCAGCAAAGGGGCCACCCTGGTTTTCCCCCGCCgg                                               >  1:780407/1‑149 (MQ=255)                                                                                                    gCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCCCGGGGGGGTTTCTTTTCACCACTGAGAAGGGGCACAACTGATTGCCCCTCCCCGCCCGGCCCCGAGAGGGGTGCCCCAAACGGGCCACGCTGGTTTGcccccac                                                  >  1:362613/1‑141 (MQ=255)                                                                                                    gCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaa                                           <  1:189400/151‑1 (MQ=255)                                                                                                     ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                          <  2:627411/151‑1 (MQ=255)                                                                                                       aaCGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                        <  1:149006/151‑1 (MQ=255)                                                                                                          gcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                     <  2:180242/151‑1 (MQ=255)                                                                                                           cgcgGGGAGAGGCGGTTGGCGTA‑TTGGGCGCCAGGGTGGTTTTTCGTGGCACCAGTGAGACGGGCAACAGCGGAGGGCCCGGCACCGCCGGGCCCGGAGAGAGGGGCAGCAAGCGGTCCACGCTGGTGTGCCCCAGCAGGCGaaaagcat                                     <  1:556410/150‑5 (MQ=255)                                                                                                                ggAGAGGCGGGTTTGGTA‑TTTGGCGCCAGGGGGGTTTTTTTTTTTACCAGGGAGACGGGCAACAGCTGATTGCCCTTCCCCCCCTTGCCCTGAGAGAGGTGCAGCAAGCGGGCCACGCTGGTTTTCCCCCGCAGGCGAAAAACCTGTTTg                                >  1:251772/1‑150 (MQ=255)                                                                                                                 gagaGGCGGTTTGCGTA‑TTGGGCGCCAGGTTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                               <  1:827216/150‑1 (MQ=255)                                                                                                                     ggCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                          <  1:106093/151‑1 (MQ=255)                                                                                                                       cGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                         >  1:30973/1‑150 (MQ=255)                                                                                                                        ggTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                       <  2:705656/151‑1 (MQ=255)                                                                                                                         gTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      <  2:979704/151‑1 (MQ=255)                                                                                                                          tttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGCGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCGCCAGCAGGTGAGAATCGTGTTTGATGGTGGTTa                      >  1:206918/1‑150 (MQ=255)                                                                                                                            tGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATACTGTTTGATGGTGGTTAAcg                   >  1:566138/1‑151 (MQ=255)                                                                                                                              cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                 >  1:963621/1‑151 (MQ=255)                                                                                                                              cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                 <  1:967058/151‑1 (MQ=255)                                                                                                                              cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGATGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  >  1:550740/1‑150 (MQ=255)                                                                                                                                    tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              >  1:849228/1‑149 (MQ=255)                                                                                                                                     gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat            <  2:29513/150‑1 (MQ=255)                                                                                                                                       gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata         <  1:28412/151‑1 (MQ=255)                                                                                                                                          ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc       >  1:626534/1‑150 (MQ=255)                                                                                                                                          ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa      <  2:670695/151‑1 (MQ=255)                                                                                                                                          ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa      <  2:335544/151‑1 (MQ=255)                                                                                                                                             gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa   <  2:528683/151‑1 (MQ=255)                                                                                                                                              ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  1:427504/151‑1 (MQ=255)                                                                                                                                                tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAACTGAg  >  1:867010/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     GGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366384‑366674

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