Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,030	G→A	7.3%	H64Y (CAC→TAC) ‡	PP_4580 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,030	0	G	A	7.3%	145.5 / 2.2	55	H64Y (CAC→TAC) ‡	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/22);  new base A (2/2);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-02

CTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAG  >  NC_002947/5199887‑5200174                                                                                                                                                |                                                                                                                                                   cTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCAGTCGATTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCCTCGCTGCGTCGGTATtgctg                                                                                                                                                 >  5:363262/1‑147 (MQ=255)  tGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTgcag                                                                                                                                              >  5:160274/1‑149 (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                                                           <  3:68143/149‑1 (MQ=255)                                                                                    tcaGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTgc                                                           >  7:232567/1‑149 (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                                <  1:235524/149‑1 (MQ=255)                                                                                                              cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCtt                                <  3:341910/149‑1 (MQ=255)                                                                                                               tGTTGAACACCGGGTCGTCGCTGGTGGGGTTGTGGTGCAGCCGCGCCGGGAAGTGGGCCATGAAAAAACGCGCGG‑CCGGTTGGCGGGACCCCACGAGCGCCGCGGGCGCTTGGAACAGGTGC‑GCCGCGACCCGCGGGCGATAGtgcgt                                >  2:142499/1‑144 (MQ=255)                                                                                                                 ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTg                               <  5:224377/147‑1 (MQ=255)                                                                                                                 ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTg                               <  5:308176/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)                                                                                                                                                |                                                                                                                                                   CTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAG  >  NC_002947/5199887‑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: ‑