Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,454	T→A	7.0%	S136S (TCA→TCT)	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,454	0	T	A	7.0%	203.6 / 4.9	71	S136S (TCA→TCT)	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (2/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAA  >  NC_002947/1017306‑1017594                                                                                                                                                     |                                                                                                                                             gttaTCATTCTGATCTGGTTCTTTATGCATTGAGTTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACTTCTCTTCGCTACAAAAGGCGTGTGCAATGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgat                                                                                                                                              <  2:268225/146‑1 (MQ=255) gTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgat                                                                                                                                              <  5:337679/149‑1 (MQ=255)           tGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCa                                                                                                                                      >  7:32692/1‑147 (MQ=255)                tGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCt                                                                                                                               >  6:170563/1‑149 (MQ=255)                 ggTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTc                                                                                                                              >  3:318574/1‑149 (MQ=255)                  gTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCGTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGAAGAACACCAAAGGCCTc                                                                                                                              <  4:80160/148‑1 (MQ=255)                       ttATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACa                                                                                                                        >  7:91597/1‑149 (MQ=255)                         aTGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACAtt                                                                                                                      >  2:266730/1‑149 (MQ=255)                          tGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACAtt                                                                                                                      <  7:265692/148‑1 (MQ=255)                            cATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCGACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCATGAGATAAATAAATCCAGTTTGGTGTTCTGATGAACACAAAACGCCACGGGACATTttc                                                                                                                   >  3:227479/1‑149 (MQ=255)                             aTTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTttct                                                                                                                  >  4:62875/1‑149 (MQ=255)                              ttGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTttctt                                                                                                                 >  8:21514/1‑149 (MQ=255)                                gagaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTttcttc                                                                                                                <  2:333785/148‑1 (MQ=255)                                 agaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCa                                                                                                                                      <  3:224255/125‑1 (MQ=255)                                 agaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCa                                                                                                                                      >  4:224255/1‑125 (MQ=255)                                            gCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTg                                                                                                   >  3:63700/1‑149 (MQ=255)                                                   tGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGt                                                                                                  <  3:181399/143‑1 (MQ=255)                                                   tGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGt                                                                                                  >  4:181399/1‑143 (MQ=255)                                                     cGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGAttttt                                                                                          <  8:56496/149‑1 (MQ=255)                                                          cTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGAAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGATCTGATGAACAACAAAGGCCTCGTGACAATTTCTTCGGTATATAAATTGTAGATTTTTgtgg                                                                                      >  8:100799/1‑148 (MQ=255)                                                           ttttGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTgggg                                                                                      <  4:339021/147‑1 (MQ=255)                                                           ttttGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTgggg                                                                                      >  3:339021/1‑147 (MQ=255)                                                                  aCAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTCAGATCAATAAATCCAGTTAGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGc                                                                              >  1:229636/1‑148 (MQ=255)                                                                   cAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGc                                                                              <  8:91597/147‑1 (MQ=255)                                                                     aCTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAAt                                                                          <  4:230339/149‑1 (MQ=255)                                                                     aCTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAAt                                                                          >  3:176803/1‑149 (MQ=255)                                                                      ctctTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATg                                                                         <  3:290239/149‑1 (MQ=255)                                                                       tctTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgt                                                                                  >  3:34571/1‑139 (MQ=255)                                                                       tctTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgt                                                                                  <  4:34571/139‑1 (MQ=255)                                                                           ctcaaCAAACGGCGTGTGCAAGGCTGGGCCCGGATGTCGGGCTTGATATCAATAAATCCAGTTTGGTTTTCTGATCAACACAAAAGGCCTCGTGACATTTTCTTCGTTATATAAATTGTAGAt                                                                                              <  6:205102/119‑1 (MQ=255)                                                                           cGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGt                                                                     <  7:21514/148‑1 (MQ=255)                                                                                cAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAttttt                                                               <  7:100971/149‑1 (MQ=255)                                                                                 aaaaGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgt                                                                                  >  8:32922/1‑129 (MQ=255)                                                                                 aaaaGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgt                                                                                  <  7:32922/129‑1 (MQ=255)                                                                                 aaaaGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTc                                                              >  1:75074/1‑149 (MQ=255)                                                                                 aaaaGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTc                                                              >  4:127244/1‑149 (MQ=255)                                                                                  aaaGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCg                                                             <  7:141320/149‑1 (MQ=255)                                                                                         tgtgCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCggg                                                      >  6:164497/1‑149 (MQ=255)                                                                                          gtgCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                     >  7:70/1‑149 (MQ=255)                                                                                              aaGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTg                                                  >  2:99324/1‑148 (MQ=255)                                                                                               aGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                     <  3:306855/144‑1 (MQ=255)                                                                                               aGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                     >  4:306855/1‑144 (MQ=255)                                                                                                ggCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAg                                               <  5:164497/149‑1 (MQ=255)                                                                                                   aGGGCCCGGATGTCGTGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGCTGATCACCAAAGGGCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGGCTGTAGata                                            <  2:229636/149‑1 (MQ=255)                                                                                                   aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGc                                                                       >  5:123129/1‑122 (MQ=255)                                                                                                   aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGc                                                                       <  6:123129/122‑1 (MQ=255)                                                                                                     ggCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGatat                                           >  4:284127/1‑148 (MQ=255)                                                                                                        ccGGATGTCGGGCTTGATTTCAATAAATGCAGTTTTGTGTTCTGATTTACACCAAAGGCCTCGTTTCATTTTCTTCGGTTTATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGt                                       <  7:100799/149‑1 (MQ=255)                                                                                                        ccGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAg                                        <  6:58726/148‑1 (MQ=255)                                                                                                        ccGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAg                                        >  5:58726/1‑148 (MQ=255)                                                                                                        ccGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGt                                       >  3:180273/1‑149 (MQ=255)                                                                                                          ggATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGc                                     >  2:140039/1‑149 (MQ=255)                                                                                                          ggATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGc                                     <  4:318574/149‑1 (MQ=255)                                                                                                           gATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGATTCGGGCCTGTAGATATGAGTGc                                     >  3:274517/1‑148 (MQ=255)                                                                                                                 gtgCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAg                               <  4:63700/146‑1 (MQ=255)                                                                                                                   gCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTaa                            <  5:170563/149‑1 (MQ=255)                                                                                                                       gagaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGa                                         >  2:134010/1‑132 (MQ=255)                                                                                                                       gagaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGa                                         <  1:134010/132‑1 (MQ=255)                                                                                                                       gagaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATg                          <  2:75074/147‑1 (MQ=255)                                                                                                                           tCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTa                     <  3:62875/148‑1 (MQ=255)                                                                                                                           tCAATAAATCCAGTTTGGTGTTCTGATGAAAACAAAAGGCATCGTGACATTTTCTTCGGTATATAAATTGTATATTTTTGGGGGAGTTTGCCAATGGAGTATTTTTCGATGCGGGCCTGTAGATATGAGTGATCATATTAATGGaatt                     >  7:50533/1‑145 (MQ=255)                                                                                                                              aTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATATTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCa                  >  2:276073/1‑148 (MQ=255)                                                                                                                                 aaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTTTTTATATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAat               <  4:227479/148‑1 (MQ=255)                                                                                                                                     cAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCTGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatctatc          <  4:180273/149‑1 (MQ=255)                                                                                                                                     cAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatctatc          <  3:284127/149‑1 (MQ=255)                                                                                                                                        tttGGTGTTCATCAGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATg                          >  5:74043/1‑130 (MQ=255)                                                                                                                                        tttGGTGTTCATCAGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATg                          <  6:74043/130‑1 (MQ=255)                                                                                                                                          tGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGc      <  3:127244/148‑1 (MQ=255)                                                                                                                                             ttttCTGATTAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCaaaa  <  8:32692/147‑1 (MQ=255)                                                                                                                                             tgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCaaaa  <  1:140039/149‑1 (MQ=255)                                                                                                                                              tttCTGATGAACAACAAAGTCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCaaaa  <  4:176803/147‑1 (MQ=255)                                                                                                                                               ttCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGa                                         >  4:111669/1‑108 (MQ=255)                                                                                                                                               ttCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGa                                         <  3:111669/108‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAA  >  NC_002947/1017306‑1017594

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