Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,475,794	A→G	6.5%	A47A (GCT→GCC)	PP_3965 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,475,794	0	A	G	6.5%	174.3 / 2.4	61	A47A (GCT→GCC)	PP_3965	transposase
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base G (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATT  >  NC_002947/4475646‑4475935                                                                                                                                                     |                                                                                                                                               gccgggcaagtgatcacggtgaaatggccagcctcggatcccgagggttgtgcgcaatttattcgggctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTa                                                                                                                                                >  5:278470‑M1/147‑149 (MQ=255) gccgggcaagtgatcacggtgaaatggccagcctcggatcccgagggttgtgcacaatttattcgggctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTa                                                                                                                                                >  5:293669‑M1/147‑149 (MQ=255)  ccgggcaagtgatcacggtgaaatggccagcctcggatcccgagggttgtgcgcaatttattcgggctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTAg                                                                                                                                               >  7:61412‑M1/146‑149 (MQ=255)    ccAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGcc                                                                                                                                             >  3:214476/1‑149 (MQ=255)     cAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGccc                                                                                                                                            <  4:79191/149‑1 (MQ=255)      aGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCg                                                                                                                                           <  2:262278/149‑1 (MQ=255)       gCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCgg                                                                                                                                          <  7:52414/149‑1 (MQ=255)        cACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGc                                                                                                                                         <  8:228758/149‑1 (MQ=255)           cAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGAt                                                                                                                                      >  1:215872/1‑149 (MQ=255)            aaCGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGAtt                                                                                                                                     <  6:211421/149‑1 (MQ=255)             aCGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCGTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTg                                                                                                                                    <  1:66219/149‑1 (MQ=255)                acggtgaaatggccagcctcggatcccgagggttgtgcgcaatttattcgggctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTGGCCCGGCGATTGGc                                                                                                                                  <  6:217120‑M1/17‑1 (MQ=255)                           tCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAg‑c                                                                                                                     >  7:1221/1‑149 (MQ=255)                             gCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACACCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAg‑cg                                                                                                                    >  5:110833/1‑148 (MQ=255)                                gTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCCCTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGGTAAG‑CGCCGt                                                                                                                >  8:49848/1‑149 (MQ=255)                                gTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGt                                                                                                                >  4:80880/1‑149 (MQ=255)                                 gcGAGTTGGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGt                                                                                                                <  7:296744/147‑1 (MQ=255)                                   gAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCCGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCg                                                                                                             >  5:163877/1‑149 (MQ=255)                                     gTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCAGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGc                                                                                                           >  1:11507/1‑149 (MQ=255)                                      ttCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCt                                                                                                          >  6:135144/1‑149 (MQ=255)                                         gACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGACCAGATAAG‑CGCAGTGCGGCTTc                                                                                                        <  7:49848/148‑1 (MQ=255)                                         gACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTcg                                                                                                       >  8:248250/1‑149 (MQ=255)                                         gACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTcg                                                                                                       >  1:132736/1‑149 (MQ=255)                                          aCCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCCGCACTGCgcttatc                                                                                                       <  8:1221/149‑7 (MQ=17)                                          aCCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTctc                                                                                                      <  4:214476/149‑3 (MQ=255)                                          aCCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTcgc                                                                                                      <  8:27540/149‑1 (MQ=255)                                               gCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCAc                                                                                                 <  5:321759/149‑1 (MQ=255)                                               gCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCAc                                                                                                 >  3:204116/1‑149 (MQ=255)                                                     cgcaatttattcgggctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTGGCCCGGCGATTGGCAAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCAAaca                                                                                           <  4:299271‑M1/55‑1 (MQ=255)                                                        gCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACAcc                                                                                         <  2:215872/148‑1 (MQ=255)                                                        gCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACAcc                                                                                         <  1:294310/148‑1 (MQ=255)                                                        gCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCg                                                                                        <  5:194437/149‑1 (MQ=255)                                                        aatttattcgggctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCg                                                                                        >  1:101264‑M1/92‑149 (MQ=255)                                                                 gggCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGACCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCAccc                                                                               >  3:277790/1‑149 (MQ=255)                                                                    ctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTGGCCCGGCGATTGGCAAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCAAACACCGCAATCACCCggg                                                                            >  3:299271‑M1/80‑149 (MQ=255)                                                                    ctttcgctcaatgatccgcatcgatgccatctggctggccaccgagccgatggacatgcgcgctggcaccgagacagcaTTGGCCCGGCGATTGGCAAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCAAACACCGCAATCACCCggg                                                                            >  5:217120‑M1/80‑149 (MQ=255)                                                                          gAACTTGCCTTGGTTCAGTCTACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATg                                                                      <  3:80880/149‑1 (MQ=255)                                                                          gAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAAt                                                                       <  3:165366/148‑1 (MQ=255)                                                                          gAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACACCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAAt                                                                       <  6:110833/148‑1 (MQ=255)                                                                             cTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTg                                                                   >  7:127521/1‑149 (MQ=255)                                                                                  cTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCg                                                              >  4:105461/1‑149 (MQ=255)                                                                                   ttGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCgg                                                             >  6:64744/1‑149 (MQ=255)                                                                                     ggTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGt                                                            <  5:135144/148‑1 (MQ=255)                                                                                     ggTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGt                                                            <  7:248250/148‑1 (MQ=255)                                                                                                cgcgCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACTGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCAt                                                <  5:134832/149‑1 (MQ=255)                                                                                                   gcTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTc                                             >  2:247834/1‑149 (MQ=255)                                                                                                       ccagccagACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCa                                           <  1:71714/147‑1 (MQ=255)                                                                                                           ccagACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCt                                     >  7:201878/1‑149 (MQ=255)                                                                                                                cacCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTg                                <  3:6036/149‑1 (MQ=255)                                                                                                                 acCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTgg                               >  6:144677/1‑149 (MQ=255)                                                                                                                 acCGACGCCATCGTGCACCAGCACTTTCATACGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCCCCCGGGCAAATGCGGTCTCCGTGCCACCGCGCATGTCCATCGGCTCGGTgg                               >  3:257594/1‑149 (MQ=14)                                                                                                                    gACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGgcca                            <  5:144677/149‑1 (MQ=255)                                                                                                                         cATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGgccagc                          >  8:183945/1‑146 (MQ=255)                                                                                                                         cATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGgccagc                          >  1:225813/1‑146 (MQ=255)                                                                                                                         cATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGccagccag                       <  3:105461/149‑1 (MQ=255)                                                                                                                                cACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCAt                <  6:163877/149‑1 (MQ=255)                                                                                                                                cACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCAt                <  2:11507/149‑1 (MQ=255)                                                                                                                                        cTTTCATCCGGTTAGACCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCgg        <  4:277790/149‑1 (MQ=255)                                                                                                                                         tttCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGt                                                            >  1:54872/1‑96 (MQ=255)                                                                                                                                         tttCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGt                                                            <  2:54872/96‑1 (MQ=255)                                                                                                                                         tttCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGa       >  5:10206/1‑149 (MQ=39)                                                                                                                                           tCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGa       <  2:225813/147‑1 (MQ=39)                                                                                                                                              tCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCAtt  <  6:10206/149‑1 (MQ=38)                                                                                                                                                     |                                                                                                                                               GGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAG‑CGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATT  >  NC_002947/4475646‑4475935

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