Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,840,799	G→C	11.1%	S202W (TCG→TGG) ‡	ynjF ←	putative phosphatidyl transferase, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,840,799	0	G	C	11.1%	220.4 / 15.2	72	S202W (TCG→TGG) ‡	ynjF	putative phosphatidyl transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/36);  new base C (2/6);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTT  >  NC_000913/1840652‑1840942                                                                                                                                                     |                                                                                                                                                aTGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                                                                                                                                               >  2:1168674/1‑150 (MQ=255)  tGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACt                                                                                                                                              <  1:1038529/150‑1 (MQ=255)  tGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACt                                                                                                                                              >  1:1210100/1‑150 (MQ=255)    tACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                                                                                                                                            <  2:954917/150‑1 (MQ=255)    tACGCACGCGCCATAGGTTCGAAGAATCTTTCAGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGCCGCTGGCGCGACTTc                                                                                                                                            >  2:528602/1‑150 (MQ=255)     acgcaggcGGCATGGGTTGTAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAGTTCa                                                                                                                                           <  2:1520121/150‑1 (MQ=255)        cacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCg                                                                                                                                        <  2:52318/150‑1 (MQ=255)        cacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCg                                                                                                                                        <  1:705833/150‑1 (MQ=255)        cacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCg                                                                                                                                        <  1:955121/150‑1 (MQ=255)          cgcgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                      >  1:394417‑M2/8‑150 (MQ=255)          cgcgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTc                                                                                                                                      <  2:331227‑M2/143‑1 (MQ=255)          cgcgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTc                                                                                                                                      <  1:32306‑M2/143‑1 (MQ=255)            cgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAg                                                                                                                                    <  1:977339‑M2/145‑1 (MQ=255)            cgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAg                                                                                                                                    <  1:1352882‑M2/145‑1 (MQ=255)            cgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAg                                                                                                                                    <  1:1561962‑M2/145‑1 (MQ=255)            cgccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCAGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAg                                                                                                                                    >  1:1046970‑M2/6‑150 (MQ=255)             gccaTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCagt                                                                                                                                   <  1:500966‑M2/146‑2 (MQ=255)                 tGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc                                                                                                                               <  1:188214/150‑1 (MQ=255)                       ggAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTcc                                                                                                                         <  1:4829/150‑1 (MQ=255)                       ggAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCAGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAATTCAGCGTCAGATACCCACTcc                                                                                                                         <  2:1046970/150‑1 (MQ=255)                         aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcac                                                                                                                       >  2:700541/1‑150 (MQ=255)                           gaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCacac                                                                                                                     >  1:1030925/1‑150 (MQ=255)                              tGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATCGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcg                                                                                                                  <  2:809144/150‑1 (MQ=255)                                tttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg                                                                                                                >  1:467492/1‑150 (MQ=255)                                  tCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTa                                                                                                              <  2:886520/150‑1 (MQ=255)                                  tCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTa                                                                                                              <  1:835597/150‑1 (MQ=255)                                       gtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAgg                                                                                                          <  1:320895/150‑1 (MQ=255)                                        tATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                                        >  2:228143/1‑150 (MQ=255)                                           gACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgtc                                                                                                     <  2:481137/150‑1 (MQ=255)                                            aCGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCa                                                                                                    <  1:1541140/150‑1 (MQ=255)                                             cggcggCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCAt                                                                                                   <  1:399722/150‑1 (MQ=255)                                             cggcggCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCAt                                                                                                   >  2:161582/1‑150 (MQ=255)                                               gcggcTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATcc                                                                                                 <  2:1529206/150‑1 (MQ=255)                                                 ggcTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAg                                                                                               <  2:635558/150‑1 (MQ=255)                                                   cTGGTATGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGca                                                                                             <  1:649645/150‑1 (MQ=255)                                                   cTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGca                                                                                             >  1:570277/1‑150 (MQ=255)                                                      gTACGAATCGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACaa                                                                                          >  2:542443/1‑150 (MQ=255)                                                       tACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAAc                                                                                         >  1:1454256/1‑150 (MQ=255)                                                                    gcGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcc                                                                            <  2:65309/150‑1 (MQ=255)                                                                      gATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccac                                                                          >  1:328455/1‑150 (MQ=255)                                                                         cc‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgca                                                                       >  1:842957/1‑150 (MQ=255)                                                                              aaaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaccc                                                                   >  1:425016/1‑150 (MQ=255)                                                                              aaaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaacc                                                                   <  1:499203/150‑1 (MQ=255)                                                                               aaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaacca                                                                  <  1:912703/150‑1 (MQ=255)                                                                                          cAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATg                                                       >  2:757172/1‑150 (MQ=255)                                                                                             ccGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAg                                                    <  1:662742/150‑1 (MQ=255)                                                                                              cGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAgg                                                   <  1:1168674/150‑1 (MQ=255)                                                                                              cGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAgg                                                   <  1:274356/150‑1 (MQ=255)                                                                                                  gAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAGCCACGCAAACCATGCAGGGAAt                                               >  2:1231839/1‑150 (MQ=255)                                                                                                   aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATa                                              >  1:204475/1‑150 (MQ=255)                                                                                                             ggACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAg                                    <  2:394417/150‑1 (MQ=255)                                                                                                              gACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACCCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGc                                   >  1:838456/1‑150 (MQ=255)                                                                                                                    aGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAAc                             <  2:878770/150‑1 (MQ=255)                                                                                                                        aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCCCTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTa                         >  1:63871/1‑150 (MQ=255)                                                                                                                        aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTa                         >  2:124414/1‑150 (MQ=255)                                                                                                                              aTTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCg                   >  2:1574062/1‑150 (MQ=255)                                                                                                                                tGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGtt                 <  1:572267/150‑1 (MQ=255)                                                                                                                                 gCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGttt                >  1:1132825/1‑150 (MQ=255)                                                                                                                                 gCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGttt                >  1:1211157/1‑150 (MQ=255)                                                                                                                                  cAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATATACATCCCAGCACAAACAGTAAGATGGTTTc               <  2:199914/150‑1 (MQ=255)                                                                                                                                   ggtatcTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCt              >  2:620464/6‑150 (MQ=255)                                                                                                                                      tACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtgtg           <  2:1211157/150‑1 (MQ=255)                                                                                                                                       tcTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACCCCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGc          >  1:502995/2‑150 (MQ=255)                                                                                                                                        cTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGcc         <  2:333481/150‑1 (MQ=255)                                                                                                                                         tGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAggg                                                  >  1:1229199/1‑108 (MQ=255)                                                                                                                                         tGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAggg                                                  <  2:1229199/108‑1 (MQ=255)                                                                                                                                            cGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCg     >  2:653541/1‑150 (MQ=255)                                                                                                                                            cGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCg     >  1:1256851/1‑150 (MQ=255)                                                                                                                                             gCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAATGTATTCATTGATCAGAACTCAACGTAGATCTACTCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCgg    <  2:1210100/150‑1 (MQ=255)                                                                                                                                               tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt  <  2:467492/150‑1 (MQ=255)                                                                                                                                               tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt  <  1:228143/150‑1 (MQ=255)                                                                                                                                               tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt  <  1:1481457/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCC‑AAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTT  >  NC_000913/1840652‑1840942

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