Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	42.4%	112.1 / 14.5	66	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 (31/7);  new base C (0/28);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 2.29e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAC  >  NC_000913/3350706‑3350997                                                                                                                                                     |                                                                                                                                                 ttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcc                                                                                                                                                <  2:978685/150‑1 (MQ=255)     cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                           <  1:583245/151‑1 (MQ=255)           cGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCt                                                                                                                                     >  1:920908/1‑151 (MQ=255)             cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                   >  2:958997/1‑151 (MQ=255)                    ccctcccaCGTCGTGCCCCCCCTCTTCTTTCTTCTCTCC‑AACCCATTCCCCGCCGTTTTCTCCCCAGGCGGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                                                                             <  2:913994/143‑1 (MQ=255)                    cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                             >  2:153621/1‑150 (MQ=255)                        tCTACGTCGTGCCCCCCCTTTTCTTTACCCCCTCC‑ATCCCATCCCCTGCCGTAATCTCCCCGGCCGGGAGGGCCGGGAGACTAATTTGGCTCACCCAGCTTTTGCAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATg                                                                                                                          <  1:824386/149‑1 (MQ=255)                        tCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtc                                                                                                                        >  1:38327/1‑151 (MQ=255)                          tACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                      >  2:909558/1‑151 (MQ=255)                           ccgtcgtGCCCCCCCTCTTTTTTCATCTCTCC‑AACCCATCCACCGACTTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                     <  2:508502/150‑1 (MQ=255)                            cgtcgtGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                                    >  2:207068/1‑151 (MQ=255)                               cgtGCCGCCACTCTTCTTTCATCTCTTC‑AATCCATCCCCCGACGTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTGCCCCATCTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                  <  2:116211/150‑1 (MQ=255)                               cgtGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                 >  1:188319/1‑151 (MQ=255)                                gtGACGCCCCTCTTCTTTCATCTCTTA‑AATCCATCCCCCGACGTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                  <  2:46123/149‑1 (MQ=255)                                gtGACCCAACTCTTCTTAC‑TCTTTTCAAATCCTTCCCCCGATGTTATTTTCCGAGGCGGGAGGGACAGGAATTTGAATTTGTAACCCCGGCTTTTGTAGATGGCGAACCCCCACGACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                                <  1:868043/151‑1 (MQ=255)                                   aCGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                              >  2:659869/1‑150 (MQ=255)                                    cccccctctTTTTCCATCTCTTA‑AATCCATCCCCCGCCTTTATCTCCCCAGGCCGGAGGGTCAGGAGACTAATTTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                        <  1:28848/134‑1 (MQ=255)                                    cccccctctTCTTTCATCTCTCC‑ACCCCATCCCCCGCCTTTATCTCCCCGGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                                            <  1:1072306/146‑1 (MQ=255)                                       ccctctTCTTCCCTCTCTCC‑AACCCATTCCCCGCCGTTATCTCCCCAGCCGGGAAGGTCAGGAGACTGATTTTGCTGCCCCTTCTTTTGCAGATGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                         <  1:5083/149‑1 (MQ=255)                                            tttttccttccctac‑AATCCATCCCCCGCCTTTCTCTCCCCAGGCGGGAAGGGCAGGAGACTGATTTGGCTGCCCCCGCGGTTGCAGATGTGGTACCCCCACTGACGCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                     <  2:692534/136‑1 (MQ=255)                                               tttCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                  >  1:220165/1‑150 (MQ=255)                                               tttCATCTCTCC‑AACCCATCCACCGACTTTATCTCCCCAGGCGGGAGGGTCAGGAGACTGATTTGGCTGCCCCAGTTGTTGAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                     <  1:1066816/147‑1 (MQ=255)                                                ttCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                                >  1:652065/1‑151 (MQ=255)                                                    tctcTTC‑AACCCATCCCCCGACTTTATCTCCCCGGGCGGGAGGGTCGGGAGACTGATTTTGCTGCCCAAGTTGTTGCAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                               <  2:674261/148‑1 (MQ=255)                                                     ctctTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                           >  1:1112389/1‑151 (MQ=255)                                                       ctTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc                                                                                         >  1:127525/1‑151 (MQ=255)                                                         tC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                       >  2:945750/1‑151 (MQ=255)                                                          c‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGt                                                                                      >  2:838136/1‑151 (MQ=255)                                                               ccATTCCCCGACGTTATCTCCCCAGGCCGGAAGGCCAGGAGACTGATTTTGCTGACCCAGTTGTTGCAGGTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGa                                                                                  <  1:528580/151‑1 (MQ=255)                                                                cATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                     >  2:872790/1‑147 (MQ=255)                                                                 aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                                >  2:306624/1‑151 (MQ=255)                                                                          cGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                         >  1:817240/1‑149 (MQ=255)                                                                           gTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                      >  1:548104/1‑151 (MQ=255)                                                                            ttctCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTCATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                      >  1:781222/4‑150 (MQ=255)                                                                             tATCTCCCCAGGCCGAAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCa                                                                     <  1:899024/150‑1 (MQ=255)                                                                              aTCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACg                                                                   >  2:213806/1‑151 (MQ=255)                                                                                   cccGGGCCAGAAGGTCGGGAAACTGATTTTGCGCCCCCGGGTGTTGCAGAGGGCGTAACCCCTCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                    <  2:1549/145‑1 (MQ=255)                                                                                     cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                            >  2:908996/1‑151 (MQ=255)                                                                                       ggCCGAAAGGACAGGAAAGTGAATTTGCTGCCCGACCTGTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAg                                                          <  2:620488/151‑1 (MQ=255)                                                                                          cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                       >  1:177881/1‑151 (MQ=255)                                                                                            gAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                     >  1:391795/1‑151 (MQ=255)                                                                                                  cgggAGACTGAATTTGCAGACCCAGCTGTTGCAGAGGGCGTAACCCCCATGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                    <  2:923054/144‑1 (MQ=255)                                                                                                  cAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                               >  1:604190/1‑151 (MQ=255)                                                                                                    ggAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                               <  2:261694/149‑1 (MQ=255)                                                                                                    ggAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGAAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                              <  1:112856/150‑1 (MQ=255)                                                                                                    ggAGACTAAATTTGCTGACCCAGCTGTTGCAGATGGCGAACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                              <  2:1170744/150‑1 (MQ=255)                                                                                                     gagaCTGAATTTGCTGACCCAGTTGTGGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                            <  2:1107922/151‑1 (MQ=255)                                                                                                         ctaattttGCTGCCCCAGTTTTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                          <  2:823734/144‑1 (MQ=255)                                                                                                         atGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                         >  1:836295/2‑150 (MQ=255)                                                                                                           gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                                       >  1:648667/1‑150 (MQ=255)                                                                                                              tttGCTGACCCAGCTTTTGCAGAGGGTTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca                                    <  2:937933/150‑1 (MQ=255)                                                                                                              ttggCTCACCCACTTTTTGCAGATGACGTACCCCCTCTACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                                     <  2:224216/146‑1 (MQ=255)                                                                                                                      cccgttttttGCAGAGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc                           <  1:933752/143‑1 (MQ=255)                                                                                                                         aGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                        >  2:244016/1‑151 (MQ=255)                                                                                                                            tGTTGCAGATGCGGAACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                      <  2:2101/150‑1 (MQ=255)                                                                                                                                    aTGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGc             >  1:376663/1‑151 (MQ=255)                                                                                                                                     tGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt            >  1:888651/1‑151 (MQ=255)                                                                                                                                      ggCGCACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt            <  2:778459/150‑1 (MQ=255)                                                                                                                                      ggCGAAACCCCCCGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc                <  2:325980/146‑1 (MQ=255)                                                                                                                                        cGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc          <  1:700469/150‑1 (MQ=255)                                                                                                                                           acccccACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg         <  2:403369/146‑1 (MQ=255)                                                                                                                                           acccccACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg         <  1:331965/146‑1 (MQ=255)                                                                                                                                           aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc      >  1:625296/1‑151 (MQ=255)                                                                                                                                             ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt    <  2:451812/151‑1 (MQ=255)                                                                                                                                               ccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc  <  1:125480/148‑1 (MQ=255)                                                                                                                                                cACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc  >  1:877364/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTC‑AATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAC  >  NC_000913/3350706‑3350997

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