Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,579,569	G→A	5.5%	A302T (GCC→ACC)	talA →	transaldolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,569	0	G	A	5.5%	239.6 / 3.3	73	A302T (GCC→ACC)	talA	transaldolase A
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base A (2/2);  total (37/36)
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.30e-01

AAAGTTTCGCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTGGTCATCC  >  NC_000913/2579425‑2579712                                                                                                                                                 |                                                                                                                                                aaaGTTTCGCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGt                                                                                                                                             >  3:330840/1‑149 (MQ=255) aaaGTTTCGCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGt                                                                                                                                             <  4:128553/149‑1 (MQ=255)    gTTTCGCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGa                                                                                                                                          >  1:814519/1‑149 (MQ=255)    gTTTCGCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGa                                                                                                                                          <  2:370824/149‑1 (MQ=255)         gCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAAc                                                                                                                                     <  4:532700/149‑1 (MQ=255)                gTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAAc                                                                                                                                     <  1:221281/142‑1 (MQ=255)                gTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAAc                                                                                                                                     >  2:221281/1‑142 (MQ=255)                    gTAAATTAATCCCACCTTCTCAGACGTTCCCACTCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGaa                                                                                                                          <  4:529007/149‑1 (MQ=255)                    gTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGaa                                                                                                                          <  2:301643/149‑1 (MQ=255)                        aTTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGAtc                                                                                                                      >  1:792317/1‑149 (MQ=255)                         ttAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTg                                                                                                                                           >  4:282521/1‑127 (MQ=255)                         ttAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTg                                                                                                                                           <  3:282521/127‑1 (MQ=255)                         ttAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTg                                                                                                                                           <  3:653125/127‑1 (MQ=255)                         ttAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTg                                                                                                                                           >  4:653125/1‑127 (MQ=255)                            aTCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGAtcttct                                                                                                                  >  4:334353/1‑149 (MQ=255)                             tCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATcttctt                                                                                                                 >  2:330991/1‑149 (MQ=255)                             tCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATcttctt                                                                                                                 >  4:674030/1‑149 (MQ=255)                               ccACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTgc                                                                                                               >  2:98097/1‑149 (MQ=255)                               ccACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTgc                                                                                                               >  4:594344/1‑149 (MQ=255)                                cACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTgcc                                                                                                              >  2:378659/1‑149 (MQ=255)                                  ccTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTgccgc                                                                                                            <  3:594344/149‑1 (MQ=255)                                       tcAGACGTTCCCACGCCCAGCTCCCATGAGAGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAAc                                                                                                       <  1:101725/149‑1 (MQ=255)                                          gACGTTCCCACGCACAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAAACAGGAGGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGGCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGACGCCAAACtat                                                                                                    <  1:602888/149‑1 (MQ=255)                                              ttCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAAc                                                                                                <  2:675619/149‑1 (MQ=255)                                                 ccACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAg                                                                                             >  2:810464/1‑149 (MQ=255)                                                  cACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGc                                                                                            >  4:283410/1‑149 (MQ=255)                                                   aCGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGcc                                                                                           >  2:708143/1‑149 (MQ=255)                                                           cTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAgtg                                                                                   <  4:85539/149‑1 (MQ=255)                                                             cccATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGtt                                                                                 <  1:788968/149‑1 (MQ=255)                                                             cccATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGtt                                                                                 <  1:710938/149‑1 (MQ=255)                                                              ccATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTta                                                                                >  1:41796/1‑149 (MQ=255)                                                              ccAGGAGCGAAGCGGAGTTCCGTTGGGATCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTta                                                                                >  2:773602/1‑149 (MQ=255)                                                                 tGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTatat                                                                             <  2:792317/149‑1 (MQ=255)                                                                   aGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGt                                                                           >  4:751681/1‑149 (MQ=255)                                                                    gCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACtata                                                                                                   <  1:807718/124‑1 (MQ=255)                                                                    gCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACtata                                                                                                   >  2:807718/1‑124 (MQ=255)                                                                         gCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGaa                                                                     >  1:450419/1‑149 (MQ=255)                                                                         gCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCACCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGaa                                                                     >  4:580043/1‑149 (MQ=255)                                                                         gCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCACCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGaa                                                                     >  4:33391/1‑149 (MQ=255)                                                                            gAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGCAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTACAAACCAGCCACGGAGTGTTATATGTCCCGAAAAg                                                                  >  1:665465/1‑149 (MQ=255)                                                                                 ccGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAGGACCtt                                                             <  4:717771/149‑1 (MQ=255)                                                                                 ccGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCACCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATGAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCtt                                                             <  3:580043/149‑1 (MQ=255)                                                                                 ccGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCACCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCtt                                                             <  3:33391/149‑1 (MQ=255)                                                                                    ttGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGcc                                                          >  4:617497/1‑149 (MQ=255)                                                                                        gAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCa                                                                                              <  3:279348/109‑1 (MQ=255)                                                                                        gAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCa                                                                                              >  4:279348/1‑109 (MQ=255)                                                                                              aaTCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTc                                                >  2:448254/1‑149 (MQ=255)                                                                                                tCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTcgc                                              <  1:378659/149‑1 (MQ=255)                                                                                                    gATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTtat                                                                               <  1:290522/112‑1 (MQ=255)                                                                                                    gATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTtat                                                                               >  2:290522/1‑112 (MQ=255)                                                                                                    gATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCAc                                          >  4:45754/1‑149 (MQ=255)                                                                                                    gATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCAc                                          >  1:147112/1‑149 (MQ=255)                                                                                                      tGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTc                                        <  2:41796/149‑1 (MQ=255)                                                                                                        cGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAg                                      <  1:810464/149‑1 (MQ=255)                                                                                                        cGATGGCGGAAGAAAAACTGTCTGAAGGCATTCGTCCGATCGCAGTAGATCAACGCAAACTGGAAGATCTCCTTGCCGCCAAACTACAAACCAGCCACGGAGTTTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAg                                      <  2:657297/149‑1 (MQ=255)                                                                                                          aTGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAg                                                                                     >  1:115330/1‑100 (MQ=255)                                                                                                          aTGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAg                                                                                     <  2:115330/100‑1 (MQ=255)                                                                                                              cGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGa                                <  2:814519/149‑1 (MQ=255)                                                                                                                  aGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCg                            <  3:334353/149‑1 (MQ=255)                                                                                                                  aGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGCAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCg                            <  1:330991/149‑1 (MQ=255)                                                                                                                     aaaaCTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTa                         >  2:789169/1‑149 (MQ=255)                                                                                                                           gTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACaaaaa                   <  1:98097/149‑1 (MQ=255)                                                                                                                             cTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGc                 >  3:21952/1‑149 (MQ=255)                                                                                                                              tGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGcc                >  1:667722/1‑149 (MQ=255)                                                                                                                              tGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGCATGGATGCGGTACAAAAAGcc                <  2:450419/149‑1 (MQ=255)                                                                                                                                aaGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCaa              <  3:538258/149‑1 (MQ=255)                                                                                                                                 aGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAAc             <  3:674030/149‑1 (MQ=255)                                                                                                                                     aTTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTg         >  4:487835/1‑149 (MQ=255)                                                                                                                                     aTTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTg         >  4:733081/1‑149 (MQ=255)                                                                                                                                     aTTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTg         >  2:742234/1‑149 (MQ=255)                                                                                                                                      ttCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTgg        <  1:448254/149‑1 (MQ=255)                                                                                                                                       tCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTGGt       <  4:440590/149‑1 (MQ=255)                                                                                                                                            tGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCAGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTGGTCATcc  >  4:727975/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                AAAGTTTCGCCAGTGGTACGTAAATTAATCCCACCTTCTCAGACGTTCCCACGCCCAGCTCCCATGAGCGAAGCGGAGTTCCGTTGGGAGCACAATCAGGATGCGATGGCGGTAGAAAAACTGTCTGAAGGCATTCGTCTGTTCGCCGTTGATCAACGCAAACTGGAAGATCTTCTTGCCGCCAAACTATAAACCAGCCACGGAGTGTTATATGTCCCGAAAAGACCTTGCCAATGCGATTCGCGCACTCAGTATGGATGCGGTACAAAAAGCCAACTCTGGTCATCC  >  NC_000913/2579425‑2579712

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