Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,495,653	C→A	5.8%	L28M (CTG→ATG)	mdeA →	methionine gamma‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,495,653	0	C	A	5.8%	246.4 / 3.8	78	L28M (CTG→ATG)	mdeA	methionine gamma‑lyase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base A (2/3);  total (38/40)
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.99e-01

GCTATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCT  >  NC_002947/1495505‑1495796                                                                                                                                                      |                                                                                                                                                     gCTATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGccc                                                                                                                                                      >  4:109194/1‑149 (MQ=255)     ttttGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTg                                                                                                                                                    <  3:371778/147‑1 (MQ=255)       ttGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTg                                                                                                                                                 <  3:419587/148‑1 (MQ=255)           cgatcttCGGGGAGCTCGGTGGCATAAGCTTATAGCATACCACTAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCAGGGCGCCAcc                                                                                                                                            <  3:312433/143‑1 (MQ=255)            aaaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCg                                                                                                                                           >  5:200102/1‑149 (MQ=255)            aaaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATAATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCg                                                                                                                                           >  5:355565/1‑149 (MQ=255)             taaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCAGGGCGCCAccat                                                                                                                                          <  3:411387/148‑3 (MQ=255)                aTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgtgt                                                                                                                                       <  3:232148/149‑1 (MQ=255)                aTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCAGGGCGCCAcc                                                                                                                                            >  3:36478/1‑144 (MQ=255)                aTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCAGGGCGCCAcc                                                                                                                                            <  4:36478/144‑1 (MQ=255)                 tCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGTCGCCCTGGTGCCACCGGTGTa                                                                                                                                      <  7:156541/149‑1 (MQ=255)                     gCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGt                                                                                                                                                  >  1:284240/1‑133 (MQ=255)                     gCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGt                                                                                                                                                  <  2:284240/133‑1 (MQ=255)                      cAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑a                                                                                                                              <  6:52307/149‑1 (MQ=255)                         cTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACg                                                                                                                            <  1:126167/149‑1 (MQ=255)                         cTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACgg                                                                                                                           <  3:109194/149‑1 (MQ=255)                            ggTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGc                                                                                                                          <  5:115388/147‑1 (MQ=255)                             gTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTAc                                                                                                                       <  1:122306/149‑1 (MQ=255)                                 cATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑At                                                                                                                  <  8:392411/149‑1 (MQ=255)                                  aTAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGTAGGCGCCCTGGTGCCACCGGTGTAACAG‑‑‑ACGGCTACCT‑aag                                                                                                                 >  4:298150/1‑147 (MQ=255)                                        ttATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTcc                                                                                                          >  2:354001/1‑149 (MQ=255)                                          ataGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTcc                                                                                                          <  1:241252/147‑1 (MQ=255)                                             gCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACt                                                                                                     >  4:368827/1‑149 (MQ=255)                                                aaCCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACTGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTc                                                                                                  <  5:2970/149‑1 (MQ=255)                                                aaCCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTc                                                                                                  >  1:285410/1‑149 (MQ=255)                                                       agagGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACg                                                                                           <  1:378919/149‑1 (MQ=255)                                                       agagGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACg                                                                                           >  4:351571/1‑149 (MQ=255)                                                       agagGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACg                                                                                           >  4:345768/1‑149 (MQ=255)                                                       agagGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACg                                                                                           <  3:368827/149‑1 (MQ=255)                                                             tGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgt                                                                                                                                         >  4:408375/1‑102 (MQ=255)                                                             tGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgt                                                                                                                                         <  3:408375/102‑1 (MQ=255)                                                               ttGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGtgctgc                                                                                    <  2:161314/148‑1 (MQ=255)                                                                 gCCATGCGCGACTCCCACCACAACACTGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGctgcctg                                                                                 >  1:232443/1‑149 (MQ=255)                                                                   cATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCt                                                                               >  6:389030/1‑149 (MQ=255)                                                                     tGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTc                                                                             >  7:353449/1‑149 (MQ=255)                                                                     tGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTc                                                                             >  2:45926/1‑149 (MQ=255)                                                                     tGCGCGACTCCCACCACAACAACGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTTTGACCCGCTTTCCCATGGTGGCGCCATGGTGCCACCGGTGTACCAG‑‑‑CCTGCTACCT‑AT‑GCCTTCCCCAATGTCGAGTACGGTGCTGCCTGCtgc                                                                             >  8:89527/1‑147 (MQ=255)                                                                        gcgACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGcc                                                                          <  6:200102/149‑1 (MQ=255)                                                                             tCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCggg                                                                       >  6:15147/1‑147 (MQ=255)                                                                              cccACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCACTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑CCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGTggat                                                                    >  8:389086/1‑148 (MQ=255)                                                                               ccaccaCAACCCCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTGCCATGGTGGCGCCCTGGTGCCACGGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCGCACTGTCGAGTACGGTGCTGCCTTCTTCGCCGGgga                                                                     <  2:354702/147‑1 (MQ=255)                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGCCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggag                                                                    >  3:405135/1‑148 (MQ=255)                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTa                                                                                                                        >  8:377326/1‑98 (MQ=255)                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTa                                                                                                                        <  7:377326/98‑1 (MQ=255)                                                                               ccaccaCAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggag                                                                    >  6:187938/1‑148 (MQ=255)                                                                                caccacAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggag                                                                    >  4:295634/1‑147 (MQ=255)                                                                                   cacaACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTc                                                                                                  >  1:157300/1‑114 (MQ=255)                                                                                   cacaACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTc                                                                                                  <  2:157300/114‑1 (MQ=255)                                                                                   cacaACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGcc                                                               >  4:193014/1‑149 (MQ=255)                                                                                    acaacaCCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGGGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCt                                                              <  2:307065/149‑1 (MQ=255)                                                                                         acCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTAGGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCa                                                         >  4:31794/1‑149 (MQ=255)                                                                                         acCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCa                                                         >  6:336079/1‑149 (MQ=255)                                                                                             gTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTCCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGc                                                          <  5:370858/144‑1 (MQ=255)                                                                                             gTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTCCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGc                                                          >  6:370858/1‑144 (MQ=255)                                                                                                    cACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGcc                                              <  7:137408/149‑1 (MQ=255)                                                                                                      cgcgGGCAA‑TTAATCATGGCTATGACCCGCTTTCCCATTGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGATTAAGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGc                                            <  2:149845/149‑1 (MQ=255)                                                                                                        cgGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCAt                                          <  2:285410/149‑1 (MQ=255)                                                                                                                ttCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAc                                   >  7:348870/1‑149 (MQ=255)                                                                                                                  catcatGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAcc                                  <  3:345768/148‑1 (MQ=255)                                                                                                                  catcatGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAcc                                  <  3:351571/148‑1 (MQ=255)                                                                                                                   atcatGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAccc                                 <  3:177322/148‑1 (MQ=255)                                                                                                                   atcatGGATTGGACCCGCTTCCCCCTGGTGGCGCTGTGGTGCCACGGGTGTACCAGCCCACGGCTACATCCTCCCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTcc                                      <  1:354001/148‑1 (MQ=255)                                                                                                                    tcatGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCa                               >  8:22464/1‑149 (MQ=255)                                                                                                                     catGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCAc                              >  5:418196/1‑149 (MQ=255)                                                                                                                       tGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGc                            >  6:432843/1‑149 (MQ=255)                                                                                                                             tGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGc                        <  5:432843/147‑1 (MQ=255)                                                                                                                             tGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGc                        <  2:232443/147‑1 (MQ=255)                                                                                                                              gACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCtt                     >  1:80811/1‑149 (MQ=255)                                                                                                                               aCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTg                    >  6:266872/1‑149 (MQ=255)                                                                                                                                 ccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCt                  >  2:3102/1‑149 (MQ=255)                                                                                                                                    cTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                 >  4:186177/1‑147 (MQ=255)                                                                                                                                       tCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGag              >  2:282472/1‑147 (MQ=255)                                                                                                                                         ccATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagcagc          <  5:266872/149‑1 (MQ=255)                                                                                                                                           aTGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAgcg         >  7:380500/1‑148 (MQ=255)                                                                                                                                              gtggCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCAt      <  5:336079/148‑1 (MQ=255)                                                                                                                                               tggCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATgg    <  8:348870/149‑1 (MQ=255)                                                                                                                                               tggCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATg     <  3:158894/148‑1 (MQ=255)                                                                                                                                               tggCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATg     >  4:158894/1‑148 (MQ=255)                                                                                                                                                 gcgcCCTGGTGCGAACGGTGTACCAG‑‑‑ACGGCTACGT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGACTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCAGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCt  <  3:366755/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCTATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCA‑TTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAG‑‑‑ACGGCTACCT‑AT‑GCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCT  >  NC_002947/1495505‑1495796

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