Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,432,875	C→A	6.7%	L319M (CTG→ATG)	xenA →	protein XenA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,432,875	0	C	A	6.7%	219.3 / 4.3	69	L319M (CTG→ATG)	xenA	protein XenA
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base A (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

ACCATTCCCGACACCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGAGGCTTGGAGGG  >  NC_002947/1432728‑1433023                                                                                                                                                    |                                                                                                                                                      aCCATTCCCGACACCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGAGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCt                                                                                                                                                      >  1:118720/1‑149 (MQ=255)        ccGACACCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTg                                                                                                                                               <  6:300051/149‑1 (MQ=255)        ccGACACCAACATCCCCTGGAGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCtt                                                                                                                                               <  5:354577/149‑2 (MQ=255)         cGACACCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTgg                                                                                                                                              <  3:224916/149‑1 (MQ=255)         cGACACCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTgg                                                                                                                                              >  1:319396/1‑149 (MQ=255)             acCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTc                                                                                                                                          >  2:197206/1‑149 (MQ=255)              ccAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑g                                                                                                                                        <  1:125221/149‑1 (MQ=255)              ccAACATACCGTGGGGTGCTGCGTTCATGGGGCCGTTTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGTCTTTGGTAGGCCGAAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑g                                                                                                                                        <  8:71175/149‑1 (MQ=255)               cAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTGGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGATCCAGCTGGttggcc                                                                                                                                        <  5:105429/149‑5 (MQ=255)               cAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTc                                                                                                                                          >  7:309173/1‑147 (MQ=255)               cAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTc                                                                                                                                          <  8:309173/147‑1 (MQ=255)               cAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑gg                                                                                                                                       >  7:52709/1‑149 (MQ=255)                aaCATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑gg                                                                                                                                       >  2:341598/1‑148 (MQ=255)                    tCCCCTGGGGTCCAGCGTTCATAGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAggg                                                                                                                                  >  6:178046/1‑149 (MQ=255)                              tCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGGGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCAcc                                                                                                                        <  3:147528/149‑1 (MQ=255)                                cAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTg                                                                                                                      >  2:260833/1‑149 (MQ=255)                                 aGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGGGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGGAGGGCGCGCGCACCTgg                                                                                                                     <  7:169187/149‑1 (MQ=255)                                     ttCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGa                                                                                                                 >  5:42296/1‑149 (MQ=255)                                       cATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGAc                                                                                                                <  1:341598/148‑1 (MQ=255)                                        aTGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGAccc                                                                                                              >  5:217355/1‑149 (MQ=255)                                        aTGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGAccc                                                                                                              >  8:164385/1‑149 (MQ=255)                                         tGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCg                                                                                                             <  7:164385/149‑1 (MQ=255)                                             gCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGGTTC‑GGTAGGGCGCGCGCACCTGGCCGCCCCGCACt                                                                                                         >  4:285796/1‑149 (MQ=255)                                                   tGCCGAGCGCGTGCGCCGGGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCt                                                                                                    >  3:287213/1‑148 (MQ=255)                                                     ccGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACTTGGCCGACCCGCACTGGGCTTa                                                                                                  <  1:197206/148‑1 (MQ=255)                                                     ccGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGATCCAGCTGGttggc                                                                                                                                         >  2:23127/1‑107 (MQ=255)                                                     ccGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGATCCAGCTGGttggc                                                                                                                                         <  1:23127/110‑4 (MQ=255)                                                       gAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACtt                                                                                               <  5:178046/149‑1 (MQ=255)                                                         gcgcGTGCGCCGTGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTtg                                                                                             <  4:287213/149‑1 (MQ=255)                                                          cgcgTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTtgc                                                                                            <  1:260833/149‑1 (MQ=255)                                                            cgTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTtgctg                                                                                          >  5:340601/1‑149 (MQ=255)                                                             gTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGATCCAGCTGGttggc                                                                                                                                         >  2:99324/1‑99 (MQ=255)                                                             gTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGATCCAGCTGGttggc                                                                                                                                         <  1:99324/102‑4 (MQ=255)                                                              tGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGcc                                                                                        <  7:24660/149‑1 (MQ=255)                                                                    gcgaagcgaaGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTtgct                                                                                           >  2:216908/1‑140 (MQ=255)                                                                    gcgaagcgaaGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTtgct                                                                                           <  1:216908/140‑1 (MQ=255)                                                                       aagcgaagcTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGACTTACTTTGCTGCCAAGGAGCTg                                                                               <  1:122703/149‑1 (MQ=255)                                                                        agcgaagcTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTgg                                                                              <  5:144218/149‑1 (MQ=255)                                                                          cgaagcTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTgggg                                                                            >  7:214055/1‑149 (MQ=255)                                                                          cgaagcTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGAGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTgggg                                                                            >  7:219727/1‑149 (MQ=255)                                                                                   ccGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTg                                                                                                        <  1:202026/112‑1 (MQ=255)                                                                                   ccGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTg                                                                                                        >  2:202026/1‑112 (MQ=255)                                                                                   ccGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCaa                                                                                      >  4:292418/1‑130 (MQ=255)                                                                                   ccGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCaa                                                                                      <  3:292418/130‑1 (MQ=255)                                                                                   ccGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGaaaaa                                                                   <  5:180070/149‑1 (MQ=255)                                                                                     gTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGc                                                                 <  6:42296/149‑1 (MQ=255)                                                                                          gTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGcctcct                                                            >  4:357915/1‑149 (MQ=255)                                                                                                 tGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCtt                                                     >  2:93775/1‑149 (MQ=255)                                                                                                  ggggCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTg                                                    >  4:165039/1‑149 (MQ=255)                                                                                                          ggTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGgcgc                                             <  2:319396/148‑1 (MQ=255)                                                                                                            tACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCt                                           >  5:29398/1‑148 (MQ=255)                                                                                                               gccgcAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCGTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATg                                       >  7:172035/1‑149 (MQ=255)                                                                                                                       tGGCTGAGGACGCATTGCAGGCCAACCAGCTGGATCTGGGTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTgg                               <  8:172035/149‑1 (MQ=255)                                                                                                                         gCGGAGGCCGCATTGCATGCCAACCAGCTTGATCTGGTTTC‑GTTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGc                              <  8:52709/148‑1 (MQ=255)                                                                                                                            gAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGc                              <  7:189143/145‑1 (MQ=255)                                                                                                                            gAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGc                              >  8:189143/1‑145 (MQ=255)                                                                                                                            gAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGa                          <  6:340601/149‑1 (MQ=255)                                                                                                                              ggCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGc                        <  6:29398/149‑1 (MQ=255)                                                                                                                               gCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCg                       >  1:49266/1‑149 (MQ=255)                                                                                                                                 cGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGt                      >  5:171512/1‑148 (MQ=255)                                                                                                                                  gCATTGCAGGCCAACCAGCTGGATCTGGTTTCGGGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGG‑AAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTa                    >  1:231424/1‑149 (MQ=255)                                                                                                                                     ttGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCAGTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGa                          <  1:143599/140‑1 (MQ=255)                                                                                                                                     ttGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGa                          >  2:143599/1‑140 (MQ=255)                                                                                                                                        cAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCt               >  6:210633/1‑148 (MQ=255)                                                                                                                                         aGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGa             <  3:285796/149‑1 (MQ=255)                                                                                                                                             cAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCAAACCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGAGGCt         >  1:123337/1‑149 (MQ=255)                                                                                                                                               aCCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGGGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGAGGCTTg       <  6:217355/149‑1 (MQ=255)                                                                                                                                                   gCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGAGGCTTGGAgg   <  2:118720/149‑1 (MQ=255)                                                                                                                                                   gCTGGATCTGGTTTC‑GGTAGGGCCCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCCCTGATGCTTGCAgg   >  7:257965/1‑149 (MQ=255)                                                                                                                                                    cTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGAGGCTTGGAggg  <  3:357915/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ACCATTCCCGACACCAACATCCCCTGGGGTCCAGCGTTCATGGGGCCGATTGCCGAGCGCGTGCGCCGCGAAGCGAAGCTGCCCGTGACGTCGGCGTGGGGCTTTGGTACGCCGCAGTTGGCGGAGGCCGCATTGCAGGCCAACCAGCTGGATCTGGTTTC‑GGTAGGGCGCGCGCACCTGGCCGACCCGCACTGGGCTTACTTTGCTGCCAAGGAGCTGGGGGTGGAAAAAGCCTCCTGGACCTTGCCGGCGCCTTATGCACACTGGCTCGAGCGTTATCGCTGAGGCTTGGAGGG  >  NC_002947/1432728‑1433023

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