Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,375,446	T→G	11.3%	R30R (CGA→CGC)	PP_5637 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,375,446	0	T	G	11.3%	160.7 / 13.8	62	R30R (CGA→CGC)	PP_5637	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base G (5/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

ACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGT  >  NC_002947/4375298‑4375587                                                                                                                                                       |                                                                                                                                                  acGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACt                                                                                                                                                   >  7:89342/1‑149 (MQ=255) acGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACt                                                                                                                                                   <  7:195805/149‑1 (MQ=255)   gACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCg                                                                                                                                                 >  7:21158/1‑149 (MQ=255)    aCTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑a                                                                                                                                            <  4:176060/149‑1 (MQ=255)     cTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑Ac                                                                                                                                           >  7:131186/1‑149 (MQ=255)        cgcATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑Acc                                                                                                                                          <  8:223129/147‑1 (MQ=255)           atatAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCagag                                                                                                                                     <  5:161113/149‑1 (MQ=255)               aGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTc                                                                                                                                 >  1:186410/1‑149 (MQ=255)                gCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCg                                                                                                                                <  8:21158/149‑1 (MQ=255)                    gCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCAAGGCCTGGGTCGCGACTc                                                                                                                                                  >  5:30836/1‑131 (MQ=255)                      cTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGa                                                                                                                          <  5:2587/149‑1 (MQ=255)                       ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTATGACCTTCGCATCGTGCAGCCGGGGGATAATCTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGa                                                                                                                          >  1:276023/1‑149 (MQ=255)                       ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                                 <  7:299302/141‑7 (MQ=255)                       ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                                 >  8:299302/1‑135 (MQ=255)                          gACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                                 >  7:315873/1‑132 (MQ=255)                          gACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                                 <  8:315873/138‑7 (MQ=255)                          gACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑CGCCAGAGGCTCGCTGGGATGGt                                                                                                                      >  4:109456/1‑149 (MQ=255)                           aCCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTc                                                                                                                     <  8:89342/149‑1 (MQ=255)                           aCCTCTGCTACCCCGTTGATGTTGATGCCGAGCGGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATACTGTCTGG‑CCCAGGCGTAGGTTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTc                                                                                                                     <  4:107921/149‑1 (MQ=255)                                  ctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                               <  6:170923/148‑1 (MQ=255)                                     ccccTTTGATTGTGATGCGTTGCGGGGCCTTGCGGT‑TTACCTTCTCATCGTGCAGCCGGGGGATAATGTGGGG‑CCCAGGCGTTGGGTTAAACACCCAGGCCTTGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                            <  8:131186/148‑1 (MQ=255)                                       ccGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGAc                                                                                                         >  7:302769/1‑149 (MQ=255)                                             atgttgatgCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgc                                                                                                   >  4:135923/1‑149 (MQ=255)                                                ttgatgCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgcg                                                                                                <  1:148182/149‑1 (MQ=255)                                                   atgCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAg                                                                                             <  4:43746/149‑1 (MQ=255)                                                     gCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                 >  1:278776/1‑143 (MQ=255)                                                     gCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                 <  2:278776/143‑1 (MQ=255)                                                     gCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAggg                                                                                           >  7:212732/1‑149 (MQ=255)                                                      ccGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTTGGGCCCCCGGCGTAGGGTTAAACACCCAGGCATGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAggg                                                                                           >  1:96417/1‑149 (MQ=255)                                                      ccGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGTGGTGGAGAGAGGCTCGCTGTGATGGTCGACAATGTGGACGTCcgcgcg                                                                                                >  1:308123/1‑147 (MQ=255)                                                            tGGGCCTTGCCGT‑TGACATTGGCATGGTGCTGCGGGGGTATAATGTCGGG‑CCCAGGCTTTGGTTAAAACACCTAGGCGTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGAc                                                                                    <  2:186410/149‑1 (MQ=255)                                                             gggCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACTCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACt                                                                                   >  8:47563/1‑149 (MQ=255)                                                               gCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACttt                                                                                 >  1:28975/1‑149 (MQ=255)                                                                ccATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTc                                                                                <  1:340647/149‑1 (MQ=255)                                                                     ccGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACTCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                            <  7:47563/148‑1 (MQ=255)                                                                      cGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAgg                                                                          <  8:302769/149‑1 (MQ=255)                                                                            aCCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACttt                                                                                 <  5:90962/137‑1 (MQ=255)                                                                            aCCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACttt                                                                                 >  6:90962/1‑137 (MQ=255)                                                                                 cGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                 <  8:212732/148‑1 (MQ=255)                                                                                  gCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAgg                                                                <  2:96417/148‑1 (MQ=255)                                                                                    aTCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGc                                                              >  6:6707/1‑148 (MQ=255)                                                                                          cAGCCGGGGGATAATGTCGGG‑CCCAGGCGTTGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCa                                                       >  1:30712/1‑149 (MQ=255)                                                                                                        tGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCgcgg                                         >  2:34037/1‑149 (MQ=255)                                                                                                           cGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGt                                      <  7:94763/149‑1 (MQ=255)                                                                                                              g‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGa                                   <  6:320910/149‑1 (MQ=255)                                                                                                                 ccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTTGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAAc                                 <  5:178601/149‑1 (MQ=255)                                                                                                                   aGGCGTAGTGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACga                               <  2:308123/149‑1 (MQ=255)                                                                                                                       gTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaag                           <  5:6707/149‑1 (MQ=255)                                                                                                                          gggTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaagg                          >  4:24445/1‑147 (MQ=255)                                                                                                                           tgTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                       <  2:276023/148‑1 (MQ=255)                                                                                                                           ggTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                       >  6:13023/1‑149 (MQ=255)                                                                                                                           ggTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACCTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                       >  4:231533/1‑149 (MQ=255)                                                                                                                            gTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                              <  6:201308/141‑1 (MQ=255)                                                                                                                            gTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                              >  5:201308/1‑141 (MQ=255)                                                                                                                                        cctcTGGGTCGAGTCGCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCAc           >  7:206795/4‑148 (MQ=255)                                                                                                                                        cctcTGGGTCGAGTCGCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCAc           >  7:201491/4‑148 (MQ=255)                                                                                                                                          tcTGGGTCGAGTCGCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCaa        >  3:57646/2‑149 (MQ=255)                                                                                                                                          ccTGGGTCGCGACTCG‑‑‑‑ACCGAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCaa        <  3:109456/149‑1 (MQ=255)                                                                                                                                             gggTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg              <  6:196842/140‑1 (MQ=255)                                                                                                                                             gggTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg              >  5:196842/1‑140 (MQ=255)                                                                                                                                               gTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTg   >  8:46685/1‑149 (MQ=255)                                                                                                                                                tCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGt  <  1:359246/149‑1 (MQ=255)                                                                                                                                                   cgACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGcc                                                   <  7:227581/97‑1 (MQ=255)                                                                                                                                                   cgACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGcc                                                   >  8:227581/1‑97 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGT  >  NC_002947/4375298‑4375587

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