Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,443,502	T→A	5.9%	F6L (TTT→TTA) ‡	PP_2140 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,443,502	0	T	A	5.9%	309.6 / 5.3	99	F6L (TTT→TTA) ‡	PP_2140	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/47);  new base A (4/2);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAA  >  NC_002947/2443354‑2443643                                                                                                                                                     |                                                                                                                                               tGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCttt                                                                                                                                                >  4:113293/1‑149 (MQ=255)   aTCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTTCAGTCCTATGTCGATGGCAAGCCGTCCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGc                                                                                                                                              <  2:123451/149‑1 (MQ=255)    tCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGc                                                                                                                                              <  4:127942/148‑1 (MQ=255)      cGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGACGACGATGGCAAGCCGACCGGCTGGAAGGCATTCGACGACGGCAAAGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGc                                                                                                                                           >  8:282383/1‑149 (MQ=255)           cAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCaa                                                                                                                                      >  8:86463/1‑149 (MQ=255)                gCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCAtt                                                                                                                                 <  3:63977/149‑1 (MQ=255)                  aaGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTa                                                                                                                               >  4:183521/1‑149 (MQ=255)                     acgaaGGAGCAGTATGTGCAGTCCGTGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTa                                                                                                                               <  6:207003/143‑1 (MQ=255)                     aCCAAGGAGCAGTAAGTGCAGTCCGAGGTCGATGGCAAGCCGACCCGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTCTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTa                                                                                                                               >  5:207003/1‑146 (MQ=255)                       cAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTatgatg                                                                                                                          >  3:97893/1‑149 (MQ=255)                       cAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTatgatg                                                                                                                          >  8:51556/1‑149 (MQ=255)                         aGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCa                                                                                                                        <  3:183521/149‑1 (MQ=255)                                 tATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGc                                                                                                                         >  1:229740/1‑140 (MQ=255)                                 tATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGc                                                                                                                         <  2:229740/140‑1 (MQ=255)                                 tATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCg                                                                                                                >  2:179660/1‑149 (MQ=255)                                  aTGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCAAAGCGGccc                                                                                                                                         >  4:334718/1‑123 (MQ=255)                                    gtgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCt                                                                                                             >  4:258553/1‑149 (MQ=255)                                     tgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAATATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCtt                                                                                                            <  2:239067/149‑1 (MQ=255)                                     tgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCtt                                                                                                            <  5:217481/149‑1 (MQ=255)                                      gCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTg                                                                                                           >  8:268055/1‑149 (MQ=255)                                       cAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTg                                                                                                           >  1:320334/1‑148 (MQ=255)                                        aGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCTGCCCCAATCATTTATGATGCAGATATCCGCCTTGTa                                                                                                         <  2:89198/149‑1 (MQ=255)                                        aGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTa                                                                                                         <  7:86463/149‑1 (MQ=255)                                         gTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAg                                                                                                        >  8:168470/1‑149 (MQ=255)                                          tCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGc                                                                                                       >  7:261859/1‑149 (MQ=255)                                                 tCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAATGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                <  5:185729/149‑1 (MQ=255)                                                   gATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                >  4:177636/1‑147 (MQ=255)                                                   gATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                >  2:347552/1‑147 (MQ=255)                                                    aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATTTCCGCCTTGTAGCCTTATGGGAg                                                                                             >  1:44904/1‑149 (MQ=255)                                                    aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                >  7:26015/1‑146 (MQ=255)                                                    aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGa                                                                                              <  8:26015/148‑1 (MQ=255)                                                    aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAg                                                                                             <  1:6986/149‑1 (MQ=255)                                                     tGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAgg                                                                                            <  2:320334/149‑1 (MQ=255)                                                      ggCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCAAATCATTTATGATGCac                                                                                                                       >  2:119507/1‑120 (MQ=255)                                                       gCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAggg                                                                                           >  5:318481/1‑148 (MQ=255)                                                          aGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACt                                                                                       <  8:261859/149‑1 (MQ=255)                                                               aCCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTCGTAGCCTTATGGGAGGGCACTGAAAt                                                                                  <  1:179660/149‑1 (MQ=255)                                                                ccGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATg                                                                                 <  7:168470/149‑1 (MQ=255)                                                                ccGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATg                                                                                 <  3:177636/149‑1 (MQ=255)                                                                     tGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                            <  3:113293/149‑1 (MQ=255)                                                                       gAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑g                                                                          >  5:79426/1‑148 (MQ=255)                                                                         aGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGaa                                                                       <  4:97893/149‑1 (MQ=255)                                                                          ggCATTCTACGACGGCAATGCCTGGAAGGTTAAAGATAAGCGCTGATCCAATGGTGCTGTGATTTGGGGCCGCATTGCGGACCCAATAATTAAAAATGAAGAAAAACGCCTAGTAGCCATATGGAAGGGCAAGGAAATGGCCCA‑GGaa                                                                       >  7:221189/1‑148 (MQ=255)                                                                            cATTCTACGACGGCAATGCCTGGAATTTTGAAGATAAGCGGTGATCCATTTGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGATGGCACTGAAATGTCCCA‑GGAActc                                                                    <  2:44904/149‑1 (MQ=255)                                                                             aTTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAActc                                                                    >  2:91916/1‑148 (MQ=255)                                                                               tCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTa                                                                  >  2:4075/1‑148 (MQ=255)                                                                                 tACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTacac                                                               >  1:133058/1‑149 (MQ=255)                                                                                  acgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAcc                                                              >  7:200296/1‑149 (MQ=255)                                                                                   cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAccc                                                             <  8:200296/149‑1 (MQ=255)                                                                                   cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAccc                                                             <  6:318481/149‑1 (MQ=255)                                                                                   cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAccc                                                             >  2:330959/1‑149 (MQ=255)                                                                                    gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGcc                                                                                                      <  3:237446/108‑1 (MQ=255)                                                                                    gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGcc                                                                                                      >  4:237446/1‑108 (MQ=255)                                                                                    gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCg                                                            >  2:167428/1‑149 (MQ=255)                                                                                      cgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTa                                                          >  1:102386/1‑149 (MQ=255)                                                                                      cgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTa                                                          <  1:167428/149‑1 (MQ=255)                                                                                              ccTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   <  1:347552/148‑1 (MQ=255)                                                                                               cTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTac                                                                 <  7:71727/133‑1 (MQ=255)                                                                                               cTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTac                                                                 >  8:71727/1‑133 (MQ=255)                                                                                                tGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGa                                                 >  3:250477/1‑148 (MQ=255)                                                                                                tGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGa                                                 >  7:171161/1‑148 (MQ=255)                                                                                                 ggAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGaa                                                >  8:275305/1‑148 (MQ=255)                                                                                                  gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGATCTCTACACCCGTACCAACCAGaaa                                               <  7:51556/148‑1 (MQ=255)                                                                                                  gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGaaa                                               <  4:250477/148‑1 (MQ=255)                                                                                                  gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGaaa                                               <  8:171161/148‑1 (MQ=255)                                                                                                  gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAAc                                              >  1:62610/1‑149 (MQ=255)                                                                                                      gTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGt                                           >  7:286214/1‑148 (MQ=255)                                                                                                        tGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   <  6:18584/138‑1 (MQ=255)                                                                                                        tGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAACCGTACCAACCAGagatcggg                                          >  5:18584/1‑140 (MQ=255)                                                                                                            gATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGc                                    >  4:14450/1‑149 (MQ=255)                                                                                                              tAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCg                                  <  2:62610/149‑1 (MQ=255)                                                                                                              tAAGCGCTGATCCAATGGTGCTGTGAATTGGGGCCGCTTTGCGGCCCCAATCATTTAAGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGAACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTATACGCCg                                  >  7:335207/1‑149 (MQ=255)                                                                                                                 gcgcTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   <  2:135549/129‑1 (MQ=255)                                                                                                                 gcgcTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   >  1:135549/1‑129 (MQ=255)                                                                                                                     tGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCg                                  >  3:180145/1‑142 (MQ=255)                                                                                                                     tGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCg                                  <  4:180145/142‑1 (MQ=255)                                                                                                                          cATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       <  3:14450/148‑1 (MQ=255)                                                                                                                          cATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       <  2:102386/148‑1 (MQ=255)                                                                                                                           aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       >  6:63150/1‑147 (MQ=255)                                                                                                                           aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  6:79426/149‑1 (MQ=255)                                                                                                                           aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  2:133058/149‑1 (MQ=255)                                                                                                                           aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  7:275305/149‑1 (MQ=255)                                                                                                                            ttGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  1:91916/149‑1 (MQ=255)                                                                                                                               gTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTc                  <  1:330959/148‑1 (MQ=255)                                                                                                                                tGCTGTGATTTGGTGTCGCTTTGCGGCCCAAATTATTTATGTTGGAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGTAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAAGTGTTCTTCGCGGGCCTCGCCCTCGAGTcc                 <  8:335207/148‑1 (MQ=255)                                                                                                                                   tgtgATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATgg             >  6:162402/1‑149 (MQ=255)                                                                                                                                   aaatgattGGGGCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     >  1:225098/7‑141 (MQ=255)                                                                                                                                   aaatgattGGGGCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  2:225098/135‑1 (MQ=255)                                                                                                                                     tgATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGc            <  3:258553/148‑1 (MQ=255)                                                                                                                                       aTTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATgg             >  1:79042/1‑145 (MQ=255)                                                                                                                                       aTTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATgg             <  2:79042/145‑1 (MQ=255)                                                                                                                                          tGGGGCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGc      >  2:104446/1‑149 (MQ=255)                                                                                                                                            gggCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCt     <  8:286214/148‑1 (MQ=255)                                                                                                                                            gggCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  1:93034/132‑1 (MQ=255)                                                                                                                                            gggCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     >  2:93034/1‑132 (MQ=255)                                                                                                                                             ggCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTg    >  2:4124/1‑148 (MQ=255)                                                                                                                                              gCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGaa  >  5:297590/1‑149 (MQ=255)                                                                                                                                                 gCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       <  6:153422/125‑1 (MQ=255)                                                                                                                                                  cTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑gg                                                                         >  2:330136/1‑74 (MQ=255)                                                                                                                                                  cTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑gg                                                                         <  1:330136/74‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               TGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAA  >  NC_002947/2443354‑2443643

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