Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,982,373	A→T	5.2%	intergenic (+101/‑2)	yifK → / → argX	putative APC family amino acid transporter/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,982,373	0	A	T	5.2%	268.4 / 5.0	77	intergenic (+101/‑2)	yifK/argX	putative APC family amino acid transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (30/43);  new base T (2/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

TTTGGCCTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAGTGGTGGCT  >  NC_000913/3982228‑3982515                                                                                                                                                  |                                                                                                                                               tttGGCCTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAgc                                                                                                                                             <  2:289113/149‑1 (MQ=255)  ttGGCCTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAgcg                                                                                                                                            <  4:635275/149‑1 (MQ=255)      ccctAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCg                                                                                                                                        >  1:82034/4‑149 (MQ=255)      ccTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTTAGCGccgt                                                                                                                                        <  2:529921/149‑3 (MQ=255)       cTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTTAGCGccgtt                                                                                                                                       <  1:269651/149‑4 (MQ=255)       cTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGt                                                                                                                                       >  4:571287/1‑149 (MQ=255)        ttAGTCGCCACGGGAAAGTGCATAAACTGGAGGAATAAGCAGCAAAACGCACCAGCCGGAACCAAACGCGTAATTTATTTAAAAAGGGACTAGACAGAGGGGAGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTa                                                                                                                                      <  2:26884/149‑1 (MQ=255)         tAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAg                                                                                                                                     <  3:409039/149‑1 (MQ=255)           aTCGCCACGGGAAAGCGCATAAACTTGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCt                                                                                                                                   <  3:562466/149‑1 (MQ=255)                cACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCt                                                                                                                              >  1:49240/1‑149 (MQ=255)                 aCGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTg                                                                                                                             <  3:570225/149‑1 (MQ=255)                   gggAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGa                                                                                                                           >  2:366394/1‑149 (MQ=255)                    ggAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGAt                                                                                                                          <  3:635858/149‑1 (MQ=255)                    ggAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGAt                                                                                                                          >  2:233871/1‑149 (MQ=255)                                aaCTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGccc                                                                                                              >  1:24155/1‑149 (MQ=255)                                    ggaggaATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCg                                                                                                          <  3:314948/149‑1 (MQ=255)                                      aggaATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGa                                                                                                        <  3:512036/149‑1 (MQ=255)                                      aggaATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGa                                                                                                        <  1:478605/149‑1 (MQ=255)                                      aggaATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGa                                                                                                        <  2:659614/149‑1 (MQ=255)                                         aaTAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGc                                                                                                     <  2:455824/149‑1 (MQ=255)                                            aaGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCaga                                                                                                  <  2:544929/149‑1 (MQ=255)                                             agcagcAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCagag                                                                                                 <  3:270990/149‑1 (MQ=255)                                             agcagcAAAACGCACAAACCGGAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGATGGATAGAGCGCTGCCCTCCGGAGGCagag                                                                                                 >  2:481889/1‑149 (MQ=255)                                                  caAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGtctc                                                                                            <  3:326799/149‑1 (MQ=255)                                                    aaaCGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAg                                                                                          <  1:233871/149‑1 (MQ=255)                                                    aaaCGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAg                                                                                          <  3:99303/149‑1 (MQ=255)                                                    aaaCGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAg                                                                                          >  3:676344/1‑149 (MQ=255)                                                     aaCGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAgg                                                                                         >  2:265587/1‑149 (MQ=255)                                                      aCGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGt                                                                                        <  1:366394/149‑1 (MQ=255)                                                        gCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTc                                                                                      <  4:641820/149‑1 (MQ=255)                                                          acaAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGa                                                                                    >  3:17019/1‑149 (MQ=255)                                                            aaaCCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAAt                                                                                  >  4:739326/1‑149 (MQ=255)                                                                cGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTaa                                                                                                                                               >  1:492459/1‑84 (MQ=255)                                                                cGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTaa                                                                                                                                               <  2:492459/84‑1 (MQ=255)                                                                  tAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTc                                                                                      >  4:642394/1‑139 (MQ=255)                                                                  tAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTc                                                                                      <  3:642394/139‑1 (MQ=255)                                                                  tAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTc                                                                            <  2:49240/149‑1 (MQ=255)                                                                  tAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTc                                                                            <  1:265587/149‑1 (MQ=255)                                                                  tAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTc                                                                            <  3:384213/149‑1 (MQ=255)                                                                           gcgcAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGcc                                                                   >  1:422148/1‑149 (MQ=255)                                                                              cAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCAtt                                                                <  1:493395/149‑1 (MQ=255)                                                                               aaTTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCAttt                                                               >  1:349346/1‑149 (MQ=255)                                                                               aaTTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCAttt                                                               >  2:682268/1‑149 (MQ=255)                                                                                 tttatttaAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAg                                                             <  1:717140/149‑1 (MQ=255)                                                                                       taAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCgg                                                       >  2:14837/1‑149 (MQ=255)                                                                                        aaaaaGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGgc                                                      >  2:469819/1‑149 (MQ=255)                                                                                         aaaaGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGgcg                                                     >  2:55056/1‑149 (MQ=255)                                                                                          aaaGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGgcgc                                                    >  4:271791/1‑149 (MQ=255)                                                                                            aGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCtt                                                  <  2:349346/149‑1 (MQ=255)                                                                                            aGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCtt                                                  >  4:130811/1‑149 (MQ=255)                                                                                             gggACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGTAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTg                                                 <  2:752978/149‑1 (MQ=255)                                                                                              ggACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGa                                                <  2:24155/149‑1 (MQ=255)                                                                                               gACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAg                                               <  1:682268/149‑1 (MQ=255)                                                                                                 cTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCt                                             <  2:82034/149‑1 (MQ=255)                                                                                                  tAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTg                                            >  3:105428/1‑149 (MQ=255)                                                                                                  tAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTg                                            >  1:478910/1‑149 (MQ=255)                                                                                                   aGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGc                                           <  4:702723/149‑1 (MQ=255)                                                                                                      cAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCggt                                        <  3:749203/149‑1 (MQ=255)                                                                                                       agagGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCggtg                                       >  1:158661/1‑149 (MQ=255)                                                                                                       agagGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCggtg                                       >  1:407710/1‑149 (MQ=255)                                                                                                        gagGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCggtgg                                      <  1:481889/149‑1 (MQ=255)                                                                                                              tGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTaa                                <  2:651749/149‑1 (MQ=255)                                                                                                                ggAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGgtagt                                  <  3:731986/145‑1 (MQ=255)                                                                                                                ggAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGgtagt                                  >  4:731986/1‑145 (MQ=255)                                                                                                                 gAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATAc                             >  2:708727/1‑149 (MQ=255)                                                                                                                    gTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCggt                                        >  3:577519/1‑135 (MQ=255)                                                                                                                    gTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCggt                                        <  4:577519/135‑1 (MQ=255)                                                                                                                    gTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACcgc                          <  2:401151/149‑1 (MQ=255)                                                                                                                     tCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGGTAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACcgcg                         >  4:510681/1‑149 (MQ=255)                                                                                                                           ttatCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACaa                   <  1:14837/149‑1 (MQ=255)                                                                                                                              tCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTAGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGAt                >  2:405598/1‑149 (MQ=255)                                                                                                                               ccACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGTGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGAtt               <  1:505679/149‑1 (MQ=255)                                                                                                                                     ccGCAACGGCGCTAAGCGGGCGTAGCTCAGGTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAgt         <  2:478910/149‑1 (MQ=255)                                                                                                                                     ccGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAgt         <  4:84713/149‑1 (MQ=255)                                                                                                                                     ccGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAgt         >  2:168717/1‑149 (MQ=255)                                                                                                                                       gCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAgtgg       <  1:716704/149‑1 (MQ=255)                                                                                                                                        tacgggCGCTTAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAgtggt      >  3:466205/5‑149 (MQ=255)                                                                                                                                        cAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAgtggt      <  4:336126/149‑1 (MQ=255)                                                                                                                                            ggCGCTTAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAGTGGTGGCt  >  3:163974/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                               TTTGGCCTTAATCGCCACGGGAAAGCGCATAAACTGGAGGAATAAGCAGCAAAACGCACAAACCGTAACCAAACGCGCAATTTATTTAAAAAGGGACTAGACAGAGGGGTGGGAAGTCCGTATTATCCACCCCCGCAACGGCGCTAAGCGCCCGTAGCTCAGCTGGATAGAGCGCTGCCCTCCGGAGGCAGAGGTCTCAGGTTCGAATCCTGTCGGGCGCGCCATTTAGTCCCGGCGCTTGAGCTGCGGTGGTAGTAATACCGCGTAACAAGATTTGTAGTGGTGGCT  >  NC_000913/3982228‑3982515

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