Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,844	0	A	C	47.8%	93.7 / 14.2	67	L728V (TTA→GTA)	arcB	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (24/11);  new base C (0/32);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 4.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTC  >  NC_000913/3350703‑3350989                                                                                                                                               |                                                                                                                                                  cTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTg                                                                                                                                           >  2:587869/1‑150 (MQ=255)  ttttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGAc                                                                                                                                         >  1:1062283/1‑151 (MQ=255)      gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGAcccc                                                                                                                                      <  2:32770/150‑1 (MQ=255)        cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                   <  2:75320/151‑1 (MQ=255)         acccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCccgccg                                                                                                                                  >  1:300785/1‑151 (MQ=255)                       cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                     >  1:831487/1‑150 (MQ=255)                          ttCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGt                                                                                                                 >  2:639750/1‑151 (MQ=255)                           tCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtc                                                                                                                >  1:27299/1‑151 (MQ=255)                              aCGTCGTGACGCCCCTCTTCTTTCATCTCTAAAATCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGTCTGAATTTGCTGCCCCAGCTGTT‑GCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGtcct                                                                                                              <  2:318939/150‑1 (MQ=255)                                 tcgtGCCCCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGCCGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTT‑GCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                   <  2:355645/142‑1 (MQ=255)                                   gtgccccccctctTCTTTCATCTCTCAAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTT‑GAAGATGGCGTAACCCCCCTGCCCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                            <  2:1086410/138‑1 (MQ=255)                                    tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAAccc                                                                                                                                                <  1:724880/110‑1 (MQ=255)                                    tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAAccc                                                                                                                                                >  2:724890/1‑110 (MQ=255)                                       ccccccttttcttcCATCTCTCCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAATTTTGCTGACCCAGCTGTT‑GAAGATGGCGAAACCCCCCGGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                           <  2:1179378/136‑1 (MQ=255)                                       cccccctctTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGAGGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                                    <  2:898273/146‑1 (MQ=255)                                       cccccctctTCTTTCATCTCTACAACCCATCCCCCGCCGTTATCTCCCCAGGCCGAAAGGTCAGGAGACTGATTTTGCTGCCCCATTTTTT‑GCAGATGGCTCACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaa                                                                                                       <  1:130364/143‑1 (MQ=255)                                        ccccctctTCTTTCATCTCTACATCCCATCCCCCGACTTTTTCTCCCCAGGCGGGAGGGTCAGGAGACTGATTTTGCTTCCCAATCTGTT‑GCAGATGCCGTAACCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                                    <  1:144996/146‑1 (MQ=255)                                         cccctctTCTTCCATCCCTTCATCCCATCCCCCGCCTTTATCTCCCCAGGCCGGAAGGACGGGAGACTGATTTTTCTGCCCCATCTGTT‑GAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                  <  1:142906/148‑1 (MQ=255)                                           aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                >  2:806888/1‑151 (MQ=255)                                                tcttTCCTCTCTCCAACCCATCCACGAACTTTTTCTCCCCAGGCCGGAGGGTCGGGAGACTAAATTTGCTGACCCAGTTGTT‑GCAGATGCCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                             <  2:1157949/149‑1 (MQ=255)                                                  tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                          >  2:377814/1‑150 (MQ=255)                                                  tttCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGTTCAGGAGACTAATTTTGCTCCCCCAGTTTTT‑GCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCGtttt                                                                                            <  1:975246/148‑1 (MQ=255)                                                  tttCATCTCTCCTCCCCATCCCCCGCCGTTATCTTCCCGGGCCGAAAGTTCAGGAGATTGATTTGGCTGACCCATCTTTT‑GCAGATGGCTAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                              <  1:977394/146‑1 (MQ=255)                                                   ttCATCTCTTCAATCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTCCCCCATTTTTT‑GCAGATGGCTTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                          <  1:189323/149‑1 (MQ=255)                                                          ctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc                                                                                 >  2:19820/1‑151 (MQ=255)                                                          ctTCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAGGGCCAGGAGACTGAATTTGCTGCCCCAGCTGTT‑GCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                   <  2:1058962/149‑1 (MQ=255)                                                              aacccATTCCCGGACTTTATCTCCCAGGGCGGGAAGGTCGGGAGACGAAATGTGCTCCCCCAGTTTTT‑GCAGTGGCCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                <  1:432149/145‑1 (MQ=255)                                                              aacccATCCCCCGACGTTATCTCCCCGGCCCGGAAGGTCAGGAGACTAAATTTGCTGCCCCAGTTGTT‑GAAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                               <  1:441907/146‑1 (MQ=255)                                                              aacccATCCACCGACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTTTT‑GAAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                               <  2:945295/146‑1 (MQ=255)                                                               accctttCACCGACGTTATCTACCAAGACCGGAGTGAAAGGAGCATAAATTTCCTCCCCTAGCTTTT‑GCAGGGGGCATACCCCCATTTCCTCCGCCGCCCCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                            <  1:945692/146‑1 (MQ=255)                                                                      ccccGACTTTTTTTTCCCAGGCGGGAGGGTCGGGAGACTAAATTTGCTACCCCAGTTTTT‑GCAGATGGCGTAACCCCTCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGa                                                                     <  1:496800/149‑1 (MQ=255)                                                                      ccccGACTTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGTTTTT‑GCAGAGGGCGTAACCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                        <  1:621159/146‑1 (MQ=255)                                                                      ccccGACGTTATCTCCCCAGGCCGGAGGGTCAGGAGACTAAATTTCCTCCCCCAGTTGTT‑GCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                       <  1:643235/147‑1 (MQ=255)                                                                      cACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                      >  2:911037/1‑150 (MQ=255)                                                                               tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                            >  2:769413/1‑151 (MQ=255)                                                                                 tctcCCCGGGCCGGAGGGTCAGGAGACTGAATTTGCTGACCCAGTTGTT‑GCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACg                                                           <  1:443376/150‑1 (MQ=255)                                                                                     cccAGGCGGGAAGGTAAGGAGACTAAATTTGCTGACCCAGCTGTG‑GCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                         <  2:833210/148‑1 (MQ=255)                                                                                     cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                      >  2:743395/1‑151 (MQ=255)                                                                                         ggCCGGAAGGTCAGAAGACTGATTTTCCTCCCCCAGTTTTT‑GAAGATGGCGAACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                                   <  1:21397/150‑1 (MQ=255)                                                                                            cGGAGGGTCAGGAGACTTAATTTGCTGCCCCAGCTGTT‑GCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                                   <  1:1092655/147‑1 (MQ=255)                                                                                            cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                               >  1:1060618/1‑151 (MQ=255)                                                                                             ggAGGGTCGGGAGACTAATTTTCCTGCCCCTTCTTTT‑CAAGATGGCATACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                    <  1:735368/145‑1 (MQ=255)                                                                                               aaGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                            >  2:121389/1‑151 (MQ=255)                                                                                                      ggAGACTGATTTTGCTCCCCCAGTTGTT‑GCAGAGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatca                                        <  1:728453/148‑1 (MQ=255)                                                                                                       gagaCTGAATTTCCTGACCCAGCTGTT‑GCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct                                     <  1:926959/150‑1 (MQ=255)                                                                                                        agaCTGAATTTGCTGCCCCAGCTGTT‑GCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct                                     <  2:663418/149‑1 (MQ=255)                                                                                                        agaCTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                   >  2:363579/1‑151 (MQ=255)                                                                                                           cTGAATTGGCTGCCCCAGCTGTT‑GCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                  <  2:896221/149‑1 (MQ=255)                                                                                                             taattttcccaCCCAGCTTTTGGCAGTGGCCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                  <  1:587860/138‑1 (MQ=255)                                                                                                             gcaTTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                               >  1:748747/3‑150 (MQ=255)                                                                                                             gAATTTGCTAACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                               >  2:868832/1‑150 (MQ=255)                                                                                                              attttGCTGACCCAGCTGTT‑GCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                                <  1:564747/146‑1 (MQ=255)                                                                                                              aaTTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                             >  2:394983/1‑151 (MQ=255)                                                                                                               aTTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTAACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                             >  1:667986/1‑150 (MQ=255)                                                                                                                tttGCTCCCCCAGCTGTT‑GCAGAGGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                  <  1:61991/144‑1 (MQ=255)                                                                                                                    cTGCCCCAGCTTTT‑GCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                         <  2:847461/149‑1 (MQ=255)                                                                                                                     tgccccAGTTGTT‑GCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                         <  1:732751/145‑1 (MQ=255)                                                                                                                        cccAGCTGTT‑GCAGAGGGCTTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                          <  2:661277/144‑1 (MQ=255)                                                                                                                          cAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCg                 >  1:115499/1‑151 (MQ=255)                                                                                                                              tGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg             >  2:1119959/1‑151 (MQ=255)                                                                                                                               gTT‑GCAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg              <  2:168213/149‑1 (MQ=255)                                                                                                                                tt‑GCAGAGGGCGAAACCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt            <  2:1143894/150‑1 (MQ=255)                                                                                                                                 t‑GCGGATGCCTTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg           <  1:450482/150‑1 (MQ=255)                                                                                                                                      gATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa       <  1:790093/150‑1 (MQ=255)                                                                                                                                       aTGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGc     >  1:869048/1‑151 (MQ=255)                                                                                                                                         ggCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCtt   <  1:947781/151‑1 (MQ=255)                                                                                                                                          gCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc  >  1:798110/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                  CTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTC  >  NC_000913/3350703‑3350989

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