Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,794,692	A→T	5.0%	L602Q (CTG→CAG)	fcs ←	feruloyl‑CoA‑synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,794,692	0	A	T	5.0%	244.1 / 3.2	81	L602Q (CTG→CAG)	fcs	feruloyl‑CoA‑synthetase
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base T (2/2);  total (40/41)
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.79e-01

TCGACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACT  >  NC_002947/3794543‑3794840                                                                                                                                                      |                                                                                                                                                      tCGACCCAGGGCGTGCGAACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgca                                                                                                                                                       >  1:913115/1‑150 (MQ=255)     cccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑c                                                                                                                                                  <  2:293481/150‑1 (MQ=255)     cccAGGGCGTGCGCACCATCTGCACGATGGCCGTGTCCTCGAAGACGATGCGCAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGAGGTCTGCGCCACGGTACAGCGCCTTGACCTGAGGGGAGCGCCGATgaagaaa                                                                                                                                                   >  1:713298/1‑147 (MQ=255)      ccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAATGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cg                                                                                                                                                 >  2:1102448/1‑150 (MQ=255)      ccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cg                                                                                                                                                 >  2:1867416/1‑150 (MQ=255)      ccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cg                                                                                                                                                 <  2:1651202/150‑1 (MQ=255)      ccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCAAGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cg                                                                                                                                                 >  1:1451782/1‑150 (MQ=255)        aGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGgc                                                                                                                                               <  1:517683/150‑1 (MQ=255)              tGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCAGCGCGCCG‑TGCTGCAGCGgcgcgccg                                                                                                                                         <  2:858708/150‑3 (MQ=255)               gcgcACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCtga                                                                                                                                        <  1:789481/150‑1 (MQ=255)                  cACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCtgattg                                                                                                                                     >  2:893862/1‑150 (MQ=255)                   aCCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCtgattga                                                                                                                                    <  2:24766/150‑1 (MQ=255)                    ccATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCtgattgat                                                                                                                                   <  2:913115/150‑1 (MQ=255)                    ccATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCtgattgat                                                                                                                                   <  2:1185398/150‑1 (MQ=255)                     cATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCCCGCCGCTGCAGCA‑CGGCGCGCTGATTGATc                                                                                                                                  >  1:200866/1‑150 (MQ=255)                       tCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACATCCGGGCCGCACTCAAGGCCGCACCTTGGGATAAAGTCCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGa                                                                                                                                <  1:1773235/150‑1 (MQ=255)                         tGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGc                                                                                                                              >  1:92923/1‑150 (MQ=255)                          ggACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGcc                                                                                                                             >  2:1792666/1‑150 (MQ=255)                           gACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGccc                                                                                                                            >  2:1918407/1‑150 (MQ=255)                             cGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgcag                                                                                                                                                      >  2:1446300/1‑123 (MQ=255)                             cGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgcag                                                                                                                                                      <  1:1446300/123‑1 (MQ=255)                               aTGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGt                                                                                                                        >  1:1310539/1‑150 (MQ=255)                                tGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTc                                                                                                                       >  1:1036238/1‑150 (MQ=255)                                   ccGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGAGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGt                                                                                                                    <  2:1029902/150‑1 (MQ=255)                                          cTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCa                                                                                                             >  1:1450489/1‑150 (MQ=255)                                           tCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCa                                                                                                             <  1:1741144/149‑1 (MQ=255)                                           tCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCa                                                                                                             >  2:1741144/1‑149 (MQ=255)                                           tCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCAc                                                                                                            >  2:929989/1‑150 (MQ=255)                                            cGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgca                                                                                                                                                       <  2:1862997/107‑1 (MQ=255)                                            cGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgca                                                                                                                                                       >  1:1862997/1‑107 (MQ=255)                                             gAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCg                                                                                                          <  1:1128680/150‑1 (MQ=255)                                              aaGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCgg                                                                                                         >  2:891812/1‑150 (MQ=255)                                                  cGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtc                                                                                                     <  1:83251/150‑1 (MQ=255)                                                   ggTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcg                                                                                                    >  2:811679/1‑150 (MQ=255)                                                   ggTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcg                                                                                                    >  2:1673337/1‑150 (MQ=255)                                                      gCACCGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatc                                                                                                 <  1:195852/150‑1 (MQ=255)                                                        acaGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatcga                                                                                               >  1:263357/1‑150 (MQ=255)                                                         caGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAc                                                                                              <  2:1310539/150‑1 (MQ=255)                                                            cccGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAcggc                                                                                           <  1:752258/150‑1 (MQ=255)                                                            cccGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAcggc                                                                                           >  2:901676/1‑150 (MQ=255)                                                             ccGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCC‑CGGCGCGCTGATTGATCGAGCCCGTGTCGGGGATCTCACCGGCGTCGATCGAcggcg                                                                                          >  2:117048/1‑150 (MQ=255)                                                             ccGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAcggcg                                                                                          <  1:1569664/150‑1 (MQ=255)                                                                gCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCt                                                                                       >  1:719394/1‑150 (MQ=255)                                                                 ccGCGCTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTc                                                                                      <  2:422841/150‑1 (MQ=255)                                                                           ggCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCagcagc                                                                            >  1:257089/1‑150 (MQ=255)                                                                             ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCAcc                                                                                                           >  2:1443189/1‑117 (MQ=255)                                                                             ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCAcc                                                                                                           <  1:1443189/117‑1 (MQ=255)                                                                              cGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGAc                                                                         <  1:893862/150‑1 (MQ=255)                                                                                 aCCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGc                                                                      >  1:692948/1‑150 (MQ=255)                                                                                 aCCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCACTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGc                                                                      >  2:53134/1‑150 (MQ=255)                                                                                     tGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACGCGAAGGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCAtt                                                                  <  1:117048/150‑1 (MQ=255)                                                                                         gTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGAt                                                              <  1:1492846/150‑1 (MQ=255)                                                                                              aTGCGTCGGGGTCTTCGCCACGGTACAGCGCAACGACCGGAGCGTCGCGCCGCGGCGGGA‑CGCCGCGCTGATTGCTGGGGCGCGTAGCCCATATTCCACCCGCGTCGATCGACGGCGGCTCGGCCGGCAGCGACAGCCAATCGTTACGGc                                                         >  1:433121/1‑150 (MQ=255)                                                                                                 cGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTgg                                                      <  2:1846644/150‑1 (MQ=255)                                                                                                  gTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGc                                                     >  1:1225549/1‑150 (MQ=255)                                                                                                   tCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCg                                                    >  1:77070/1‑150 (MQ=255)                                                                                                     ggggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGTA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGtt                                                  <  1:1057100/150‑1 (MQ=255)                                                                                                     ggggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACCGCCATTCGATACGGCTGGCGtt                                                  <  2:719394/150‑1 (MQ=255)                                                                                                     ggggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGtt                                                  >  1:182208/1‑150 (MQ=255)                                                                                                       ggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGc                                                >  2:1353163/1‑150 (MQ=255)                                                                                                        gTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGcc                                               <  1:749150/150‑1 (MQ=255)                                                                                                         tCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGtt                                                  >  1:1299705/1‑146 (MQ=255)                                                                                                         tCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGtt                                                  <  2:1299705/146‑1 (MQ=255)                                                                                                          cTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCtt                                             <  1:1792666/150‑1 (MQ=255)                                                                                                               ccACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACa                                                                        >  2:1711465/1‑118 (MQ=255)                                                                                                               ccACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACa                                                                        <  1:1711465/118‑1 (MQ=255)                                                                                                               ccACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCa                                        >  1:1210520/1‑150 (MQ=255)                                                                                                                 aCGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATc                                      <  2:692948/150‑1 (MQ=255)                                                                                                                 aCGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTATTCTCACCGGCGTCGATCGCCGGAGACTAGAACACCAGCTACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATc                                      >  1:1868393/1‑150 (MQ=255)                                                                                                                      aCAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGt                                 <  1:929989/150‑1 (MQ=255)                                                                                                                        aGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTc                               >  1:605036/1‑150 (MQ=255)                                                                                                                          cgcCTCGACCTGAGCGGGGGGCCGCTGCAGCA‑CGGCGGGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                             <  2:716050/150‑1 (MQ=255)                                                                                                                          cgcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                             <  2:257089/150‑1 (MQ=255)                                                                                                                          cgcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                             >  1:1002209/1‑150 (MQ=255)                                                                                                                             cTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAgcg                          <  1:811679/150‑1 (MQ=255)                                                                                                                                   tcaatcagcgcgcCG‑TGCTGCAGCGGCGCGCTGATTGATCTAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccag                    >  2:1844009/8‑150 (MQ=255)                                                                                                                                    tGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagc                   <  2:92923/150‑1 (MQ=255)                                                                                                                                    tGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagc                   <  1:1200265/150‑1 (MQ=255)                                                                                                                                      aGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagcca                 >  1:714665/1‑150 (MQ=255)                                                                                                                                          gcgcgcCG‑TGCTGCAGCGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagc                   <  2:169577/144‑1 (MQ=255)                                                                                                                                          gcgcgcCG‑TGCTGCAGCGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagc                   >  1:169577/1‑144 (MQ=255)                                                                                                                                            gcgcCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGc           >  2:554754/1‑150 (MQ=255)                                                                                                                                              gccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                             <  2:663077/130‑1 (MQ=255)                                                                                                                                              gccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                             >  1:663077/1‑130 (MQ=255)                                                                                                                                               ccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGaa        <  2:1225549/150‑1 (MQ=255)                                                                                                                                               ccgcGGCCGCC‑CGGCGCGCTGATTGAATGAGCCCTTGCCTGTGATTCCTCTGGCGTTGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGaa        <  2:1451782/150‑1 (MQ=255)                                                                                                                                                   tGCAGCA‑CGGCGCGCTGATTGAGCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCaa    <  2:812825/150‑1 (MQ=255)                                                                                                                                                     cagca‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGTGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACt  >  1:1760910/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TCGACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACT  >  NC_002947/3794543‑3794840

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