Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,026	T→C	6.5%	Y65C (TAC→TGC)	PP_4580 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,026	0	T	C	6.5%	175.3 / 2.3	62	Y65C (TAC→TGC)	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base C (2/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

GGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAG  >  NC_002947/5199880‑5200174                                                                                                                                                   |                                                                                                                                                       ggTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAg                                                                                                                                                      <  4:270029/149‑1 (MQ=255)        cTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCAGTCGATTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCCTCGCTGCGTCGGTATtgctg                                                                                                                                                 >  5:363262/1‑147 (MQ=255)         tGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTgcag                                                                                                                                              >  5:160274/1‑149 (MQ=255)             ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTa                                                                                                                                                       <  5:38218/136‑1 (MQ=255)             ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTa                                                                                                                                                       >  6:38218/1‑136 (MQ=255)             ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCagcagc                                                                                                                                          >  4:68143/1‑149 (MQ=255)             ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCAGCACTACCGcagcag                                                                                                                                           <  6:328951/148‑1 (MQ=255)              cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCAGCACTACCGCAGCAgcg                                                                                                                                         <  7:232837/149‑1 (MQ=255)               aGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcgc                                                                                                                                        <  3:83925/149‑1 (MQ=255)                    aTGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGc                                                                                                                                    >  4:119080/1‑148 (MQ=255)                                 ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctg                                                                                                                                                 <  1:309888/122‑1 (MQ=255)                                 ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctg                                                                                                                                                 >  2:309888/1‑122 (MQ=255)                                    gCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGaaaaa                                                                                                                   <  1:158344/149‑1 (MQ=255)                                     cccATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAg                                                                                                                  <  5:10866/149‑1 (MQ=255)                                        aTGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGaaaaa                                                                                                                   <  8:319347/145‑1 (MQ=255)                                        aTGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTAGCCCGTGAAAAAgag                                                                                                                >  7:319347/1‑146 (MQ=255)                                                 gcagcaCCTGCCATTCGTTTTCCAGCAGTGCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGGGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGt                                                                                                      <  6:363262/148‑1 (MQ=255)                                                             aTTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtg                                                                                                                                                    >  2:93477/1‑91 (MQ=255)                                                             aTTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtg                                                                                                                                                    <  1:93477/91‑1 (MQ=255)                                                                cGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATc                                                                                       <  6:234779/148‑1 (MQ=255)                                                                 gTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTTCCGGAACCCATc                                                                                       >  7:168879/1‑147 (MQ=255)                                                                                    cGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGa                                                                   >  2:95908/1‑148 (MQ=255)                                                                                    cGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCATCGGGCGTAGTGCAGCAGCAGCGCCTGCAAGTCGCCCATGAAAAAACCCACCG‑CCGTTTGCCCGAAACCCACGAGCGCCACGCCCGAATAGaa                                                                  >  4:341910/1‑149 (MQ=255)                                                                                     ggATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGaa                                                                  >  5:356304/1‑148 (MQ=255)                                                                                      gATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGGTGCGGTATTGCTGCAGCAGCGCCGGCAAGTCGCCCGTGTAAAAGCGCTCGG‑GCGGTTTGCGGCACCCCTAGAGCGCCACGCGCGATTCGcaca                                                                >  6:224377/1‑149 (MQ=255)                                                                                           tcaGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTgc                                                           >  7:232567/1‑149 (MQ=255)                                                                                           tcaGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTgc                                                           <  3:68143/149‑1 (MQ=255)                                                                                                 tttCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCG‑CCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCCCGATTCGAACAGGTGC‑GCCTCGa                                                   >  6:83747/1‑149 (MQ=255)                                                                                                   tCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGAc                                                  >  6:311903/1‑149 (MQ=255)                                                                                                   tCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCATGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCg                                                    <  8:342592/147‑1 (MQ=255)                                                                                                      aTGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTcgcgc                                                                          <  7:275302/123‑1 (MQ=255)                                                                                                      aTGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCGTGAAAAAGCCCTCGG‑CCCTGTGCCGGACCCCATCGCGCGccgcccgc                                                                          >  8:275302/1‑123 (MQ=255)                                                                                                        gCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCt                                             >  5:73410/1‑149 (MQ=255)                                                                                                        gCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCt                                             >  6:308176/1‑149 (MQ=255)                                                                                                              ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCGTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGCACAGGTGC‑GCCTCGACCCGCGGGCta                                        >  7:324280/1‑148 (MQ=255)                                                                                                                 aTCACTGTTGAACACCGGGTCGTCGCTGGTGCGGTTGTGGTGCAGCCGCGCCGGCAAGTAGCCCTTGAAAAAGCGCTCGG‑CCGTTTTGCGGAACCCAACGAGCGCCTCGCGCGATTAGAACAGGTGC‑GCCTCGACCCGCTGGGTATTGt                                    >  2:235524/1‑149 (MQ=255)                                                                                                                 aTCACTGTTGAACACCGGGTCGTCCCAGCGGCGGGAGTGGTGCAGCCCCGCCGGGAAGTAGCCCCTGAAAAAACCCGCGGCCCGTTTG‑CGGGACCCAAACAGCCCCACCCGCGCTTCCAACAGCGGC‑GCCACGCCCCCCAGGCTATAg                                     >  6:235168/1‑148 (MQ=255)                                                                                                                     cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCtt                                <  3:341910/149‑1 (MQ=255)                                                                                                                     cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCtt                                <  1:235524/149‑1 (MQ=255)                                                                                                                      tGTTGAACACCGGGTCGTCGCTGGTGGGGTTGTGGTGCAGCCGCGCCGGGAAGTGGGCCATGAAAAAACGCGCGG‑CCGGTTGGCGGGACCCCACGAGCGCCGCGGGCGCTTGGAACAGGTGC‑GCCGCGACCCGCGGGCGATAGtgcgt                                >  2:142499/1‑144 (MQ=255)                                                                                                                        ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTg                               <  5:308176/147‑1 (MQ=255)                                                                                                                        ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTg                               <  5:224377/147‑1 (MQ=255)                                                                                                                        ttGAACACCGGGTCGTCCCTGCAGCGGTAGTGCTGCAGCAGCGCCAGCAAGTCGCCCATGAAAAAACGCTCGG‑CCGTTTGCCGGAACCCATAGAACGCCTCGCGCGATTAGAAAAGGTGC‑GCCTCAACCCCCAGGCTATAGTCCTTGa                              >  4:325507/1‑148 (MQ=255)                                                                                                                          gAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGGTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCGCGACCCGCTGGCGATAGTGCTTGACCa                           >  2:107998/1‑149 (MQ=255)                                                                                                                          gAACACCGGGTCGTCGCTGCGGCGGTAGTGCTGCAGCCGCGCCTGCAAGTCGCCCTTTAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCAACGAGCGCCTCCCGCCATTCGAACAGGTTC‑GCCTAGACCCCCTGGCTATAGTCCTTGACCa                           >  5:94631/1‑149 (MQ=255)                                                                                                                             cacCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGcc                        <  1:142499/149‑1 (MQ=255)                                                                                                                               ccGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGcccc                      <  1:82583/149‑1 (MQ=255)                                                                                                                                cGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCATGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCATGAACAGccccc                     >  8:251344/1‑148 (MQ=255)                                                                                                                                  ggTCGTCGCTGCTGCGGTAGTGCTGCAGCACCGCCTGCAAGTCGCCCATGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg                    >  2:284370/1‑148 (MQ=255)                                                                                                                                   gcaggcgCTGCTGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGA‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGa                  >  4:305952/6‑149 (MQ=255)                                                                                                                                    tcgtcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc                 <  3:119080/149‑1 (MQ=255)                                                                                                                                         gctgctGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTgg             <  7:251344/148‑1 (MQ=255)                                                                                                                                         gctgctGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa            <  6:160274/149‑1 (MQ=255)                                                                                                                                         gctgctGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa            >  5:268961/1‑149 (MQ=255)                                                                                                                                         gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGAATCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTgg             >  2:100682/1‑148 (MQ=255)                                                                                                                                            gctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGAt           <  4:323646/147‑1 (MQ=255)                                                                                                                                            gctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGAt           >  3:323646/1‑147 (MQ=255)                                                                                                                                             ctggGGTTGTGCTGATGAAGCGCCTGCATGTCGCCCTTCAAAAAGCGCTCGG‑CCGTTTCCCGGAACCCATCTCGCGCCTCGCGCTATTTGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATa          <  1:107998/147‑1 (MQ=255)                                                                                                                                              tgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg       >  8:145889/1‑149 (MQ=255)                                                                                                                                                cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCAGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg       <  1:284370/148‑1 (MQ=255)                                                                                                                                                cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCgg      <  8:168879/148‑1 (MQ=255)                                                                                                                                                 ggTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGc     >  6:7409/1‑148 (MQ=255)                                                                                                                                                   tAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCcag  >  8:223679/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAG  >  NC_002947/5199880‑5200174

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