Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,932,096	C→G	16.0%	intergenic (+37/+19)	purT → / ← eda	phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,932,096	0	C	G	16.0%	243.0 / 37.8	94	intergenic (+37/+19)	purT/eda	phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (41/38);  new base G (4/11);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 9.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGG  >  NC_000913/1931997‑1932195                                                                                                    |                                                                                                    agTGTTGTTGACGCCATTGAACGCGCAAAGCACGCCGCCGGACAGATAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCAACAAAAATGCCCGATCCTc                                                                                                     <  1:1786946/100‑1 (MQ=255)  gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCg                                                                                                    >  2:28453/1‑100 (MQ=255)     ttgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATc                                                                                                 >  2:1169622/1‑100 (MQ=255)       gttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                               >  2:1490823/1‑100 (MQ=255)        ttgACGCCATTGAACGCGCGAAGCATGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCggg                                                                                              <  2:722334/100‑1 (MQ=255)         tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgggn                                                                                             <  2:72445/100‑2 (MQ=255)          gACGGCATGGAACGCGCGAAGCACTCCGCCGGACAGGTAAAAGTACAGGGTTAAAGCCGGCAAAAAAACGCTACAATAATGCCCGATCGAGGATCGGGCa                                                                                            <  1:1142515/100‑1 (MQ=255)          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                            <  1:1039702/100‑1 (MQ=255)          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                            <  2:730650/100‑1 (MQ=255)          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                            <  1:1719904/100‑1 (MQ=255)          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                                                            <  2:1183631/100‑1 (MQ=255)          gACGCCATTGAACGCGCGAAGCACGCCGACGGACAGGTAAACGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                            <  1:115702/100‑1 (MQ=255)           aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                            <  2:2193577/99‑1 (MQ=39)           aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                                                           <  2:282192/100‑1 (MQ=255)           aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGATAAACCAGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                                                           >  2:270782/1‑100 (MQ=255)            cGCTATTGACCGCGCGCACCCCGCCGCCAGCTAGATAAAAGTACAGGCTTACACCCGGCAAAAAACCGCTACCAAAATGACCGATCCTCGATCGGGc                                                                                             <  1:1677965/97‑1 (MQ=255)            cGCCATTGAACGCGGGAGGCACGCCGCCGGAAAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAATTGCCCGATCGAGGATCGGGCAt                                                                                           <  1:783538/99‑1 (MQ=255)            cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                          <  2:1187452/100‑1 (MQ=39)             gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                                                         <  1:587840/100‑1 (MQ=255)               cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAttttt                                                                                       <  2:2289165/100‑2 (MQ=38)               cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTg                                                                                       <  1:1117690/100‑1 (MQ=255)                atggAACGCGCGAAGCGCGCCGCCGGACAGGGAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAACAATGCCCGATCGAGGATCGGGCAtttttg                                                                                      <  1:1977102/97‑3 (MQ=37)                aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtttttg                                                                                      <  2:1780023/100‑3 (MQ=38)                  tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAAAGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGAc                                                                                     >  2:1509498/1‑99 (MQ=255)                   gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtt                                                                                   >  2:2023881/1‑100 (MQ=255)                    accgcgccTAGCACGCAGCCGCACAGGTAAACGTACACGCTAAAACAAGGCAACAAAACGCTACAAAAATGCCCTATCCTCGATCGGGCATTTTGACtt                                                                                   <  1:1543587/97‑1 (MQ=255)                    aaCGAGCAAAGTACGACGCAGATCAGGGAAAGCAACATGGTTAAAGCCGACAAAAATACGCTTCAAAAATGCCCGGTCCTCGATCGGGCATTTTGACttt                                                                                  <  1:1801443/100‑1 (MQ=255)                      cgcgcgAAGCACGCCGCCGGACAGGTACAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGCTCCTCGATCGGGCATTTTGACTTTTa                                                                                >  2:1281420/1‑100 (MQ=255)                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                                                                 <  1:2442588/99‑1 (MQ=255)                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                                >  2:264048/1‑100 (MQ=255)                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                                <  2:2380419/100‑1 (MQ=255)                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                                >  2:113885/1‑100 (MQ=255)                      cgcgcgAAGCACGCAGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAAAGCCCGATACTCGATCGGGCATTTAGACTTTTa                                                                                >  1:2297881/1‑100 (MQ=255)                       gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                               >  1:2250764/1‑100 (MQ=255)                       gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                               >  2:1257689/1‑100 (MQ=255)                         gcgAAGCACGCCGGCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGAGACAAAAATGCCCGATCCTCGATCGGGCATTTTGGCTTTTACa                                                                              >  1:518219/1‑99 (MQ=255)                          cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGc                                                                            >  1:288554/1‑100 (MQ=255)                          cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGc                                                                            >  1:1723009/1‑100 (MQ=255)                          cgAAGCACGCCGCAGGACAGGTAAAAGCACAGGGTTAAACACGGCAAAAAAACGCTACAACACTGCCCGCTCCTCGATCGGGCATTTTGACTTTTACAGc                                                                            <  1:1509489/100‑1 (MQ=255)                            aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCtt                                                                          >  2:2237240/1‑100 (MQ=255)                            aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCtt                                                                          >  2:2446728/1‑100 (MQ=255)                              gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                        >  1:1265028/1‑100 (MQ=255)                               cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                       >  2:2222412/1‑100 (MQ=255)                               cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                       >  1:832174/1‑100 (MQ=255)                                aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                      >  2:2301230/1‑100 (MQ=255)                                 cgccgcCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcgc                                                                     >  1:755381/1‑100 (MQ=255)                                     gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCtt                                                                  <  1:2287104/99‑1 (MQ=255)                                     gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                                 <  1:2475731/100‑1 (MQ=255)                                      ccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                                                                >  1:1850827/1‑100 (MQ=255)                                       cgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                                                               <  1:2393080/100‑1 (MQ=255)                                        ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                              <  2:2250779/100‑1 (MQ=255)                                         gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                             <  2:2332132/100‑1 (MQ=255)                                          aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAg                                                            <  1:586697/100‑1 (MQ=255)                                           cAGATACAAGTAGAGGATTGATCCCGGAGAGAAAACGCTACAAAAATTCCCGATCCTCGTTCGTGCATTTTGACTTTTACAGCTTa                                                                         <  1:2493353/86‑1 (MQ=255)                                             ggTAAAAGTACAGTGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                         <  2:184712/100‑1 (MQ=255)                                             ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                         >  1:2565473/1‑100 (MQ=255)                                             ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                         >  1:221866/1‑100 (MQ=255)                                             ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                         >  2:890070/1‑100 (MQ=255)                                              gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                        <  1:1281892/100‑1 (MQ=255)                                               tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa                                                       >  2:291892/1‑100 (MQ=255)                                                aaaaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc                                                      >  1:1864527/1‑100 (MQ=255)                                                   aGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                   <  1:291889/100‑1 (MQ=255)                                                    gTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc                                                  >  1:755673/1‑100 (MQ=255)                                                     tACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGACCGGGCACTTTGACCATTACCGCTTAGCGCCTTCTACAGACTCACCCGcc                                                 >  1:2269390/1‑100 (MQ=255)                                                      aCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                                >  2:1901732/1‑100 (MQ=255)                                                       cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                               >  2:978504/1‑100 (MQ=255)                                                       cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                               >  2:753714/1‑100 (MQ=255)                                                       cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                               >  2:906269/1‑100 (MQ=255)                                                         gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                             >  1:251918/1‑100 (MQ=255)                                                          ggTTAAACCCGGCAAAAAAACGCTACAAATATGCCCGATCCTCGATTGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                            <  1:2347273/100‑1 (MQ=255)                                                           gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                           <  2:86918/100‑1 (MQ=255)                                                               aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATCTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTaa                                       <  2:2461663/100‑1 (MQ=255)                                                                  ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTATCGCCTTCTACAGCTTCACGCCCCAGCTTAGTAATg                                     >  1:1085414/1‑99 (MQ=255)                                                                  ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                    <  2:44914/100‑1 (MQ=255)                                                                    ggCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                                  <  2:1723018/100‑1 (MQ=255)                                                                    ggCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                                  <  1:843243/100‑1 (MQ=255)                                                                     gCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTATAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGt                                 <  2:1810601/100‑1 (MQ=255)                                                                      cAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                                >  2:1024901/1‑100 (MQ=255)                                                                       aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                                <  1:1807566/99‑1 (MQ=255)                                                                           aaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAAt                           >  2:1480376/1‑100 (MQ=255)                                                                            aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                          <  2:668147/100‑1 (MQ=255)                                                                                cTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGccc                      <  1:1173224/100‑1 (MQ=255)                                                                                 tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                     <  2:1500909/100‑1 (MQ=255)                                                                                  aCAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                    >  1:231390/1‑100 (MQ=255)                                                                                   cAAAAATGCCCGATCCTCGATCGGGCATCTTCACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                    <  1:514463/99‑1 (MQ=255)                                                                                    caaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                  >  2:1034210/2‑100 (MQ=38)                                                                                     aaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                 >  2:2012378/1‑100 (MQ=38)                                                                                     aaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                 >  1:2378197/1‑100 (MQ=255)                                                                                       aaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa               <  2:2226189/100‑1 (MQ=255)                                                                                        aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg              >  2:880161/1‑100 (MQ=39)                                                                                        aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg              <  1:1281415/100‑1 (MQ=255)                                                                                         tGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc             >  2:1886242/1‑100 (MQ=255)                                                                                            ccGATCTTCGATCGGGCATTTTGACTTTTACAGCTTTGCGCCTTCTACAGCTTCACGAGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa          <  1:1042410/100‑1 (MQ=255)                                                                                            ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa          >  2:752588/1‑100 (MQ=255)                                                                                                 ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTATAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCTCGCTTCCAGCGCATCTg      <  1:131802/99‑1 (MQ=255)                                                                                                    cGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCgg  <  1:113884/100‑1 (MQ=255)                                                                                                    |                                                                                                    AGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGG  >  NC_000913/1931997‑1932195

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