Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,632,097	A→G	7.2%	E284G (GAG→GGG)	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,632,097	0	A	G	7.2%	242.4 / 6.1	83	E284G (GAG→GGG)	lepA	elongation factor 4
Reads supporting (aligned to +/- strand):  ref base A (42/35);  new base G (4/2);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

CACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGA  >  NC_002947/1631949‑1632241                                                                                                                                                     |                                                                                                                                                  cACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   <  4:133576/149‑1 (MQ=255)  aCTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAg                                                                                                                                                  >  1:183005/1‑149 (MQ=255)  aCTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAg                                                                                                                                                  >  3:261232/1‑149 (MQ=255)       gTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGa                                                                                                                                            <  2:352647/149‑1 (MQ=255)          gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGT‑CGAgga                                                                                                                                         >  6:106800/1‑148 (MQ=255)           aCAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGgtg                                                                                                                                        <  6:77379/149‑1 (MQ=255)            cAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGT‑CGAggaca                                                                                                                                       <  4:292393/149‑4 (MQ=255)                       tCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCTCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGT‑CGAggacagg                                                                                                                                     <  4:239795/140‑6 (MQ=255)                       tCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGT‑CGAggacagg                                                                                                                                     >  3:239795/1‑135 (MQ=255)                       tCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCtt                                                                                                                             <  4:117120/148‑1 (MQ=255)                         ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGCCCCTCTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCa                                                                                                                           >  3:45793/1‑148 (MQ=255)                         ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTc                                                                                                                            >  4:298496/1‑147 (MQ=255)                            aCCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCaagaa                                                                                                                       >  2:134742/1‑149 (MQ=255)                            aCCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCaaga                                                                                                                        <  3:301711/148‑1 (MQ=255)                              cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCaag                                                                                                                         >  7:62695/1‑145 (MQ=255)                              cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                     <  3:81035/149‑1 (MQ=255)                              cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                     <  3:63166/149‑1 (MQ=255)                                cGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAg                                                                                                                                                  <  8:304641/119‑1 (MQ=255)                                cGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAg                                                                                                                                                  >  7:304641/1‑119 (MQ=255)                                        acCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGCCCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCg                                                                                                             >  2:295687/1‑147 (MQ=255)                                         cccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGATGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAg                                                                                                          <  3:30798/149‑1 (MQ=255)                                           cAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGGATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGAGCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                        <  1:6117/149‑1 (MQ=255)                                                cGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACg                                                                                                   >  8:248317/1‑149 (MQ=255)                                                 gCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGc                                                                                                  >  6:183352/1‑149 (MQ=255)                                                   tGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCg                                                                                                >  5:100164/1‑149 (MQ=255)                                                         gAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                          >  4:26269/1‑149 (MQ=255)                                                         gAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                          >  5:146622/1‑149 (MQ=255)                                                          aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGcc                                                                                             <  4:261232/145‑1 (MQ=255)                                                          aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGtt                                                                                         >  1:47757/1‑149 (MQ=255)                                                             gacGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCT‑‑ACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCgg                                                                                    >  5:53235/3‑149 (MQ=255)                                                             gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTccc                                                                                      <  3:26269/149‑1 (MQ=255)                                                                      gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTc                                                                             <  1:134742/149‑1 (MQ=255)                                                                      gTAGGCTTCATCACCGCCACCATCAAGGACACCCACGGTGCCCCGGTTGGTCACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAACATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGGCAGccc                                                                             <  2:320083/149‑3 (MQ=255)                                                                        aGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                    <  3:220691/108‑1 (MQ=255)                                                                        aGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                    >  4:220691/1‑108 (MQ=255)                                                                        aGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                            >  1:364189/1‑148 (MQ=255)                                                                         ggCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcg                                                                           <  3:298496/148‑1 (MQ=255)                                                                          gCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTccc                                                                                      >  6:257710/1‑136 (MQ=255)                                                                          gCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTccc                                                                                      <  5:257710/136‑1 (MQ=255)                                                                           cTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgac                                                                         >  6:228794/1‑148 (MQ=255)                                                                             tcatcaTCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCgacgac                                                                      >  7:193220/1‑149 (MQ=255)                                                                                tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCg                                                                                                             <  6:112196/107‑1 (MQ=255)                                                                                tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCg                                                                                                             >  5:112196/1‑107 (MQ=255)                                                                                tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCT‑‑ACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  6:53235/149‑1 (MQ=255)                                                                                  atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  8:62695/149‑1 (MQ=255)                                                                                  atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  2:183005/149‑1 (MQ=255)                                                                                   tcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAg                                                                >  7:306153/1‑149 (MQ=255)                                                                                          cATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTGTTGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTcc                                                         <  1:295687/149‑1 (MQ=255)                                                                                          cATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTcc                                                         <  5:228794/149‑1 (MQ=255)                                                                                            tCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCcgc                                                       >  1:204097/1‑149 (MQ=255)                                                                                                      ccACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTg                                              >  7:257330/1‑148 (MQ=255)                                                                                                      ccACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTg                                              >  1:254000/1‑148 (MQ=255)                                                                                                          ggTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGaa                                         <  6:100164/149‑1 (MQ=255)                                                                                                           gTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTGCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAg                                           >  1:82335/1‑147 (MQ=255)                                                                                                            tGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGc                                       <  5:183352/149‑1 (MQ=255)                                                                                                               cccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                    <  2:82335/149‑1 (MQ=255)                                                                                                                ccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAc                                   >  3:116013/1‑149 (MQ=255)                                                                                                                 cGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAcc                                  <  7:248317/149‑1 (MQ=255)                                                                                                                  ggTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAccc                                 <  4:116013/149‑1 (MQ=255)                                                                                                                       gTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGaa                             <  8:306153/148‑1 (MQ=255)                                                                                                                       gTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAAc                            <  7:245731/149‑1 (MQ=255)                                                                                                                           cacCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGAct                        <  8:257330/149‑1 (MQ=255)                                                                                                                           cacCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGAct                        <  2:254000/149‑1 (MQ=255)                                                                                                                            accctgaccctgTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGGCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 >  1:118010/1‑107 (MQ=255)                                                                                                                            accctgaccctgTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  2:118010/107‑1 (MQ=255)                                                                                                                            accctgaccctgTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGaa                             >  6:102692/1‑143 (MQ=255)                                                                                                                            accctgaccctgTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGaa                             <  5:102692/143‑1 (MQ=255)                                                                                                                            accctgaccctgTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActc                       >  8:270004/1‑149 (MQ=255)                                                                                                                                     ctgTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTg                >  3:58751/1‑147 (MQ=255)                                                                                                                                        accTCGACCTCAGGGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTa           >  2:192390/2‑149 (MQ=255)                                                                                                                                          cTCGACCTCAGGGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         >  3:344556/1‑149 (MQ=255)                                                                                                                                          cTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         >  6:47976/1‑149 (MQ=255)                                                                                                                                          cTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         >  2:75756/1‑149 (MQ=255)                                                                                                                                          cTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         >  2:189029/1‑149 (MQ=255)                                                                                                                                          cTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTACGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTAc          <  4:45793/148‑1 (MQ=255)                                                                                                                                           tCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACAt        >  2:18079/1‑149 (MQ=255)                                                                                                                                                cccTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGc      >  1:272214/1‑146 (MQ=255)                                                                                                                                                cccTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGc      <  2:272214/146‑1 (MQ=255)                                                                                                                                                cccTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCg    >  7:294638/1‑148 (MQ=255)                                                                                                                                                cccTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCg    >  4:198653/1‑148 (MQ=255)                                                                                                                                                 ccTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCgg   >  4:355558/1‑148 (MQ=255)                                                                                                                                                  cTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa  >  4:84476/1‑148 (MQ=255)                                                                                                                                                   tGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa  >  6:313504/1‑147 (MQ=255)                                                                                                                                                   tGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa  >  5:42215/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGT‑CGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGA  >  NC_002947/1631949‑1632241

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