Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,375,444	A→T	10.4%	V31D (GTC→GAC)	PP_5637 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,375,444	0	A	T	10.4%	172.0 / 9.2	65	V31D (GTC→GAC)	PP_5637	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base T (5/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGT  >  NC_002947/4375296‑4375587                                                                                                                                                       |                                                                                                                                                    ccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGa                                                                                                                                                     >  2:359246/1‑149 (MQ=255)  cacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGAc                                                                                                                                                    <  8:217233/149‑1 (MQ=255)   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                              >  5:201308/1‑141 (MQ=255)                                                                                                                              gTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                              <  6:201308/141‑1 (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                                                   >  8:227581/1‑97 (MQ=255)                                                                                                                                                     cgACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGcc                                                   <  7:227581/97‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    CCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGT  >  NC_002947/4375296‑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: ‑