Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,024	G→A	7.1%	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	7.1%	194.2 / 3.9	70	R66C (CGC→TGC)	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/41);  new base A (3/2);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-02

GCGGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199877‑5200169                                                                                                                                                    |                                                                                                                                                   gCGGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCgg                                                                                                                                                   >  7:100836/1‑149 (MQ=255)  cGGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCgg                                                                                                                                                   <  7:4929/148‑1 (MQ=255)  cGGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTAGAACACCGGGTCGTCGCTGCTGCg                                                                                                                                                    <  4:254553/147‑1 (MQ=255)   gggTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTa                                                                                                                                                 <  7:42926/149‑1 (MQ=255)   gggTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTa                                                                                                                                                 <  5:73380/149‑1 (MQ=255)     gTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGt                                                                                                                                               <  7:173661/149‑1 (MQ=255)      tCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtg                                                                                                                                              <  4:219475/149‑1 (MQ=255)        ggCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgct                                                                                                                                            <  3:57198/149‑1 (MQ=255)         gCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctg                                                                                                                                           >  3:239522/1‑149 (MQ=255)         gCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctg                                                                                                                                           <  3:146394/149‑1 (MQ=255)         gCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATAATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTAGCATCTTCGGTaattaat                                                                                                                                           >  6:261584/1‑143 (MQ=255)          ccTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctgc                                                                                                                                          <  6:4929/149‑1 (MQ=255)           cTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTTGGGTAgtggtg                                                                                                                                           >  5:245020/1‑147 (MQ=255)                ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCagcagc                                                                                                                                    <  5:96859/149‑1 (MQ=255)                ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCAGCACTACCGcagcag                                                                                                                                     <  5:183497/148‑1 (MQ=255)                 cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcg                                                                                                                                   >  5:139834/1‑149 (MQ=255)                 cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcg                                                                                                                                   >  4:283209/1‑149 (MQ=255)                          ccAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGGGGTTGTGCTGCAGCAGCGCCTGCaact                                                                                                                          >  5:19569/1‑147 (MQ=255)                            aGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCg                                                                                                                        <  2:222168/149‑1 (MQ=255)                            aGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCAGTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGGTGCGGGTGTGCTGGAGCGGCGGCTggacttgg                                                                                                                        >  8:81886/1‑142 (MQ=255)                              gTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGcc                                                                                                                      <  3:283209/149‑1 (MQ=255)                               tACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGccc                                                                                                                     <  4:239522/149‑1 (MQ=255)                                   cccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTg                                                                                                                  >  7:91747/1‑148 (MQ=255)                                    ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTg                                                                                                                               <  3:155767/134‑1 (MQ=255)                                    ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGa                                                                                                                 >  1:187936/1‑148 (MQ=255)                                       gCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGaaaaa                                                                                                             <  2:187936/149‑1 (MQ=255)                                                ccTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGc                                                                                                    >  2:219219/1‑149 (MQ=255)                                                    gcagcaCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATAATCAGGCTTTCGATGCTGTTGGCATAACTGTTGAACACCGGGTCGTCGCTGCTGGGGTATTGCTGCAGACCCGCCTTCAATTCTCCCATGAAAAAACGGTCGGACggt                                                                                                >  3:289070/1‑147 (MQ=255)                                                      agcaCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTg                                                                                              <  8:13830/149‑1 (MQ=255)                                                        caCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGcc                                                                                            <  5:73218/149‑1 (MQ=255)                                                         aCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCg                                                                                           <  8:100836/149‑1 (MQ=255)                                                            tGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        >  1:51829/1‑149 (MQ=255)                                                              ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  6:153133/147‑1 (MQ=255)                                                              ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  6:139834/147‑1 (MQ=255)                                                              ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAAc                                                                                       <  8:91747/148‑1 (MQ=255)                                                              ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAAc                                                                                       <  1:219219/148‑1 (MQ=255)                                                              ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTATCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  5:261584/147‑1 (MQ=255)                                                                     ttttCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAg                                                                               <  2:51829/149‑1 (MQ=255)                                                                        tCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAgcg                                                                             <  6:37080/148‑1 (MQ=255)                                                                         ccAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGAGCCTGCAAGTAGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCAACGAccgc                                                                            >  6:25864/1‑148 (MQ=255)                                                                          cagcagTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCt                                                                          >  2:197782/1‑149 (MQ=255)                                                                            gcagTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCTCTGTTTAACACCGGGTCGTCGCTGGTTCGGTAGTGGAGCAGCAGCGGCCGGAAGTCGCCCTTTAAAAACCGCTAGGCCGTTTTCCGGAACCCAAGCAGCGCCTcg                                                                        >  1:59013/1‑149 (MQ=255)                                                                                        ggATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑Ac                                                           <  1:147393/149‑1 (MQ=255)                                                                                          atcatcAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAg                                                         <  7:81886/149‑1 (MQ=255)                                                                                                    tttCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCCTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTc                                                >  6:310150/1‑148 (MQ=255)                                                                                                        gATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGAccc                                           >  3:238345/1‑149 (MQ=255)                                                                                                         aTGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCg                                          >  2:55560/1‑149 (MQ=255)                                                                                                                 ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCta                                   <  2:59013/148‑1 (MQ=255)                                                                                                                 ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCta                                   <  7:16138/148‑1 (MQ=255)                                                                                                                 ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAACGCTCGGCCGTTTTCCGGAACCCATAGAGCGCCTCGCGCGAATCGA‑ACAAGTGCGCCTCGACCCCCTGGCta                                   >  8:16138/1‑148 (MQ=255)                                                                                                                    aTCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGt                               >  3:184716/1‑149 (MQ=255)                                                                                                                        cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGATACAGGTGCGCCTCGAACCGCTGGCTATAGTCCt                            <  1:197782/149‑1 (MQ=255)                                                                                                                           ttGAACACCGGGTCGTCGCTGCTGCGGTATTGCTGCAGCAGCGACTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        <  8:218232/149‑1 (MQ=255)                                                                                                                           ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  4:25755/148‑1 (MQ=255)                                                                                                                           ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  6:182712/148‑1 (MQ=255)                                                                                                                             gAACACCGGGTCGTCGCTTGTGCGGTAGTGCTGCAGCACCGCCATCAAGTACCCCGTTAAAAAACGCTCGGCCGTTTGCCGGCACCCAACGAGCGCCACGCGCGAATCGA‑AAAGGTGCGCCTCGACCCGCTGGCTATAGTcaatgaac                       >  8:15386/1‑141 (MQ=255)                                                                                                                              aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          >  3:197823/1‑144 (MQ=255)                                                                                                                              aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          <  4:197823/144‑1 (MQ=255)                                                                                                                               acacCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGc                    <  4:184716/149‑1 (MQ=255)                                                                                                                                       tcgtcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc            <  6:245020/149‑1 (MQ=255)                                                                                                                                       caggcgCTGCTGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc            >  1:291587/5‑149 (MQ=255)                                                                                                                                          tcgCTGCTGCGGTAGTGCTGCAGCAGGGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg         <  5:25864/149‑1 (MQ=255)                                                                                                                                          tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          <  3:42223/148‑1 (MQ=255)                                                                                                                                          tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          >  4:42223/1‑148 (MQ=255)                                                                                                                                          tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg         >  3:129517/1‑149 (MQ=255)                                                                                                                                            gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg               <  7:152340/141‑1 (MQ=255)                                                                                                                                            gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg               >  8:152340/1‑141 (MQ=255)                                                                                                                                            gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa       >  7:185555/1‑149 (MQ=255)                                                                                                                                                ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  7:15386/148‑1 (MQ=255)                                                                                                                                                ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  3:23979/148‑1 (MQ=255)                                                                                                                                                   cGGTAGTGCTGCAGCATCGCCTGCAAGTCGCAGTTGAAATAGCGCTCGGCCGTGTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAcc   <  7:206389/146‑1 (MQ=255)                                                                                                                                                   cGGTAGTGCTGCAGCAGCGCCTGCAAGTCTCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCTCCTCGCGCGATTCGA‑ACAGGTTCGCCTCGACCCGCTGTCTATAGTCATTGACCAGCCCCGGGGCCTGGATACCg  <  4:289070/147‑1 (MQ=255)                                                                                                                                                   cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg  <  6:19569/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199877‑5200169

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