Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	40.9%	114.6 / 11.0	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 (33/6);  new base C (0/27);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 9.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCC  >  NC_000913/3350705‑3350994                                                                                                                                                        |                                                                                                                                              tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAc                                                                                                                                              <  2:943456/150‑1 (MQ=255)   tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                           >  2:1015461/1‑151 (MQ=255)    gccgccCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                           <  2:573942/150‑1 (MQ=255)       acccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgccg                                                                                                                                       >  2:576555/1‑151 (MQ=255)           acGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAcc                                                                                                                                   >  1:602808/1‑151 (MQ=255)              cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                >  1:381800/1‑151 (MQ=255)               tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGCCTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACACCACTGACCCCGCCGCACCt                                                                                                                                  >  1:218299/1‑148 (MQ=255)                ttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTa                                                                                                                               >  1:1047429/1‑150 (MQ=255)                     cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                          >  2:235984/1‑150 (MQ=255)                          cTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcc                                                                                                                    >  2:570384/1‑151 (MQ=255)                             cgtggtGCCCCCCCTCTTCTTTCCTCTCTACAA‑CCCATCCCCCGACGTTATCTCCCCGGGCCGGAA‑GGTCAGGAGACTAAT‑TTTCCTGACCCAGTTGTTGAAGAGGGCGTACCCCCACAGCCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                     <  2:924517/147‑1 (MQ=255)                              ttcttgcccccACTCTTTTTTCACCCCTCAAA‑TCCTTCCACCCAGTTTATCTTCCCGGGCCGGAA‑GGCCAGGAGATTTAT‑TTGTCTGACCCTGCTGTTGATGTTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                        <  2:96083/134‑1 (MQ=255)                              gtggtgccccccctctTCTTTCATCT‑TTCAACTCCATCCCCCGAATTTTTCTCCCCGGCCGGGAG‑GGTCGGGAGACTAAA‑TTTGCGAACCCGGTTGTTGAAGGTGGCGTACCCCCCCGGACCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                   <  1:842845/136‑1 (MQ=255)                                 gtGCCGCCCCTTTTTTTTCTTCTCTTCAA‑CCCATCCCCGGACTTTTTCTCCCCAGGCGGGAG‑GCTCGGGAAACTGAT‑TTTGTTCCCCCGTTTTTTGAAGTGGGCGTACCCCCACGCCCCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                                <  2:216569/148‑1 (MQ=255)                                 gtGACTCCTCTCTTCTTTCCCCCCTCAAA‑CCCATCCCCCGGCTTTATCTCCCCAGGCCGAAA‑GTGCGGGCGACTGAT‑TTGGCTGCCCCGGTTTTTGCAATTGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                  <  2:599808/146‑1 (MQ=255)                                     cGCCCCTCTTCTTCCATCTCTTCAA‑CCCATCCCCCGACGTTATCTCCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGCCCCAGTTTTTGCAGATGGCCTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                          <  2:467058/150‑1 (MQ=255)                                       cccctctTCTTTCATCTCTTCAC‑TCCATCCCCCGACGTTATCTCCCCAGGCCGGAG‑GGTCAGGAGACTGAA‑TTTGCTGCCCCGGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                                        <  1:988204/147‑1 (MQ=255)                                       ccACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                       >  2:168654/1‑151 (MQ=255)                                         aCTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                     >  2:625331/1‑151 (MQ=255)                                             ttttttCATCTCTCCAA‑CCCATCCCCCGACTTTTTCTCCCCGGGCGGGAA‑GGTCAGGAGCCTGAA‑TTTCCTGCCCCATTTGTTGCAGATGGCCTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                     <  2:3945/144‑1 (MQ=255)                                             tttttccATCTCTCCAA‑CCCATCCACCGACGTTATCTCCCCAGCCGGGAG‑GGTCAGGAGACTAAT‑TTTGCTGCCCAAGTTGTTGCGGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                       <  2:725599/139‑1 (MQ=255)                                             ttcttTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                                 >  2:515175/1‑151 (MQ=255)                                              tcttTCATCTCTCCAA‑CCCATCCCCCGACTTTATCTTCCCGGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGCCCCAGCTGTTGCAGATGGCTTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                                 <  1:1164901/150‑1 (MQ=255)                                                tttCATCTCTTCAA‑CCCATCCCCCGTCTTTATCTCCCCGGGCGGGAA‑GGTCGGGATCCTAAT‑TTTGCTGCCCCAGTTTTTGCTGAGGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                           <  2:812783/138‑1 (MQ=255)                                                 ttCATCGCTTCAA‑ACCATTCACATACGGTAGCTTCACAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGCTGGCGTAACGCCACTGACCCCGCCGCACCTTTAATTTAATGTCCTTCCTTAACAATGCCTTTTTTa                                                                                              >  2:643287/1‑150 (MQ=255)                                                        ctccAA‑TCCTTCCACCGACTTATTCTCCCCAGGCGGGAG‑GGTCGGGAGA‑TGAATTTTGCTGCCCCGGTTGTTGCAGGGGGCATAACCCCACTCACTCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                         <  2:926492/145‑1 (MQ=255)                                                         ttCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                     >  2:656513/1‑151 (MQ=255)                                                         ttCAA‑CCCATTCCCCGGCTTTATCTCCCCAGGCGGAAA‑GGTCAGGAAACTAAA‑TTTGCTGCCCAGGCTGTTGCAGGTGGCGTACCCCCCCAGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt                                                                                          <  1:424105/146‑1 (MQ=255)                                                                   tccccgatcGTAATCTCCCCAGGCCGGAA‑GTACAGGAGCCTGAT‑TTTTCCCCCCCCGCGGTTGCAGAGGGAGTACCCCCGCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                              <  1:533117/141‑1 (MQ=255)                                                                     cccctactttctctcCCCGGGCGGGAG‑GGTCAGGAGCTTAAA‑TTGGCTCACCCGGTTTTGCCAGATGACATACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                           <  2:984585/139‑1 (MQ=255)                                                                        cGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                       >  2:723139/1‑151 (MQ=255)                                                                        cGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCATCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                       >  1:552353/1‑151 (MQ=255)                                                                         gACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                      >  1:197474/1‑151 (MQ=255)                                                                           cGTTATCTCCCCGGCCCGAAA‑GGTCGGGGAACTGAT‑TTTGTTGCCCCATTTGTTGAAAATGGCGAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                       <  1:1076467/148‑1 (MQ=255)                                                                            gTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                   >  2:1195156/1‑151 (MQ=255)                                                                                 cTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                              >  2:26099/1‑151 (MQ=255)                                                                                     cccgggcGGAAGGGTCGGGAGACTAAT‑TTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                 <  2:283257/139‑1 (MQ=255)                                                                                     ccAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                           >  1:62916/1‑150 (MQ=255)                                                                                          cgggAA‑GCTCAGGAGATTAAT‑TTTCCTCCCCCAGTTGTTGCAGAGGCCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                         <  2:162951/145‑1 (MQ=255)                                                                                              aa‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                 >  2:195657/1‑151 (MQ=255)                                                                                                         gaCTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                       >  2:256791/1‑151 (MQ=255)                                                                                                            tGAT‑TTTGCGGCCCCAGTTTTTGCAGATGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                      <  2:293562/149‑1 (MQ=255)                                                                                                                     cTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                             >  2:433221/1‑150 (MQ=255)                                                                                                                       gACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTa                           >  1:970092/1‑150 (MQ=255)                                                                                                                         cccggtttttGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc                        <  1:1035726/143‑1 (MQ=255)                                                                                                                         cccAGTTTTTGAAGATCCCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                            <  2:62916/147‑1 (MQ=255)                                                                                                                         cccAGCTGTGGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc                        <  1:885944/151‑1 (MQ=255)                                                                                                                          ccAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAccc                       >  1:103610/1‑151 (MQ=255)                                                                                                                            aGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                     >  1:287093/1‑151 (MQ=255)                                                                                                                             gCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc                    >  1:46874/1‑151 (MQ=255)                                                                                                                              cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                   >  2:322284/1‑151 (MQ=255)                                                                                                                              cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                   >  2:614301/1‑151 (MQ=255)                                                                                                                               tGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                  >  2:171435/1‑151 (MQ=255)                                                                                                                                gTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                  >  1:814705/1‑150 (MQ=255)                                                                                                                                 ttGCATATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTccg                  <  1:1032672/149‑3 (MQ=255)                                                                                                                                 ttGCAGATGGCTTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                  <  1:1080561/149‑1 (MQ=255)                                                                                                                                      gagggCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                <  1:574949/143‑1 (MQ=255)                                                                                                                                      gagggCGAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc                    <  1:415968/139‑1 (MQ=255)                                                                                                                                      gATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAg           <  1:1136384/151‑1 (MQ=255)                                                                                                                                      gATGGCGAACCCCCCCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt              <  1:1148639/148‑1 (MQ=255)                                                                                                                                      gATGCCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa            <  2:713984/150‑1 (MQ=255)                                                                                                                                           catccccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc       <  2:208639/145‑1 (MQ=255)                                                                                                                                           cGTACCCCCCCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg      <  1:48914/151‑1 (MQ=255)                                                                                                                                              aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc   >  1:503588/1‑151 (MQ=255)                                                                                                                                              aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc   >  2:602027/1‑151 (MQ=255)                                                                                                                                               aCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc  >  2:209607/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                              TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAA‑TCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAA‑TTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCC  >  NC_000913/3350705‑3350994

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