Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,032	T→C	7.4%	Q63R (CAG→CGG)	PP_4580 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,032	0	T	C	7.4%	151.2 / 2.7	55	Q63R (CAG→CGG)	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base C (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 = 8.28e-02

CTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAGCAGGC  >  NC_002947/5199887‑5200179                                                                                                                                                  |                                                                                                                                                      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)                                                         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                                                  >  6:308176/1‑149 (MQ=255)                                                                                                 gCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCt                                                  >  5:73410/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)                                                                                                                                                 ctgcAGCAGCGCCTGGAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAGCAgg   <  5:235168/148‑1 (MQ=255)                                                                                                                                                 ctgcAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAGCAGGc  <  5:311903/149‑1 (MQ=255)                                                                                                                                                  tgcAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCt                                              >  3:223328/1‑104 (MQ=255)                                                                                                                                                  tgcAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCt                                              <  4:223328/104‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGG‑CCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGC‑GCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAGCAGGC  >  NC_002947/5199887‑5200179

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