Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,024	G→A	6.6%	R66C (CGC→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,024	0	G	A	6.6%	170.7 / 2.2	61	R66C (CGC→TGC)	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (2/2);  total (33/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 = 4.23e-01

GGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGC  >  NC_002947/5199880‑5200171                                                                                                                                                 |                                                                                                                                                      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)                                                                                                                                                 |                                                                                                                                                      GGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGC  >  NC_002947/5199880‑5200171

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