Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,523,304	G→T	P264H (CCT→CAT)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,523,304	0	G	T	93.6%	316.6 / 10.7	94	P264H (CCT→CAT)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (43/45);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ATCGCTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAG  >  NC_000913/4523208‑4523399                                                                                                 |                                                                                                aTCGCTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCt                                                                                              >  1:13100/1‑100 (MQ=255) aTCGCTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCt                                                                                              >  1:162626/1‑100 (MQ=255)  tCGCTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTc                                                                                             <  1:1292915/100‑1 (MQ=255)    gCTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCAc                                                                                           >  1:1658500/1‑100 (MQ=255)    gCTTTATCGGTGAGAATTTCCCGGGTGGGGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCAc                                                                                           >  1:1385633/1‑100 (MQ=255)     cTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCAcc                                                                                          >  1:1815580/1‑100 (MQ=255)     cTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCAcc                                                                                          <  2:1934972/100‑1 (MQ=255)      tttATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCg                                                                                         <  1:1343085/100‑1 (MQ=255)           cGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAgg                                                                                    >  2:1116961/1‑100 (MQ=255)            ggggagaATTTCCCGGGTGGTGATGCCTTTTTGGCTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCGGGCGATCTCCCGCCACACATGCTCACCGGAAggg                                                                                   <  2:2226822/97‑1 (MQ=255)            ggTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAggg                                                                                   <  2:1458925/100‑1 (MQ=255)            ggTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAggg                                                                                   <  1:1994179/100‑1 (MQ=255)            ggTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAggg                                                                                   >  2:2151624/1‑100 (MQ=255)            ggTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACAGGCTCACCGGAAggg                                                                                   >  2:786122/1‑100 (MQ=255)               gagaATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGcc                                                                                <  2:1934993/100‑1 (MQ=255)                 gaATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAg                                                                              <  2:1523006/100‑1 (MQ=255)                 gaATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAg                                                                              >  2:2139465/1‑100 (MQ=255)                 gaATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAg                                                                              >  2:80146/1‑100 (MQ=255)                 gaATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAg                                                                              >  2:1145277/1‑100 (MQ=255)                   aTTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACAGGCTCACCGGAAGGGGCCAggg                                                                            >  1:1110189/1‑100 (MQ=255)                      tCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTg                                                                         >  2:453211/1‑100 (MQ=255)                         cGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGTGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAAt                                                                      <  2:1491282/100‑1 (MQ=255)                          gggTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAAtg                                                                     <  2:1761270/100‑1 (MQ=255)                          gggTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCGCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAAt                                                                      <  2:897459/99‑1 (MQ=255)                             tggtgATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTgg                                                                  <  1:1321432/100‑1 (MQ=255)                              ggtgATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTggg                                                                 <  1:2176584/100‑1 (MQ=255)                              ggtgATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTggg                                                                 <  2:1385630/100‑1 (MQ=255)                                tgATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGAt                                                               <  1:1793006/100‑1 (MQ=255)                                tgATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGAt                                                               <  1:720053/100‑1 (MQ=255)                                tgATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGAt                                                               >  1:1379541/1‑100 (MQ=255)                                   tGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGc                                                            >  1:1816530/1‑100 (MQ=255)                                   tGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGc                                                            >  1:274292/1‑100 (MQ=255)                                   ngCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGc                                                            >  2:2274532/2‑100 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  2:988783/100‑1 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  2:274290/100‑1 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  2:163497/100‑1 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  1:81382/100‑1 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  2:1816525/100‑1 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           >  2:676945/1‑100 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  1:1397519/100‑1 (MQ=255)                                    gCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGcc                                                           <  2:1062267/100‑1 (MQ=255)                                      cTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAg                                                         <  1:1532451/100‑1 (MQ=255)                                      cTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAg                                                         >  1:1881563/1‑100 (MQ=255)                                      cTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAg                                                         >  2:568770/1‑100 (MQ=255)                                          ttGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTccc                                                     <  2:1815575/100‑1 (MQ=255)                                          ttGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTccc                                                     >  2:1670670/1‑100 (MQ=255)                                           tGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCa                                                    >  1:640431/1‑100 (MQ=255)                                           tGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCa                                                    >  1:1081069/1‑100 (MQ=255)                                            gACTCAGGTTCAGCGCAGCTCGCGCGGCAGCTCGGGCGACCTCCGCGCACGCATGCCCACCGGAAGGGCCCGGGGCTGACTGTGGGCTTGCCAGTCCCa                                                    >  2:508387/1‑99 (MQ=255)                                            gACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCaa                                                   >  1:2101975/1‑100 (MQ=255)                                             aCTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAAt                                                  <  2:1856239/100‑1 (MQ=255)                                             aCTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAAt                                                  >  1:348234/1‑100 (MQ=255)                                             aCTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAAt                                                  <  1:2139472/100‑1 (MQ=255)                                               tCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATcc                                                >  2:456251/1‑100 (MQ=255)                                                  ggTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTc                                             >  2:1750869/1‑100 (MQ=255)                                                   gTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCCCCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAACCCTTCg                                            >  2:468978/1‑100 (MQ=255)                                                    ttCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCgg                                           <  1:1317147/100‑1 (MQ=255)                                                    ttCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCgg                                           <  1:1317258/100‑1 (MQ=255)                                                    ttCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCgg                                           <  1:1670675/100‑1 (MQ=255)                                                      cAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCTGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCGATCCTTCGGcc                                         <  2:1233202/100‑1 (MQ=255)                                                      cAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGcc                                         <  1:453213/100‑1 (MQ=255)                                                      cAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGcc                                         <  2:1379538/100‑1 (MQ=255)                                                       tctgcAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGcca                                        <  2:627258/97‑1 (MQ=255)                                                       aGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGcca                                        <  1:537195/100‑1 (MQ=255)                                                          gcAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGccaaga                                     >  1:627259/1‑97 (MQ=255)                                                             gCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCt                                  <  2:876153/100‑1 (MQ=255)                                                              cTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTg                                 <  2:1055942/100‑1 (MQ=255)                                                              cTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTg                                 <  2:1658495/100‑1 (MQ=255)                                                              cTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTg                                 >  1:1232110/1‑100 (MQ=255)                                                              cTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTg                                 >  1:2117554/1‑100 (MQ=255)                                                              cTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTg                                 <  2:162624/100‑1 (MQ=255)                                                                     ggAAGCTCTGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTc                          >  1:1537876/1‑100 (MQ=255)                                                                          ctctGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGc                     <  2:1601714/100‑1 (MQ=255)                                                                            ctGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAg                   >  1:199372/1‑100 (MQ=255)                                                                            ctGGCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAg                   >  2:470315/1‑100 (MQ=255)                                                                            ctGGCGATCGCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAg                   <  1:551652/100‑1 (MQ=255)                                                                              ggCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTg                 >  2:256299/1‑100 (MQ=255)                                                                              ggCGATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTg                 >  1:10443/1‑100 (MQ=255)                                                                                ngATCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGcc               >  2:670239/2‑100 (MQ=255)                                                                                cGATCTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGcc               <  2:1110187/100‑1 (MQ=255)                                                                                cGACCTCCCGCCACACAGGCTCACCGGACGGTGCCAGGGCTTAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTg                 <  2:1586613/98‑1 (MQ=255)                                                                                  aTCTCCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCa             >  1:940051/1‑100 (MQ=255)                                                                                    ctcCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCaaa           >  1:1760558/1‑100 (MQ=255)                                                                                     tcCCGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCaaaa          <  1:520761/100‑1 (MQ=255)                                                                                      cccGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCaaaaa         >  1:1011786/1‑100 (MQ=255)                                                                                      cacGCCATACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCaaaaa         >  1:2172209/3‑100 (MQ=255)                                                                                      cacGCCATACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCaaaaa         >  1:1022218/3‑100 (MQ=255)                                                                                       ccGCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAAt        <  1:1720995/100‑1 (MQ=255)                                                                                         gCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTg      <  1:2274540/100‑1 (MQ=255)                                                                                         gCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTg      <  2:675241/100‑1 (MQ=255)                                                                                         gCCACACATGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTg      <  1:1750874/100‑1 (MQ=255)                                                                                           cacacaTGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGAc    <  1:2151631/100‑1 (MQ=255)                                                                                             cacaTGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAg  >  2:1196039/1‑100 (MQ=255)                                                                                             cacaTGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAg  >  1:625403/1‑100 (MQ=255)                                                                                                 |                                                                                                ATCGCTTTATCGGTGAGAATTTCCCGGGTGGTGATGCCTTTTTGACTCAGGTTCAGCGCAGCTCGCGCGGAAGCTCTGGCGATCTCCCGCCACACAGGCTCACCGGAAGGGGCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAG  >  NC_000913/4523208‑4523399

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