Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,579,154	A→C	8.7%	A163A (GCA→GCC)	talA →	transaldolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,154	0	A	C	8.7%	256.1 / 15.2	79	A163A (GCA→GCC)	talA	transaldolase A
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base C (4/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

TGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTAC  >  NC_000913/2579006‑2579298                                                                                                                                                     |                                                                                                                                                 tggtggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCa                                                                                                                                                  >  1:512810/1‑149 (MQ=255)    tggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACgg                                                                                                                                               <  3:323441/149‑1 (MQ=255)         ttGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCtg                                                                                                                                          >  4:408751/1‑149 (MQ=255)         ttGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCtg                                                                                                                                          >  1:627251/1‑149 (MQ=255)           gTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCtgtg                                                                                                                                        <  3:637913/149‑1 (MQ=255)            tATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGc                                                                                                                                       >  4:551950/1‑149 (MQ=255)              tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGc                                                                                                                                       >  4:64955/1‑147 (MQ=255)              tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGc                                                                                                                                       <  3:64955/147‑1 (MQ=255)                 gCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAGAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                                                                  >  4:479993/1‑149 (MQ=255)                   aacaaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCa                                                                                                                                >  2:39382/1‑149 (MQ=255)                    acaaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAg                                                                                                                               <  3:408751/149‑1 (MQ=255)                    acaaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAg                                                                                                                               <  1:47884/149‑1 (MQ=255)                     caaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCa                                                                                                                                                  >  2:498610/1‑129 (MQ=255)                     caaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCa                                                                                                                                                  <  1:498610/129‑1 (MQ=255)                      aaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAATGAAGATATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGc                                                                                                                             >  3:644455/1‑149 (MQ=255)                       aGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCg                                                                                                                            >  2:42583/1‑149 (MQ=255)                             tGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTc                                                                                                                      <  3:448243/149‑1 (MQ=255)                             tGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTc                                                                                                                      <  1:156200/149‑1 (MQ=255)                                gaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAATAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGa                                                                                                                   <  4:397155/149‑1 (MQ=255)                                    tCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTc                                                                                                               >  2:271032/1‑149 (MQ=255)                                    tCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTc                                                                                                               <  2:93414/149‑1 (MQ=255)                                      aCGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGc                                                                                                             >  1:124966/1‑149 (MQ=255)                                      aCGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTATCTGATTTCGc                                                                                                             >  3:443562/1‑149 (MQ=255)                                           ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaa                                                                                                                                   >  1:645816/1‑118 (MQ=255)                                           ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaa                                                                                                                                   <  2:645816/122‑5 (MQ=255)                                           ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGttt                                                                                                        <  3:606082/149‑1 (MQ=255)                                           ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGttt                                                                                                        >  3:676408/1‑149 (MQ=255)                                                  cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGc                                                                                                 <  3:551950/149‑1 (MQ=255)                                                  cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGc                                                                                                 <  4:92492/149‑1 (MQ=255)                                                  cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGc                                                                                                 <  1:158923/149‑1 (MQ=255)                                                        ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTTGCCGTTTGTCGGGCGTAttt                                                                                           >  3:676333/1‑149 (MQ=255)                                                          cTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTAt                                                                                         <  4:141828/149‑1 (MQ=255)                                                          cTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTAt                                                                                         <  4:104659/149‑1 (MQ=255)                                                          cTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTAt                                                                                         <  1:442792/149‑1 (MQ=255)                                                           ttCGACCTGGGAAGGAGTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTTACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATg                                                                                        >  2:480896/1‑149 (MQ=255)                                                            tCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGa                                                                                       <  2:512810/149‑1 (MQ=255)                                                                      aaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAg                                                                             <  2:477941/149‑1 (MQ=255)                                                                       aggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAgg                                                                            <  4:204624/149‑1 (MQ=255)                                                                         gaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCa                                                                          >  2:85827/1‑149 (MQ=255)                                                                                    gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGAt                                                               <  3:102022/149‑1 (MQ=255)                                                                                    gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGAt                                                               <  1:525298/149‑1 (MQ=255)                                                                                         agagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAcc                                                          >  2:324829/1‑149 (MQ=255)                                                                                            gCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGAATGGTATCAGGCACGCAAGCCGATGGACCCGt                                                       <  4:379080/149‑1 (MQ=255)                                                                                              tGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAt                                                     <  1:53782/149‑1 (MQ=255)                                                                                                   aaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGAACGCAAGCCGATGGACCCGTATgtggt                                                <  2:406843/149‑1 (MQ=255)                                                                                                      gaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGga                                             >  3:63821/1‑149 (MQ=255)                                                                                                          gTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGCTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTTGTGgaagaa                                         >  2:346412/1‑149 (MQ=255)                                                                                                             ttAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaagaa                                         >  3:48459/1‑146 (MQ=255)                                                                                                             ttAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaagaa                                         <  4:48459/146‑1 (MQ=255)                                                                                                                aCTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGTTGGAAGAAGATCCg                                   >  1:427130/1‑149 (MQ=255)                                                                                                                 cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCgg                                  <  2:565214/149‑1 (MQ=255)                                                                                                                 cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCCCCCTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTCGAAGAAGATCCga                                  <  4:442956/149‑2 (MQ=255)                                                                                                                  tGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCTTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCGGTATGTGGTGGAAGAAGATCCggg                                 >  4:673894/1‑149 (MQ=255)                                                                                                                  tGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCggg                                 >  1:421239/1‑149 (MQ=255)                                                                                                                   gCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTTGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGc                                <  4:676333/149‑1 (MQ=255)                                                                                                                        cTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATcc                                    <  3:442654/140‑1 (MQ=255)                                                                                                                        cTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTaa                           >  2:184546/1‑149 (MQ=255)                                                                                                                        atGACGCTGCTGTATTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGCGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATcc                                    >  4:442654/2‑140 (MQ=255)                                                                                                                          gACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAAt                         >  4:504717/1‑149 (MQ=255)                                                                                                                            cGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGGTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGCGCGTTAAATCg                       <  4:63821/149‑1 (MQ=255)                                                                                                                             gctgctGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTaaa                          >  4:248184/1‑145 (MQ=255)                                                                                                                             gctgctGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTaaa                          <  3:248184/145‑1 (MQ=255)                                                                                                                               tgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGt                     >  3:458204/1‑148 (MQ=255)                                                                                                                               tgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGt                     <  4:458204/148‑1 (MQ=255)                                                                                                                                gctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTgc                   <  1:206606/149‑1 (MQ=255)                                                                                                                                 ctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggtg                                               <  1:352948/120‑1 (MQ=255)                                                                                                                                 ctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggtg                                               >  2:352948/1‑120 (MQ=255)                                                                                                                                 ctgCTTTCTTTTGCACAGGCACGGGCCTGTGTGGAAGCAGGCGTTTTTGTGATTTCGCCGTTTGTCCGGCGCATTTATGACTGCTATCAGGCACGCAAGCCGATGGACCAGTAATTGCTGGAAGAAGATGCTGGCGTTGAATCGGTgcg                  >  2:29149/1‑149 (MQ=255)                                                                                                                                  tgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTgggc                 >  1:56538/1‑149 (MQ=255)                                                                                                                                        cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATc           <  4:142081/149‑1 (MQ=255)                                                                                                                                         ttttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt          >  1:439718/1‑149 (MQ=255)                                                                                                                                          tttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTTATTTCGCCGTTTGTCGGGCGTATTTATGAATGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTa         <  2:56538/149‑1 (MQ=255)                                                                                                                                          tttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTa         <  2:627251/149‑1 (MQ=255)                                                                                                                                           ttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTAc        <  3:148611/149‑1 (MQ=255)                                                                                                                                            cgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                              >  4:75257/2‑94 (MQ=255)                                                                                                                                            cgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                              <  3:75257/93‑1 (MQ=255)                                                                                                                                             gCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGa      >  3:542444/1‑149 (MQ=255)                                                                                                                                                caGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGacta   <  3:548114/149‑1 (MQ=255)                                                                                                                                                 aGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGactac  >  4:585564/1‑149 (MQ=255)                                                                                                                                                 aGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGactac  <  1:324829/149‑1 (MQ=255)                                                                                                                                                   gCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaa                                            <  1:554551/105‑1 (MQ=255)                                                                                                                                                   gCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaa                                            >  2:554551/1‑105 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTAC  >  NC_000913/2579006‑2579298

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑