Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	99,900	A→G	5.9%	Q86R (CAG→CGG)	murG →	N‑acetylglucosaminyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	99,900	0	A	G	5.9%	208.8 / 3.0	68	Q86R (CAG→CGG)	murG	N‑acetylglucosaminyl transferase
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

GCTGGCTGGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCAAT  >  NC_000913/99753‑100045                                                                                                                                                    |                                                                                                                                                  gctggctgGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAg                                                                                                                                                  <  5:33567/149‑1 (MQ=255) gctggctgGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAg                                                                                                                                                  <  4:58739/149‑1 (MQ=255)     gctggGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGgcgc                                                                                                                                              >  8:248071/1‑149 (MQ=255)             cTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTGCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGAtt                                                                                                                                      <  1:254206/149‑1 (MQ=255)                             aGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAAc                                                                                                                      <  5:144317/149‑1 (MQ=255)                                ggACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCt                                                                                                                    >  4:233482/1‑148 (MQ=255)                                 gACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACaaa                                                                                                                       >  6:360921/1‑144 (MQ=255)                                 gACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACaaa                                                                                                                       <  5:360921/144‑1 (MQ=255)                                  aCTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGAc                                                                                                                 <  4:195275/149‑1 (MQ=255)                                    ttAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgt                                                                                                               >  5:181006/1‑149 (MQ=255)                                    ttAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGACCCGCTGCGTATCTTCAACGCCTGGCGTCAGGAGCGGGCGATTATGAAATCGTACAAACATGACgt                                                                                                               >  8:176018/1‑149 (MQ=255)                                      aGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgtg                                                                                                              <  7:248071/148‑1 (MQ=255)                                       gTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgtggt                                                                                                            <  5:348036/149‑1 (MQ=255)                                         gCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgtggt                                                                                                            >  8:49988/1‑147 (MQ=255)                                         gCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgtggt                                                                                                            <  7:49988/147‑1 (MQ=255)                                            aaaaCATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCg                                                                                                       >  3:41882/1‑149 (MQ=255)                                            aaaaCATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCg                                                                                                       >  2:268763/1‑149 (MQ=255)                                                  tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATgg                                                                                                 >  7:181503/1‑149 (MQ=255)                                                  tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATgg                                                                                                 <  5:165423/149‑1 (MQ=255)                                                    gCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATgg                                                                                                 >  3:12147/1‑147 (MQ=255)                                                    gCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGa                                                                                               <  7:176018/149‑1 (MQ=255)                                                    gCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGa                                                                                               <  1:268763/149‑1 (MQ=255)                                                       tCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCACGCTGAGTATGTTCAACGCATTGAGTAATGCGAGGTCGATTATGAAAGCGTACAAACCTGACGTGGTGATCGGTATGGGAGGc                                                                                            >  3:81657/1‑149 (MQ=255)                                                         gAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCt                                                                                           <  5:129375/148‑1 (MQ=255)                                                          aaaTTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTAc                                                                                         >  5:62610/1‑149 (MQ=255)                                                            aTTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACgt                                                                                       <  3:233482/149‑1 (MQ=255)                                                                aTTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTc                                                                                    <  4:41882/148‑1 (MQ=255)                                                                 tttCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAg                                                                                  >  1:26112/1‑149 (MQ=255)                                                                  ttCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAgg                                                                                 >  5:140830/1‑149 (MQ=255)                                                                     aTTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTc                                                                               >  1:73317/1‑148 (MQ=255)                                                                           aTCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATgg                                                                                                 <  3:114156/124‑1 (MQ=255)                                                                           aTCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATgg                                                                                                 >  4:114156/1‑124 (MQ=255)                                                                                      gCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGt                                                             >  5:2148/1‑149 (MQ=255)                                                                                       cGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGt                                                             <  8:181503/148‑1 (MQ=255)                                                                                          ggAAAACGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGAGGCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTgggc                                                         <  5:106268/149‑3 (MQ=255)                                                                                           gAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTc                                                         <  6:2148/148‑1 (MQ=255)                                                                                            aaaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCGGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTg                                                                    <  8:314243/136‑1 (MQ=255)                                                                                            aaaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCGGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTg                                                                    >  7:314243/1‑136 (MQ=255)                                                                                               aGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAg                                                    <  6:181006/149‑1 (MQ=255)                                                                                                      aaaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGTTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTcc                                             >  2:98422/1‑149 (MQ=255)                                                                                                      aaaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAggtggt                                                                       <  5:231997/123‑1 (MQ=255)                                                                                                      aaaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAggtggt                                                                       >  6:231997/1‑123 (MQ=255)                                                                                                      aaaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTcc                                             >  3:266078/1‑149 (MQ=255)                                                                                                        aGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTAc                                                                                         <  8:271218/103‑1 (MQ=255)                                                                                                        aGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTAc                                                                                         >  7:271218/1‑103 (MQ=255)                                                                                                         gCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGt                                          >  5:305700/1‑149 (MQ=255)                                                                                                          cACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGtt                                         <  6:140830/149‑1 (MQ=255)                                                                                                             tGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTa                                      >  6:228266/1‑149 (MQ=255)                                                                                                             tGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTa                                      >  1:329253/1‑149 (MQ=255)                                                                                                                tAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCGGGCGCGGGCGa                                                                                                                                        <  6:11409/48‑1 (MQ=255)                                                                                                                tAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCGGGCGCGGGCGa                                                                                                                                        >  5:11409/1‑48 (MQ=255)                                                                                                                   cTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCAt                                <  2:73317/149‑1 (MQ=255)                                                                                                                    tGCCACGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTCGTGCTCGGTTTGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGAGTGGTAGTTATGCATTCAGGTGGTACTTCAt                                >  2:43554/1‑148 (MQ=255)                                                                                                                      ccccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATg                               <  1:98422/147‑1 (MQ=255)                                                                                                                       cccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGaa                             >  1:341995/1‑148 (MQ=255)                                                                                                                       cccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCCATTATGAAAGCGTACAAACCTGACGTGGTGCGCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGaa                             >  7:11223/1‑148 (MQ=255)                                                                                                                        ccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAAc                            >  4:193864/1‑148 (MQ=255)                                                                                                                              cGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGGCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAAc                       <  6:305700/147‑1 (MQ=255)                                                                                                                                      cAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTg                <  2:26112/146‑1 (MQ=255)                                                                                                                                      cAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCg              <  4:266078/148‑1 (MQ=255)                                                                                                                                       aaCGCCTGGCGTCAGGCGCGGGCGATTATGAAACCGTACAACCCTGAACTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCAGGCCGCGTGGTCGTTAGGCATTCCGGTTGTCCTTCATGAACAAAACTGTATTGCggg            >  6:301525/1‑149 (MQ=255)                                                                                                                                           ccTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTc                                              >  3:78586/1‑111 (MQ=255)                                                                                                                                           ccTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTc                                              <  4:78586/111‑1 (MQ=255)                                                                                                                                               gCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAAc      <  6:62610/147‑1 (MQ=255)                                                                                                                                                 gTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCaa   >  5:49025/1‑148 (MQ=255)                                                                                                                                                 gTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCAAt  <  2:329253/149‑1 (MQ=255)                                                                                                                                                  tCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGt                                       >  8:14953/1‑111 (MQ=255)                                                                                                                                                  tCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGt                                       <  7:14953/111‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCTGGCTGGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCAAT  >  NC_000913/99753‑100045

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