Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,995,396	A→T	8.5%	intergenic (‑117/+33)	PP_3521 ← / ← PP_3522	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,995,396	0	A	T	8.5%	207.7 / 7.9	71	intergenic (‑117/+33)	PP_3521/PP_3522	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (3/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAGGCC  >  NC_002947/3995254‑3995543                                                                                                                                               |                                                                                                                                                    gTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCcgc                                                                                                                                              <  1:709108/150‑1 (MQ=255) gTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCcgc                                                                                                                                              <  1:1517993/150‑1 (MQ=255)  tAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCTCGGGTAAAccc                                                                                                                                                >  2:191760/1‑147 (MQ=255)  tAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAAccc                                                                                                                                                <  1:191760/147‑1 (MQ=255)       gTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCAc                                                                                                                                        <  1:1695560/150‑1 (MQ=255)                  tGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                              <  1:1961384/133‑1 (MQ=255)                  tGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                              >  2:1961384/1‑133 (MQ=255)                     aCTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGc                                                                                                                                    >  2:688747/1‑140 (MQ=255)                     aCTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGc                                                                                                                                    <  1:688747/140‑1 (MQ=255)                       tCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGACTACGGAGCCGGc                                                                                                                        >  2:1152544/1‑150 (MQ=255)                            gtgtGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGCTGGCACCGGCTACGCCGGTGTTCGCTGGTCAACCCGCGCCCACGGGCGACCGAGCCGGCtttcc                                                                                                                   >  2:74021/1‑148 (MQ=255)                            gtgtGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTc                                                                                                                   >  2:1277682/1‑150 (MQ=255)                             tgtgGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCa                                                                                                                  <  1:1068138/150‑1 (MQ=255)                               tgGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAgg                                                                                                                <  2:1495319/150‑1 (MQ=255)                                 gTCACCTCCGCGCGAGGTTGGTGGTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGtt                                                                                                              >  2:1306235/1‑150 (MQ=38)                                    aCCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTg                                                                                                           >  2:1748686/1‑150 (MQ=255)                                     ccTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCTCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGa                                                                                                          <  1:1152544/150‑1 (MQ=255)                                     ccTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGa                                                                                                          >  2:1961086/1‑150 (MQ=255)                                     ccTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGa                                                                                                          <  1:893518/150‑1 (MQ=255)                                       tctGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGc                                                                                                        <  1:845283/150‑1 (MQ=255)                                            gTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTATGCCGGTGTTCGCGGGTTTACCcgcgaac                                                                                                                                          <  1:2267192/111‑4 (MQ=38)                                            gTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTTTACCcgcgaac                                                                                                                                          >  2:2267192/1‑108 (MQ=255)                                                    gttgttATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGCGTAAACCCGCCCCCACGGGCGACCGCGCCGGCTTTTCAGGTTAAGAGCCAGCCggcggcga                                                                                           >  2:124277/1‑148 (MQ=255)                                                    gtggttATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCAGCTACGCCGGTGGTCGCGGAGAAACCCGCGACCTCGGGCGACCGAGCCGCGCCGTCAGATTAAGAGCAAGACGGCGGCTa                                                                                           >  2:171322/4‑150 (MQ=255)                                                                   tttAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGa                                                                                                                              >  2:1338515/1‑100 (MQ=255)                                                                   tttAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGa                                                                                                                              <  1:1338515/100‑1 (MQ=255)                                                                        tCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGc                                                                       >  2:1547391/1‑150 (MQ=255)                                                                         cAGGACCGGCACCGATCTACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCa                                                                      <  1:124277/150‑1 (MQ=255)                                                                         cAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCATTTAAACCCGCCCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCa                                                                      <  1:171322/150‑1 (MQ=255)                                                                                 gCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTcc                                                              >  1:2275030/1‑150 (MQ=255)                                                                                     cGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCTCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTg                                                          >  2:433376/1‑150 (MQ=255)                                                                                            cAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCg                                                   >  1:2257576/1‑150 (MQ=255)                                                                                            cAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCg                                                   <  1:1105619/150‑1 (MQ=255)                                                                                             aaGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCa                                                                      >  2:795717/1‑130 (MQ=255)                                                                                             aaGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCa                                                                      <  1:795717/130‑1 (MQ=255)                                                                                             aaGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTACGGCGGTAAACCCGCGCACGCGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCGAATGCTTCGAGCTCCGCCAGCATGCCGTCCAGGGCCAGCCgg                                                  >  2:905390/1‑150 (MQ=255)                                                                                               ggAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCgggg                                                >  2:1215689/1‑150 (MQ=255)                                                                                                   gaaCCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCAcc                                            >  1:766502/1‑150 (MQ=255)                                                                                                      ccaccGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGt                                         <  1:1547391/150‑1 (MQ=255)                                                                                                       caccGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTa                                        <  1:1277682/150‑1 (MQ=255)                                                                                                       caccGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTa                                        <  1:74021/150‑1 (MQ=255)                                                                                                       caccGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTa                                        >  1:1274581/1‑150 (MQ=255)                                                                                                       caccGTGAGGAGGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTa                                        >  2:478581/1‑150 (MQ=255)                                                                                                            tGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGc                                   <  1:1961086/150‑1 (MQ=255)                                                                                                                  tGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAg                             <  2:2257576/150‑1 (MQ=255)                                                                                                                  tGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAg                             >  1:842357/1‑150 (MQ=255)                                                                                                                   ggCAACGGCTACGCCGGTGTTCGCTGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAgg                            >  2:1955932/1‑150 (MQ=255)                                                                                                                    gCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCTACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGc                                   <  2:1511950/142‑1 (MQ=255)                                                                                                                    gCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCTACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGc                                   >  1:1511950/1‑142 (MQ=255)                                                                                                                     cACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGGACCGGGATCAGCGTGCAGGTc                          >  2:2003985/1‑150 (MQ=255)                                                                                                                     cACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTc                          >  2:544747/1‑150 (MQ=255)                                                                                                                      aCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGt                           >  2:1583148/1‑148 (MQ=255)                                                                                                                      aCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGt                           <  1:1583148/148‑1 (MQ=255)                                                                                                                      aCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCg                         >  1:2015260/1‑150 (MQ=255)                                                                                                                        cggcgtagCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGt                       <  2:160417/143‑1 (MQ=255)                                                                                                                        cGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTGGTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGt                       >  1:1021064/1‑150 (MQ=255)                                                                                                                        cGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGt                       >  2:48348/1‑150 (MQ=255)                                                                                                                          gCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGttt                     <  1:1526695/150‑1 (MQ=255)                                                                                                                               gCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCAACAGACCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCAtgt                <  1:738477/150‑1 (MQ=255)                                                                                                                                       gcgcgGGTTTACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGa                                                                                  <  1:1414644/75‑1 (MQ=255)                                                                                                                                       gcgcgGGTTTACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGa                                                                                  >  2:1414644/2‑76 (MQ=255)                                                                                                                                         gcgGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTccccc      >  2:588639/1‑149 (MQ=255)                                                                                                                                         gcgGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCa      >  1:17609/1‑150 (MQ=255)                                                                                                                                         gcgGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCa      >  1:238863/1‑150 (MQ=255)                                                                                                                                          cgGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAg     <  2:2015260/150‑1 (MQ=255)                                                                                                                                          cgGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAg     <  2:1274581/150‑1 (MQ=255)                                                                                                                                             gTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGtttt                                                                                                            >  2:691673/1‑44 (MQ=255)                                                                                                                                             gTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGtttt                                                                                                            <  1:691673/44‑1 (MQ=255)                                                                                                                                             gTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTc                                                                                    <  1:1447869/68‑1 (MQ=255)                                                                                                                                             gTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTc                                                                                    >  2:1447869/1‑68 (MQ=255)                                                                                                                                             gTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAACCGGGGCCCCAGCGTAAGCAGCCAGTACGGGTGCACGTTCGCTCCCCAGGcc  >  2:1231197/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                    GTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCGGTTTCATGTGCGCTCCCCAGGCC  >  NC_002947/3995254‑3995543

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