Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,666,886	T→C	S321G (AGC→GGC)	hcaT ←	putative 3‑phenylpropionic transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,666,886	0	T	C	98.0%	348.6 / ‑2.6	100	S321G (AGC→GGC)	hcaT	putative 3‑phenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (46/52);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TATTGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGC  >  NC_000913/2666792‑2666982                                                                                               |                                                                                                 tatTGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCt                                                                                             <  1:243615/100‑1 (MQ=255)   ttGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTgg                                                                                           <  1:2759150/100‑1 (MQ=255)   ttGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTgg                                                                                           >  1:48787/1‑100 (MQ=255)    tGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGc                                                                                          >  1:2641396/1‑100 (MQ=255)    tGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGc                                                                                          >  2:237327/1‑100 (MQ=255)    tGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGc                                                                                          <  2:486993/100‑1 (MQ=255)     gATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCg                                                                                         <  2:2190304/100‑1 (MQ=255)      aTACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGa                                                                                        <  2:326017/100‑1 (MQ=255)       tACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGag                                                                                       >  1:969190/1‑100 (MQ=255)       tACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGag                                                                                       <  1:1501742/100‑1 (MQ=255)         cAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGagca                                                                                     <  1:305228/100‑1 (MQ=255)          aGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGagcag                                                                                    <  1:2006561/100‑1 (MQ=255)          aGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGagcag                                                                                    <  2:2490263/100‑1 (MQ=255)           ggAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGagcatc                                                                                   <  1:2094052/100‑3 (MQ=255)           ggAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGagcagc                                                                                   <  2:1512481/100‑1 (MQ=255)            gAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAgcagca                                                                                  <  2:1947744/100‑1 (MQ=255)              aaCCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAat                                                                                <  2:2111103/100‑1 (MQ=255)               aCCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAata                                                                               >  1:290791/1‑100 (MQ=255)                ccGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAatat                                                                              >  2:1771807/1‑100 (MQ=255)                ccGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAatat                                                                              <  2:1868760/100‑1 (MQ=255)                  ggCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGCCTTCGCCACCCTGGCGAGCAGCAATATaa                                                                            >  2:1264562/1‑100 (MQ=255)                  ggCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAAtata                                                                             <  2:3088750/99‑1 (MQ=255)                  ggCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATaa                                                                            >  2:2005043/1‑100 (MQ=255)                  ggCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATaa                                                                            >  2:2721491/1‑100 (MQ=255)                   gCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAAtatatc                                                                           <  2:426745/100‑3 (MQ=255)                    cGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACg                                                                          >  2:2539404/1‑100 (MQ=255)                    cGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACg                                                                          <  1:770683/100‑1 (MQ=255)                    cGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACg                                                                          >  1:344360/1‑100 (MQ=255)                     gAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGGTGACTTCGCCACCCTGGCGAGCAGCAATATAACGc                                                                         <  2:2934838/100‑1 (MQ=255)                     gAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGc                                                                         >  2:496610/1‑100 (MQ=255)                          cGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggc                                                                    <  1:1329687/100‑1 (MQ=255)                          cGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggc                                                                    <  2:2039778/100‑1 (MQ=255)                           ggTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggcg                                                                   >  1:1463499/1‑100 (MQ=255)                            gTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggcgg                                                                  <  2:1121067/100‑1 (MQ=255)                            gTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggcgg                                                                  <  2:913710/100‑1 (MQ=255)                             tCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggcggc                                                                 >  2:1538649/1‑100 (MQ=255)                             tCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATggcgg                                                                  <  2:1956033/99‑1 (MQ=255)                               aTGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCa                                                               <  2:689623/100‑1 (MQ=255)                                tGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAg                                                              <  1:2209322/100‑1 (MQ=255)                                tGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAg                                                              <  1:1906552/100‑1 (MQ=255)                                tGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAg                                                              <  2:1102321/100‑1 (MQ=255)                                  aTAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAACCGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGt                                                            >  2:1897809/1‑100 (MQ=255)                                  aTAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGt                                                            >  2:2963377/1‑100 (MQ=255)                                   tAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGTGGCCAGGTg                                                           <  1:2712250/100‑1 (MQ=255)                                   tAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTg                                                           >  2:2236719/1‑100 (MQ=255)                                   tAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGCATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTg                                                           <  1:1128064/100‑1 (MQ=255)                                     cgatctcTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGc                                                         <  1:2211753/94‑1 (MQ=255)                                      cGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCa                                                        <  1:1779292/100‑1 (MQ=255)                                       gATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAg                                                       >  2:2969454/1‑100 (MQ=255)                                       gATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAg                                                       <  1:637213/100‑1 (MQ=255)                                        aTACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGa                                                      >  1:419978/1‑100 (MQ=255)                                           cTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGatcg                                                   >  2:2211731/1‑97 (MQ=255)                                            tGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGt                                                  >  2:2986742/1‑100 (MQ=255)                                             gccgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTg                                                 >  1:1236650/1‑100 (MQ=255)                                              ccgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGa                                                <  1:2721517/100‑1 (MQ=255)                                               cgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGaa                                               <  2:2364462/100‑1 (MQ=255)                                                gccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAg                                              >  1:1748628/1‑100 (MQ=255)                                                gccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAg                                              >  2:561195/1‑100 (MQ=255)                                                 cccATCGCGACGGCAGAGTAACCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAgg                                             >  1:522386/1‑100 (MQ=255)                                                 cccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAgg                                             >  2:3178/1‑100 (MQ=255)                                                 cccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAgg                                             <  1:1922965/100‑1 (MQ=255)                                                 cccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACACATGGCGGCCAGGTGGCAGACCGTGAAgg                                             <  1:23567/100‑1 (MQ=255)                                                     tCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGcc                                         <  1:2815502/100‑1 (MQ=255)                                                     tCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGcc                                         <  1:1547219/100‑1 (MQ=255)                                                      cgcgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCg                                        >  2:267639/1‑100 (MQ=255)                                                       gcgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGc                                       <  1:555078/100‑1 (MQ=255)                                                       gcgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGc                                       >  1:1240697/1‑100 (MQ=255)                                                       gcgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGc                                       >  1:562194/1‑100 (MQ=255)                                                        cgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCa                                      >  1:1112140/1‑100 (MQ=255)                                                        cgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCa                                      >  1:1395544/1‑100 (MQ=255)                                                           cGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATg                                   <  2:2801986/100‑1 (MQ=255)                                                             gCAGAGTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCa                                 <  2:2536930/100‑1 (MQ=255)                                                                gagTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGaa                              <  1:2005065/100‑1 (MQ=255)                                                                 agTAAACCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAAt                             >  1:303940/1‑100 (MQ=255)                                                                    aaaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTg                          <  2:168452/100‑1 (MQ=255)                                                                    aaaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTg                          <  1:1057332/100‑1 (MQ=255)                                                                    aaaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTg                          <  1:2969484/100‑1 (MQ=255)                                                                     aaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGc                         >  1:135601/1‑100 (MQ=255)                                                                     aaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGc                         >  2:587100/1‑100 (MQ=255)                                                                     aaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGc                         >  2:587084/1‑100 (MQ=255)                                                                     aaCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGc                         <  1:2254921/100‑1 (MQ=255)                                                                      aCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGca                        >  2:2182111/1‑100 (MQ=255)                                                                      aCCGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGc                         <  1:496611/99‑1 (MQ=255)                                                                       ccGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGcac                       >  2:679446/1‑100 (MQ=255)                                                                       ccGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGcac                       >  2:678958/1‑100 (MQ=255)                                                                        cGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGcacc                      <  1:897792/100‑1 (MQ=255)                                                                        cGCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGcacc                      <  2:344358/100‑1 (MQ=255)                                                                         gCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCGGAATTTGcacca                     <  1:2031246/100‑1 (MQ=255)                                                                         gCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGcacca                     >  1:2918173/1‑100 (MQ=255)                                                                         gCCTGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGcacc                      <  2:1399141/99‑1 (MQ=255)                                                                            tGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTa                  <  2:562192/100‑1 (MQ=255)                                                                            tGTAAACGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTa                  >  2:837260/1‑100 (MQ=255)                                                                               aaaCGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCa               <  2:2881143/100‑1 (MQ=255)                                                                                aaCGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCaa              >  1:531543/1‑100 (MQ=255)                                                                                  cGGATGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAAcc            <  2:1463485/100‑1 (MQ=255)                                                                                     aTGACTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACg         >  2:1389270/1‑100 (MQ=255)                                                                                      tGACTTCGCCACCCTGGCGCGCAGCAATATAACGCATGGCTTCTAGGTGGCATACCGTGAAGGCGCCGCAACCCAGAATTTGCACCACCATCTGCCAc          >  2:2494600/1‑98 (MQ=255)                                                                                        aCTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCa      >  1:1925939/1‑100 (MQ=255)                                                                                         cTTCGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCaa     >  1:2436935/1‑100 (MQ=255)                                                                                            cGCCACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGc  >  2:112622/1‑100 (MQ=255)                                                                                               |                                                                                                 TATTGATACAGGAAACCGGCGAAAACGGTCATGATAGCGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGC  >  NC_000913/2666792‑2666982

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