Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,844	0	A	C	43.8%	107.5 / 12.8	64	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 (32/4);  new base C (0/28);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTC  >  NC_000913/3350695‑3350993                                                                                                                                                      |                                                                                                                                                      tttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTaa                                                                                                                                                       <  2:1115236/150‑1 (MQ=255)  ttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAAc                                                                                                                                                      <  2:426319/150‑1 (MQ=255)     gTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAAcccc                                                                                                                                                   >  1:342676/1‑150 (MQ=255)             tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                          >  1:692621/1‑151 (MQ=255)              gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                          <  1:557843/150‑1 (MQ=255)               ccacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgc                                                                                                                                        <  1:554418/151‑1 (MQ=255)                  cccacGCTTTCAGCACTTCTACGTCGTGCCCCCCCTCTTCTTTCATCTCTTCAATCCATTCCCCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCcgccgc                                                                                                                                     <  2:894598/151‑1 (MQ=255)                        cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                               >  1:1188160/1‑151 (MQ=255)                            cAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtt                                                                                                                           >  2:938416/1‑151 (MQ=255)                               cccTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCCCCGACGTTATCTTCCCAGCCCGGAAGGTCAGGAGACTAAATTTGCTGACCCAGTTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                         <  1:927807/148‑1 (MQ=255)                                  ttCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGt                                                                                                                     >  1:10444/1‑151 (MQ=255)                                     tACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                  >  2:1176446/1‑151 (MQ=255)                                     tACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                  >  2:728330/1‑151 (MQ=255)                                      acgcggtGCCGCCCCTCTTCTTCCATCTCTTAAACCCATCCACCGACTTTTTCTCCCCAGGCGGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGTTTTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                  <  1:400449/145‑1 (MQ=255)                                       cgtcgtGACCCCCCTCTTCTTTCATCTCTTCAACCCATCCACCGACTTTATCTCCCCAGGCCGGAGGGTCAGGAGATTAATTTTGCTCCCCCGTTTGTTGCAGTGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                 <  2:943371/150‑1 (MQ=255)                                       cgtcggGACCCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGACTTTATCTCCCCGGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGAGGCCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcc                                                                                                                   <  2:597587/148‑1 (MQ=255)                                       cgtcggCACGCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGACTTTTTCTTCCCGGGCGGGAGGGTCAGGAGACTAAATTTGCGGCCCCAGCTGTTGCAGGTGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                  <  1:196111/149‑1 (MQ=255)                                           gggccccccctctTCTTTCATCTCTTAAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCGGGAGACTAATTTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                             <  1:772860/141‑1 (MQ=255)                                            tgccccccctctTCTTTCATCTCTCCAACCAATCCCCGGACTTTTTCTCCCCGGGCGGAAAGGTCAGGAAAGTAATTTTGCTCACCCTTTTGTTGAAGTGGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                                <  2:396713/138‑1 (MQ=255)                                               cccctctctTCTTCCATCTCTGCAACCCATTCCCCGCCTTTATCTCCCCAGGCGGGAAGTCCAGGAGTCTGATTTTGCTGCCCCGGCTGTGGAAGATGGTATACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaa                                                                                                           <  1:612185/143‑1 (MQ=255)                                               cccccctctTCTTTCTTCTCTTCAACCCATCCCCCGCCTTTTTCTCCCCAGGCGGGAAGGTCAGGAGATTGATTTTGCTGCCCCGGTTGTTGAAGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                                        <  2:1063646/146‑1 (MQ=255)                                               cccccctctTCTTTCATCTCTCCAACCCATCCCCCGCCGTTTTCTCCCCGGGCCGGAAGGTCAGAAGACTGAATTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                               <  2:461836/139‑1 (MQ=255)                                                 ccACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                      >  1:884119/1‑151 (MQ=255)                                                 ccACTCTTCTTTCATCTCTTCAACCCATCCCCCGACTTTATCTCCCCAGGCGGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCTTTAATCTTATGTCCTTCCTCAACAAt                                                                                                       <  1:604221/150‑1 (MQ=255)                                                  ccctctTCTTTCATCTCTTCAATCCATTCCCCGACGTTATCTCCCCGGGCCGGAAGGTCAGGAGACTAAATTTGCTGCCCCAGTTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                     <  1:914520/149‑1 (MQ=255)                                                         cttTCATCTCTTCAACCCATTCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                  <  1:1086616/147‑1 (MQ=255)                                                          tttCATCTCTTCAACCCATTCCCCGACGTTATCTCCCCAGGCCGGAAGGTCGGGAGACTGATTTTGCTGCCCCAGTTGTTGAAGATGGCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                  <  2:1170628/146‑1 (MQ=255)                                                          tttCATCTCTTCAACCCATCCACCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                                <  2:894905/148‑1 (MQ=255)                                                           ttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                            >  2:911314/1‑151 (MQ=255)                                                           ttCATCTCTTCAACCCATCCCCCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGCTTTTGCAGATGGCTTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTa                                                                                             <  1:550593/150‑1 (MQ=255)                                                                ctctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                       >  2:699365/1‑151 (MQ=255)                                                                ctctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                       >  2:824434/1‑151 (MQ=255)                                                                 tctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                      >  1:956176/1‑151 (MQ=255)                                                                   ttCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                    >  2:1153753/1‑151 (MQ=255)                                                                        tCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                 >  2:1222547/1‑149 (MQ=255)                                                                         ccATCCCCCGACGTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTGACCCAGTTTTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                <  2:1001429/149‑1 (MQ=255)                                                                           aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                            >  2:1017265/1‑151 (MQ=255)                                                                           aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCCCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                            >  1:23185/1‑151 (MQ=255)                                                                             cccccgcctttATCTCCCCAGCCCGAAAGGTCAGGAGACTGATTTTGCTCCCCCAGTTTTGGAAGATGGCGAACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                           <  2:946761/141‑1 (MQ=255)                                                                              ccccgcctttATCTCCCCAGGCCGGAGGGTCAGGCATATAATTTGGCTGCCCCATTTTTTGCAGTGGCCGTACCCCCATGCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                             <  2:400049/139‑1 (MQ=255)                                                                              cACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                          >  2:237753/1‑150 (MQ=255)                                                                                  gACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                     >  1:662413/1‑151 (MQ=255)                                                                                   aCGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                     >  2:882380/1‑150 (MQ=255)                                                                                      ttctctcCCCGGGCGGGAGGGCCGGGAGACTAATTTTGCTGCCCCAGTTTTTGCAGATGCCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                    <  2:103796/145‑1 (MQ=255)                                                                                       tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                >  2:680389/1‑151 (MQ=255)                                                                                       tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                >  2:828148/1‑151 (MQ=255)                                                                                                 ggCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                        <  1:1036971/149‑1 (MQ=255)                                                                                                  gCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGc                                                     >  2:896213/1‑151 (MQ=255)                                                                                                    cGGGAGGGCAGGAGACTAAATTTGCTCACCCGGTTGTTGGAGGTGGCCTACCCCCCCGGACCCCGCCGCACCTTTAAATTTAAGTCCTTCCCCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGcc                                                    <  2:903626/150‑1 (MQ=255)                                                                                                    cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                   >  2:921968/1‑151 (MQ=255)                                                                                                      gAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                 >  1:73954/1‑151 (MQ=255)                                                                                                         ggtcggggtaCTGATTTTGCTCCCCCAGTTGTTGCAGATGGCATACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGca                                               <  2:10444/141‑1 (MQ=255)                                                                                                         ggTCAGGAAATTGAATTGGCTGCCCCGGCTGTGGCAAAGGCCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                 <  2:37485/148‑1 (MQ=255)                                                                                                            cAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCC‑TTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                          >  2:327001/1‑151 (MQ=255)                                                                                                               gagaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                        >  2:688825/1‑151 (MQ=255)                                                                                                                     gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                                   >  2:95488/1‑150 (MQ=255)                                                                                                                        tttGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAacac                               >  1:851360/1‑151 (MQ=255)                                                                                                                          tccgctcacAGCTTTTGCAGAGGCCTTAACCCCACGTCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                              <  2:766920/142‑1 (MQ=255)                                                                                                                                cccAGCTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTaa                         <  1:570378/149‑1 (MQ=255)                                                                                                                                     cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                  >  1:951891/1‑151 (MQ=255)                                                                                                                                       gTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt                >  1:529968/1‑151 (MQ=255)                                                                                                                                        ttGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg               >  1:418078/1‑151 (MQ=255)                                                                                                                                         tGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg               >  1:984665/1‑150 (MQ=255)                                                                                                                                            aattGGCGAACCCCCCCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                 <  2:320582/142‑1 (MQ=255)                                                                                                                                                ggCGTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCtt       <  1:972184/151‑1 (MQ=255)                                                                                                                                                     aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc  >  2:1220603/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTC  >  NC_000913/3350695‑3350993

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