Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	495,077	T→G	6.2%	Q319H (CAA→CAC) ‡	PP_0408 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	495,077	0	T	G	6.2%	197.8 / 3.7	64	Q319H (CAA→CAC) ‡	PP_0408	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base G (2/2);  total (33/31)
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.49e-01

GGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTGG  >  NC_002947/494932‑495225                                                                                                                                                  |                                                                                                                                                     ggTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGc                                                                                                                                                   <  3:37429/149‑1 (MQ=255)  gTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGct                                                                                                                                                  <  5:201803/149‑1 (MQ=255)         aaGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTc                                                                                                                                           >  8:157899/1‑149 (MQ=255)         aaGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAAGAGCAATTGTTc                                                                                                                                           >  3:147151/1‑149 (MQ=255)                cgcgcgTACGCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAAGAGCAATTGTTccgccggg                                                                                                                                    <  4:147151/149‑3 (MQ=255)                           gTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtgg                                                                                                                         >  3:147122/1‑149 (MQ=255)                              aGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtggcat                                                                                                                      <  8:158716/149‑1 (MQ=255)                               gTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtggcatg                                                                                                                     >  4:124119/1‑149 (MQ=255)                               gTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtg                                                                                                                          <  1:93865/144‑1 (MQ=255)                               gTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtg                                                                                                                          >  2:93865/1‑144 (MQ=255)                                    aCCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGcaat                                                                                                                <  6:233811/149‑1 (MQ=255)                                           cTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGAc                                                                                                         <  7:157899/149‑1 (MQ=255)                                            tATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGctc                                                                                                                                                 <  8:294443/108‑1 (MQ=255)                                            tATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGctc                                                                                                                                                 >  7:294443/1‑108 (MQ=255)                                            tATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGAc                                                                                                         >  1:193154/1‑148 (MQ=255)                                             aTGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACaa                                                                                                       >  8:76151/1‑149 (MQ=255)                                                   tGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCg                                                                                                 >  4:220834/1‑149 (MQ=255)                                                     aTCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCgg                                                                                                <  8:191578/148‑1 (MQ=255)                                                      tCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtggcatgg                                                                                                                    >  5:374712/1‑127 (MQ=255)                                                      tCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtggcatgg                                                                                                                    <  6:374712/127‑1 (MQ=255)                                                         gCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTAATg                                                                                           >  5:262626/1‑149 (MQ=255)                                                          cAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTgg                                                                                          <  5:47869/149‑1 (MQ=255)                                                           aaGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGc                                                                                         >  2:250911/1‑149 (MQ=255)                                                             gcggcgcggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTgg                                                                                          <  8:313316/146‑1 (MQ=255)                                                             gcggcgcggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGcag                                                                                       >  7:313316/1‑147 (MQ=255)                                                                  gcggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCt                                                                                   >  8:14712/1‑148 (MQ=255)                                                                    ggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACt                                                                                            <  3:114657/137‑1 (MQ=255)                                                                    ggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACt                                                                                            >  4:114657/1‑137 (MQ=255)                                                                        tttGACGCCCCGCCACAGGGATCAACCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCATGTTCCGCCGCTATGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACAGGCAATCCATCTGCCg                                                                            >  8:155654/1‑149 (MQ=255)                                                                            aCGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTg                                                                        >  4:260010/1‑149 (MQ=255)                                                                             cGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGc                                                                       <  7:14712/149‑1 (MQ=255)                                                                                   gCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCt                                                                 <  4:215399/149‑1 (MQ=255)                                                                                       caGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATAAATGACAATTTTCGGTACTGGCAATCATTCTGCCGTTTGCTGGCATGCAg                                                             >  6:26707/1‑149 (MQ=255)                                                                                        aGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGTCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGc                                                            >  6:287683/1‑149 (MQ=255)                                                                                        aGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGc                                                            >  4:246656/1‑149 (MQ=255)                                                                                        aGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCA‑CCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGcc                                                           >  3:298459/1‑149 (MQ=255)                                                                                        aGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCA‑CCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGcc                                                           >  3:293436/1‑149 (MQ=255)                                                                                                aGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTg                                                    <  2:193154/149‑1 (MQ=255)                                                                                                 gCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTgg                                                                                          >  6:379592/1‑110 (MQ=255)                                                                                                 gCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTgg                                                                                          <  5:379592/110‑1 (MQ=255)                                                                                                         cAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCCGCAGCCGTTGCTGGGCTGCGtt                                           >  4:131507/1‑149 (MQ=255)                                                                                                          aGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTg                                          <  2:255786/149‑1 (MQ=255)                                                                                                            gccTGAGGGCTCATCCACCCCCTCACCCTGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCg                                        <  5:26707/148‑1 (MQ=255)                                                                                                            cccTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCg                                        <  3:260010/149‑1 (MQ=255)                                                                                                            cccTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTAATGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCg                                        <  6:262626/149‑1 (MQ=255)                                                                                                             ccTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTCCAGCCGTTGCTGGGCTGCGTTGCgg                                       >  4:281095/1‑149 (MQ=255)                                                                                                              cTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGt                                      >  8:92642/1‑149 (MQ=255)                                                                                                                      tCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTg                              >  1:94768/1‑149 (MQ=255)                                                                                                                        aTCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGc                             <  1:40649/148‑1 (MQ=255)                                                                                                                        aTCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCg                            >  5:108343/1‑149 (MQ=255)                                                                                                                        aTCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCg                            >  2:171873/1‑149 (MQ=255)                                                                                                                        aTCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGAAATCAATGACAATTGTAGGTACTGGCTATCCTTCTTCCGTTTGCTGGCCTGCAGCCTTTTATGGGCTGCGTTGCGGTCAGAGCTGCg                            >  2:177535/1‑149 (MQ=255)                                                                                                                           cACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTg                                                                     >  4:197469/1‑105 (MQ=255)                                                                                                                           cACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTg                                                                     <  3:197469/105‑1 (MQ=255)                                                                                                                            aCCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCCGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTg                        <  3:131507/149‑1 (MQ=255)                                                                                                                                  gCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtag                  >  4:215259/1‑149 (MQ=255)                                                                                                                                   cACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtagt                 <  5:287683/149‑1 (MQ=255)                                                                                                                                   cACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtagt                 <  3:246656/149‑1 (MQ=255)                                                                                                                                    aCCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtagta                <  3:215259/149‑1 (MQ=255)                                                                                                                                    aCCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtagta                <  7:92642/149‑1 (MQ=255)                                                                                                                                        ggcggAACAAGAGCAATTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGata            >  3:314927/1‑149 (MQ=255)                                                                                                                                          cggAACAAGAGCAATTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTg                                                                     <  6:197416/90‑1 (MQ=39)                                                                                                                                                aaTTGCTCTTGTTCCTCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGTCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGAGGGTGGCGTAGTAGATATCGGCCGt    <  1:177535/149‑1 (MQ=255)                                                                                                                                                aaTTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGt    <  1:171873/149‑1 (MQ=255)                                                                                                                                                  ttGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTgg  <  4:147122/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTGG  >  NC_002947/494932‑495225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑