Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	56.7%	70.0 / 97.9	67	S725A (TCA→GCA)	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 (23/6);  new base C (0/38);  total (23/44)
Fisher's exact test for biased strand distribution p-value = 8.95e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-10

TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350704‑3351003                                                                                                                                                       |                                                                                                                                                       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)                                      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)                                                                                                                                              cccccACTGCCCCCGCCGCACCTTCAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc           <  2:570142/150‑1 (MQ=255)                                                                                                                                               ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGt             <  2:244134/148‑1 (MQ=255)                                                                                                                                                   aCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACa       >  2:709368/1‑150 (MQ=255)                                                                                                                                                       aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  1:93137/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTT‑GCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350704‑3351003

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