Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,763,001	Δ1 bp	7.7%	coding (201/723 nt)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,763,001	0	G	.	7.7%	243.3 / 14.8	78	coding (201/723 nt)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base . (3/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

ACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTC  >  NC_002947/3762853‑3763150                                                                                                                                                     |                                                                                                                                                       aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGACGATTACGCCTACAAGGCGGAAAAGATCGAGCGCGCCTGGGCGCTCGACGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                                       >  1:1867715/1‑150 (MQ=255)  cGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCAGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGt                                                                                                                                                      <  2:1098446/150‑1 (MQ=255)  cGCAGGCGCCCACGGCCTGTGCACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGt                                                                                                                                                      <  1:138889/150‑1 (MQ=255)  cGCAGGCGCCCACGGCCTGCGGACCGAGGAGCGACGTGGCAATATCGAGGCAGTGTGCGACCACGGGCGCTATGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGt                                                                                                                                                      <  1:224894/150‑1 (MQ=255)      ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAg                                                                                                                                                  <  1:459950/150‑1 (MQ=255)           ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                                             >  1:1788015/1‑150 (MQ=255)             aCGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTgg                                                                                                                                           >  2:894977/1‑150 (MQ=255)               ggCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGGGCCAGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑g                                                                                                                                        >  2:1159436/1‑150 (MQ=255)                   tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATc                                                                                                                                    >  2:721026/1‑150 (MQ=255)                   tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATc                                                                                                                                    >  1:1551042/1‑150 (MQ=255)                       gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGc                                                                                                                                >  2:1827976/1‑150 (MQ=255)                           gaggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGtg                                                                                                                            >  2:568086/1‑150 (MQ=255)                                gCGACGTGGCAATATCGAGGCCGTGTACGTCCACGGCGCCGCGGGAAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTCCGACGTGACCCGGAACATGATCTCGGTCACGGACTCTCGCCTGGC‑GGCTCATGCCGGGCTGCa                                                                                                                       >  2:403211/1‑150 (MQ=255)                                gCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCa                                                                                                                       >  2:84226/1‑150 (MQ=255)                                gCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCa                                                                                                                       >  1:372907/1‑150 (MQ=255)                                   aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAAcc                                                                                                                    <  1:1807226/150‑1 (MQ=255)                                    cGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCa                                                                                                                   <  1:1071589/150‑1 (MQ=255)                                        gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGCAAACCGAATGCCGGTCCTGGTCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTa                                                                                                               <  1:431508/150‑1 (MQ=255)                                          aaTATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGTC‑GGATCATGCCGTGGTGCAACCACTTAcc                                                                                                             >  1:1219960/1‑150 (MQ=255)                                                 aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                                                      >  1:947821/1‑150 (MQ=255)                                                 aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                                                      >  2:1408472/1‑150 (MQ=255)                                                    ccGTGTACGGCCACGGCGCCGAGGCTAACGCCTACACGGCGGAAAAGATCAAACGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg                                                                                                   >  1:1184737/1‑150 (MQ=255)                                                            ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGcc                                                                                           <  2:1867715/150‑1 (MQ=255)                                                              ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                                         <  2:1016965/150‑1 (MQ=255)                                                                           gATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAc                                                                            >  2:1104121/1‑150 (MQ=255)                                                                           gATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAc                                                                            >  1:1747646/1‑150 (MQ=255)                                                                             tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCGACCACTTACCCCGCCTGCGGCa                                                                                                 >  1:1877501/1‑127 (MQ=255)                                                                             tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCGACCACTTACCCCGCCTGCGGCa                                                                                                 <  2:1877501/127‑1 (MQ=255)                                                                              aaCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACggg                                                                         <  1:721026/150‑1 (MQ=255)                                                                               aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGc                                                                        >  1:608109/1‑150 (MQ=255)                                                                                cGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGcc                                                                       >  1:75231/1‑150 (MQ=255)                                                                                 gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCt                                                                      >  2:1254666/1‑150 (MQ=255)                                                                                  ccTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTgctgct                                                                                   >  2:265011/1‑136 (MQ=255)                                                                                  ccTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTgctgct                                                                                   <  1:265011/136‑1 (MQ=255)                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCAc                                                                                                                  <  2:847753/104‑1 (MQ=255)                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCAc                                                                                                                  >  1:847753/1‑104 (MQ=255)                                                                                   cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                                    <  2:840958/150‑1 (MQ=255)                                                                                      cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCGAGC‑GGCTCATGCCGTGCTGCAAGCATTGTCGCTGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGTTc                                                                 >  1:1377546/1‑150 (MQ=255)                                                                                       aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                                                >  2:603076/1‑150 (MQ=255)                                                                                         ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                              <  2:372907/150‑1 (MQ=255)                                                                                           cGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                                                            <  2:1702225/150‑1 (MQ=255)                                                                                               aaaGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                                                        >  2:1393589/1‑150 (MQ=255)                                                                                                   aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgc                                                                                                                                                 >  2:417546/1‑58 (MQ=255)                                                                                                   aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgc                                                                                                                                                 <  1:417546/58‑1 (MQ=255)                                                                                                   aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTgg                                                                                                                                           >  1:1804126/1‑64 (MQ=255)                                                                                                   aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTgg                                                                                                                                           <  2:1804126/64‑1 (MQ=255)                                                                                                   aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGAc                                                    >  1:1784843/1‑150 (MQ=255)                                                                                                     cAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCAGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                                                  <  2:264280/150‑1 (MQ=255)                                                                                                           cgcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                                            <  2:1551042/150‑1 (MQ=255)                                                                                                           cgcCTGGGCGTTCGACGTCACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                                            <  2:576491/150‑1 (MQ=255)                                                                                                            gccgcggCGTTCGACGTGACCGGGACCATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACTACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGACACAGCGCCCGGACAAGCAATgg                                           <  2:947821/145‑1 (MQ=255)                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑gg                                                                                                                                       >  2:826975/1‑58 (MQ=255)                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑gg                                                                                                                                       <  1:826975/58‑1 (MQ=255)                                                                                                            gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                                           >  1:732149/1‑150 (MQ=255)                                                                                                              cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                                         <  1:1393589/150‑1 (MQ=255)                                                                                                              cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                                         <  2:732149/150‑1 (MQ=255)                                                                                                              cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAG‑AACTGTCTGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                                         <  1:593791/150‑1 (MQ=255)                                                                                                               tGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTATCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTc                                        >  2:1183976/1‑150 (MQ=255)                                                                                                                aggCGCTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                                             >  1:1380978/2‑49 (MQ=255)                                                                                                                aggCGCTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                                             <  2:1380978/48‑1 (MQ=255)                                                                                                                   cGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATc                                    >  2:1862255/1‑150 (MQ=255)                                                                                                                     ttCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg                                  <  1:1506401/150‑1 (MQ=255)                                                                                                                          cGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCg                             <  1:1159436/150‑1 (MQ=255)                                                                                                                             gACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCCCTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc                          <  1:259543/150‑1 (MQ=255)                                                                                                                              aCCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCCGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                                            <  2:1141565/82‑1 (MQ=21)                                                                                                                              aCCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCCGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                                            >  1:1141565/1‑82 (MQ=21)                                                                                                                                cGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCa                       >  2:1314780/1‑150 (MQ=255)                                                                                                                                      cATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTATCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCaa                 <  1:1183976/150‑1 (MQ=255)                                                                                                                                      cATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCaa                 <  1:821553/150‑1 (MQ=255)                                                                                                                                       aTGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg                <  1:1104121/150‑1 (MQ=255)                                                                                                                                        tGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAgg               <  2:709993/150‑1 (MQ=255)                                                                                                                                          aTCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGAt             >  2:284751/1‑150 (MQ=255)                                                                                                                                          aTCCCGGTCACGGTCGAGCGCCCAGA‑GGATCATTCCGGTCGTCAACCGCGTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGAt             <  2:1451884/150‑1 (MQ=255)                                                                                                                                             ccGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGAGTCAATCAGGATCGCACCCATGTCAAGGATGCg          >  1:1381024/1‑150 (MQ=255)                                                                                                                                              cGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTCCCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACACATGTCAAGGATGCGa         <  2:842698/150‑1 (MQ=255)                                                                                                                                                gTCAC‑GTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                                       <  1:1696932/69‑1 (MQ=255)                                                                                                                                                gTCAC‑GTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                                       <  1:1697914/69‑1 (MQ=255)                                                                                                                                                gTCAC‑GTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                                       >  2:1697914/1‑69 (MQ=255)                                                                                                                                                gTCAC‑GTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgc                                                                                       >  2:1696932/1‑69 (MQ=255)                                                                                                                                                     ggTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTc  <  1:894977/150‑1 (MQ=255)                                                                                                                                                     ggTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTc  >  2:922797/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTC  >  NC_002947/3762853‑3763150

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