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

TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350707‑3351003                                                                                                                                                    |                                                                                                                                                       tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                                    >  2:79096/1‑151 (MQ=255)    cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                 >  2:446448/1‑151 (MQ=255)    cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                                 <  2:187262/151‑1 (MQ=255)     acccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGCCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgccg                                                                                                                                                >  2:19399/1‑151 (MQ=255)        cacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAc                                                                                                                                             <  1:543201/151‑1 (MQ=255)            cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                         >  2:790178/1‑151 (MQ=255)             tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                         >  2:625892/1‑150 (MQ=255)             ttccAGCCCTTCTCCGTCGGGACCCCACTCTTCTTTCATCTCTCCAATCCATGCCCCGACTTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTCCTCCCCCAGCTGTTGCAGAGGGCTAAACCCCACTGCCCCCGCCGCACCttt                                                                                                                                         <  2:572982/147‑1 (MQ=255)                  ccACTTCACCGTCGTGACGCCCCTCTTCTTTCTTCCCTTAAATCCCTGCCCCGTCGTTATCTCCCCAGGCC‑GAAAGGTCGGGAGACTGATTTTGCTGATCCAGTTGTTGAGGTGGGCGAACCCCCCCTGACCCCGCCGCACCTTTAAtt                                                                                                                                     <  2:643732/148‑1 (MQ=255)                   cACTTCTCCGTCGTGACCCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAAttt                                                                                                                                    <  1:338676/149‑1 (MQ=255)                      ttCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAGTTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGt                                                                                                                               >  1:827755/1‑151 (MQ=255)                      ttCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGt                                                                                                                               >  1:331119/1‑151 (MQ=255)                         tACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                            >  2:60745/1‑151 (MQ=255)                             tggtGACGCCACTCTTCTTTCCTCTCTTAACCCCATCCCCCGCCTTTATCTCCCCAGGCG‑GGAGGGCCGGGAGACTAAATTTGCTCCCCCGTTTTTTGCAGATGGCTTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcct                                                                                                                            <  2:25361/145‑1 (MQ=255)                               gtGACCCCCTTCTTCTTTCATCTCTTCAACCCATCCCCCGTCGTTATCTCCCCGGGCC‑GGAGGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                                      <  1:393642/151‑1 (MQ=255)                                   ccccctttttttttCCTCTCTTAAACCAATCCCCCGAGTTTATCTCCCCAGGCG‑GGAGGGTCAGGAGACTAAATTTGCTCCCCCAGCTGTTAAAGGGGGCGAACCCCCCCGCCCCCCGCCACACCTTTAATTTTATGTCCTTCCTca                                                                                                                      <  1:117013/136‑1 (MQ=255)                                   ccccACTCTTCTTTCATCTCTTCAATCCATCCACCGCCTTTATCTTCCCGGGCC‑GGAAGGACAGGAGACTGATTTTGCTGACCCAGTTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                        <  1:779186/143‑1 (MQ=255)                                   cGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                                   >  2:404823/1‑150 (MQ=255)                                     cccctctTTTTCCATCTCTACAACCCATTCCCCGCCGTTTTCTCCCCAGCCC‑GAAAGGTCAGGAGACTGAATTTGCTCACCCAGCTGTTGCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                                 <  1:247507/131‑1 (MQ=255)                                     ccACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                                >  2:529703/1‑151 (MQ=255)                                      cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                               >  2:349250/1‑151 (MQ=255)                                        ctctTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                                             >  2:243603/1‑151 (MQ=255)                                        ctctTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                                             >  2:761952/1‑151 (MQ=255)                                        ctctTCTTTCATCTCTTAAACCCATTCCCCGGCTTTATCTCCCAGGGCC‑GGAAGGTCAGGAGTCTTATTTTGCTGCCCCATTTTTTTCAAAGGCCTCACCCCCACAGCCGCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                                                <  2:390136/148‑1 (MQ=255)                                          cttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                           >  1:5905/1‑151 (MQ=255)                                              tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                        >  1:482856/1‑150 (MQ=255)                                               ttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                                      >  1:810050/1‑151 (MQ=255)                                                  aTCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt                                                                                                   >  1:554785/1‑151 (MQ=255)                                                     tctTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                                >  1:545336/1‑151 (MQ=255)                                                       tacAATCCATCCCCCACCGTTTTCTCCCCAGCCT‑GAAAGGACAGGAACCTAATTTTGCTCACCCAGTTTTTGCAGAGGCCGAACCCCCCCGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                         <  2:624048/138‑1 (MQ=255)                                                          aacccATCCCCCGACTTTATCTCCCCAGGCC‑GGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                                <  1:377627/143‑1 (MQ=255)                                                          aaTCCTTCCCCCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTAATTTTGCGGACCCACCTGTGGCAGATGGCGCAACCCCCCCCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                                  <  1:611267/144‑1 (MQ=255)                                                           acccATCCCCCGCCTTTATCTCCCCAGGCG‑GGAAGGTAAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGTTGGCTTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                             <  2:170186/146‑1 (MQ=255)                                                           aTCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                           >  1:436243/1‑150 (MQ=255)                                                             caaTTCCCCTACTTTACCTCCCCAGGCC‑GGAAGGTCAGGAGACTAATTTGGCTGCCCTAGCTGTTGCAGATGGCTTACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                          <  1:163903/147‑1 (MQ=255)                                                               aTTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                                      >  1:35851/1‑151 (MQ=255)                                                                  ccccggttttttctcCCCGGGCC‑GGAAGGGCAGGAGACTAATTTTGCTGCCCCGGTTGTTGGAGTGGGAGTACCCCCCCTTCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                          <  2:31100/133‑1 (MQ=255)                                                                  ccccgcctttATCTCCCCAGGCG‑GGAGGGTCAGGAGATGGATTTGGCTGCCCCATTTGTTGCAGAGGGCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                       <  2:381043/139‑1 (MQ=255)                                                                    ccaccttaccctccccAGGCG‑GAAAGGTCGGGAGCCGAAATTCGCTCCCCCATTTTTTGCAAAGGCCTTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                           <  2:484624/128‑1 (MQ=255)                                                                     cGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                                >  1:647320/1‑151 (MQ=255)                                                                          ttATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                            >  2:722485/1‑150 (MQ=255)                                                                          ttATCTCCCCAGGCC‑GGAGGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCAAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                               <  1:241552/147‑1 (MQ=255)                                                                            tttctCCCCGGCCGGGAGGGACAGGAGACTAAATTTCCTGCCCCGGTTTTTGCAGATGGCTAACCCCCATTCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                            <  2:741736/145‑1 (MQ=255)                                                                                 cccGGGCG‑GGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                                    <  1:768144/151‑1 (MQ=255)                                                                                 cccAGGCC‑GGAGGGTCAGGAGACTGAATTTGCTCCCCCAGTTGTTGCAGATGCCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGc                                                                        <  2:449503/147‑1 (MQ=255)                                                                                  ccAGGCG‑GAAAGGTCAGGACACAAATTTGGCTCCCCCCGTTGTTGAAGTGGTCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCa                                                                           <  2:23670/143‑1 (MQ=255)                                                                                   cAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTAACAt                                                                  >  1:332074/1‑151 (MQ=255)                                                                                        c‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTCTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                             >  2:219686/1‑151 (MQ=255)                                                                                          ggAGGGTCGGGACACTGATTTTGCTGCCCCTGTTTTTGCAGAGGGCTCACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                             <  2:232208/150‑1 (MQ=255)                                                                                          ggAGGGCCAGGAGACTGATTTTGCTCCCCCAGTTTTTGAAGATGGCGACCCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                             <  2:321776/150‑1 (MQ=255)                                                                                            aaGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                          >  1:301468/1‑151 (MQ=255)                                                                                              ggTCAGAGGACTGAATTTGCTCACCCAGCTGTTGCTGATGCCCTCCCCCCACTCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                           <  2:719687/148‑1 (MQ=255)                                                                                                    gagaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                                  >  1:367297/1‑151 (MQ=255)                                                                                                              ttGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAacac                                         <  1:479691/150‑1 (MQ=255)                                                                                                                  cgccccAGCCGCTGCGGATGAAGTAACCCTCCCATCCCCGCTGCCCCTTTAAGTTTATGTCCCTCCTCAACAATGCCTTTTTTATCCTGCGCCTTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcct                                       <  2:718821/145‑2 (MQ=255)                                                                                                                    aCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                                  >  1:219045/1‑151 (MQ=255)                                                                                                                     cccgccTGTTGAAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc                                 <  1:592717/146‑1 (MQ=255)                                                                                                                     cccATCTGTTGCAGATGGCTCCCCCCTCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACACCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                                  <  2:871571/150‑1 (MQ=255)                                                                                                                     ccatgCTGTTGAAGAGGGCTTCCCCCCCCAACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAacac                                         <  2:20014/139‑1 (MQ=255)                                                                                                                     ccaggCTTTTGCAGAGGCCTTACCCCCACGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTaa                                   <  2:827568/145‑1 (MQ=255)                                                                                                                          gtgttgggggtGGCATACCCCCCCGGCCTCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc                                 <  1:915114/136‑1 (MQ=255)                                                                                                                           tGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCTACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAACCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                           >  1:509143/1‑151 (MQ=255)                                                                                                                           tGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                           >  1:668546/1‑151 (MQ=255)                                                                                                                              tGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                         >  2:533004/1‑150 (MQ=255)                                                                                                                                 aGATGGCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                     <  1:558028/151‑1 (MQ=255)                                                                                                                                  gttgcctaacccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc                      <  2:331031/139‑1 (MQ=255)                                                                                                                                  gATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAg                    <  1:369564/151‑1 (MQ=255)                                                                                                                                       cGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg               >  2:26672/1‑151 (MQ=255)                                                                                                                                           cccccGCTCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc           <  1:746765/150‑1 (MQ=255)                                                                                                                                           aCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc           >  1:93388/1‑151 (MQ=255)                                                                                                                                                aCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCC‑TTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACa       >  1:159869/1‑149 (MQ=255)                                                                                                                                                    aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATCGCCCGGCATCATCTTCTCACACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  2:295296/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCC‑GGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350707‑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: ‑