Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,044,765	T→A	6.5%	L158Q (CTG→CAG)	PP_0904 →	protein InaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,044,765	0	T	A	6.5%	221.4 / 3.4	76	L158Q (CTG→CAG)	PP_0904	protein InaA
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base A (2/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCG  >  NC_002947/1044617‑1044913                                                                                                                                                     |                                                                                                                                                     gAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCt                                                                                                                                                      <  2:104774/149‑1 (MQ=255)   aGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCCGTGCTGCCAGt                                                                                                                                                    <  8:216113/149‑1 (MQ=255)   aGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCCGTGCTGCCAGt                                                                                                                                                    <  5:310022/149‑1 (MQ=255)    gcgcTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTTGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTa                                                                                                                                                   <  7:90494/149‑1 (MQ=255)          gACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCCGTGCTGCCAGt                                                                                                                                                    >  8:205074/1‑142 (MQ=255)          gACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCCGTGCTGCCAGt                                                                                                                                                    <  7:205074/142‑1 (MQ=255)          gACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCaa                                                                                                                                             >  5:7911/1‑149 (MQ=255)           aCGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAg                                                                                                                                            <  6:7911/149‑1 (MQ=255)           aCGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAg                                                                                                                                            >  1:364964/1‑149 (MQ=255)             ggTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCa                                                                                                                                          <  5:189732/149‑1 (MQ=255)                  tGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTa                                                                                                                                      <  8:279160/148‑1 (MQ=255)                  tGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCAGTATGGCAAGCACGTAt                                                                                                                                     >  7:139898/1‑149 (MQ=255)                   gTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTAt                                                                                                                                     <  6:387133/148‑1 (MQ=255)                      gAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATtcat                                                                                                                                 >  2:8494/1‑149 (MQ=255)                         cTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATtcatc                                                                                                                                >  3:31578/1‑147 (MQ=255)                          tCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAg                                                                                                                             >  3:279374/1‑149 (MQ=255)                            gACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCt                                                                                                                                                      <  4:47113/122‑1 (MQ=255)                            gACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCt                                                                                                                                                      >  3:47113/1‑122 (MQ=255)                            gACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAg                                                                                                                             <  8:180957/147‑1 (MQ=255)                             acacCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGt                                                                                                                           >  2:329786/1‑148 (MQ=255)                               acCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAAt                                                                                                                        >  4:206818/1‑149 (MQ=255)                                 cTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTg                                                                                                                      >  2:269482/1‑149 (MQ=255)                                    gCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGc                                                                                                                    >  7:27142/1‑148 (MQ=255)                                         ccGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGc                                                                                                              <  3:206818/149‑1 (MQ=255)                                         ccGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGc                                                                                                              >  5:17381/1‑149 (MQ=255)                                         ccGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGc                                                                                                              >  8:218273/1‑149 (MQ=255)                                          cGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCTTGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCg                                                                                                             <  7:128881/149‑1 (MQ=255)                                            aaGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAg                                                                                                           <  1:329786/149‑1 (MQ=255)                                                  cccGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACTTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGcc                                                                                                     >  4:72010/1‑149 (MQ=255)                                                  cccGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGcc                                                                                                     <  1:8494/149‑1 (MQ=255)                                                   ccGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGCTGCCCCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAATGTAATTGGCGAAGGCGAGCCGGccc                                                                                                    >  5:155827/1‑149 (MQ=255)                                                       gAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGt                                                                                                <  5:47184/149‑1 (MQ=255)                                                                ccGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCAc                                                                                                                                         <  6:204740/99‑1 (MQ=255)                                                                ccGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCAc                                                                                                                                         >  5:204740/1‑99 (MQ=255)                                                                      gtggtgCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCgctgct                                                                                 <  4:279374/149‑1 (MQ=255)                                                                         gtgCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGa                                                                              >  1:349861/1‑149 (MQ=255)                                                                           gCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCAGCGTCGAAGTGGCGCTGCTGGccc                                                                            >  2:125647/1‑149 (MQ=255)                                                                             aTGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTg                                                                          >  3:202512/1‑149 (MQ=255)                                                                               gAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTgg                                                                         >  7:207589/1‑148 (MQ=255)                                                                                  cgcATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGaaaa                                                                     <  6:75197/149‑1 (MQ=255)                                                                                  cgcATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGaaaa                                                                     <  2:349861/149‑1 (MQ=255)                                                                                    cATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGt                                                                   >  6:125348/1‑149 (MQ=255)                                                                                      tGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTg                                                                  >  6:180477/1‑148 (MQ=255)                                                                                      tGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGc                                                                 <  6:17381/149‑1 (MQ=255)                                                                                       gCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTg                                                                  >  2:380756/1‑147 (MQ=255)                                                                                        cTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAg                                                               >  7:163435/1‑149 (MQ=255)                                                                                        cTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAg                                                               <  7:296877/149‑1 (MQ=255)                                                                                         tCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAgg                                                              >  8:72362/1‑149 (MQ=255)                                                                                         tCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAgg                                                              >  3:255738/1‑149 (MQ=255)                                                                                            aGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCATGCTCGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGTCCCGCGTCTAAGTGGCGCTGCTGGACCTGGAAAAGTGCAggcgg                                                           <  1:125647/149‑1 (MQ=255)                                                                                             ggACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAggcggc                                                          <  8:207589/149‑1 (MQ=255)                                                                                               aCCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGgcggcg                                                         >  6:195115/1‑148 (MQ=255)                                                                                                 cTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTAt                                                      >  2:121446/1‑149 (MQ=255)                                                                                                  tGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAg                                                               <  5:267050/139‑1 (MQ=255)                                                                                                  tGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAg                                                               >  6:267050/1‑139 (MQ=255)                                                                                                      ggACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCt                                                 >  5:283542/1‑149 (MQ=255)                                                                                                      ggACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGGGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCt                                                 <  5:225399/149‑1 (MQ=255)                                                                                                       gACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTg                                                >  5:386486/1‑149 (MQ=255)                                                                                                                  gcgGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGaaaa                                                                     <  4:356769/117‑1 (MQ=255)                                                                                                                  gcgGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGaaaa                                                                     >  3:356769/1‑117 (MQ=255)                                                                                                                  gcgGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAgcgggcgg                                     <  1:22474/149‑1 (MQ=255)                                                                                                                   cgGATGCACCTGGGGCACTGGCAGCACGGCTGCCCGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGggcggc                                    <  8:139898/149‑1 (MQ=255)                                                                                                                     gATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGcggccg                                  <  3:72010/149‑1 (MQ=255)                                                                                                                     gATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGcggccg                                  <  4:255738/149‑1 (MQ=255)                                                                                                                      aTGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACTTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGcggccgg                                 <  1:269482/149‑1 (MQ=255)                                                                                                                         cACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCaa                              >  1:205145/1‑149 (MQ=255)                                                                                                                          aCCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAAt                             <  7:107064/149‑1 (MQ=255)                                                                                                                          aCCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAAt                             <  4:202512/149‑1 (MQ=255)                                                                                                                          aCCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAAt                             >  3:378916/1‑149 (MQ=255)                                                                                                                              ggggCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACt                         >  6:51073/1‑149 (MQ=255)                                                                                                                               gggCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGACGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGAc                          <  5:125348/147‑1 (MQ=255)                                                                                                                               gggCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACtt                        <  1:380756/149‑1 (MQ=255)                                                                                                                                 gCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGc                      >  4:129371/1‑149 (MQ=255)                                                                                                                                   aCTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGt                    <  6:386486/149‑1 (MQ=255)                                                                                                                                       gcagcaCGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGc                >  1:103885/1‑149 (MQ=255)                                                                                                                                          gcaCGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGTCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTAc             <  8:27142/149‑1 (MQ=255)                                                                                                                                           caCGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGTAGGCCCGCGTCGAAGTTGAGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTAcg            >  1:295543/1‑149 (MQ=255)                                                                                                                                                cTGCCTGTATTGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACgccgcc       <  4:31578/149‑1 (MQ=255)                                                                                                                                                    cTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGCCTTGCGTCAGCTACGCCGCCAGtc   >  8:361581/1‑149 (MQ=255)                                                                                                                                                     tGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATtcg  >  3:43706/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGGCTGCCTGTATGGCAAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCG  >  NC_002947/1044617‑1044913

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: ‑