Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,436,937	A→T	100%	D305E (GAT→GAA)	pdxB ←	erythronate‑4‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,436,937	0	A	T	96.7%	292.0 / 0.5	91	D305E (GAT→GAA)	pdxB	erythronate‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (37/51);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 8.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTATAA  >  NC_000913/2436841‑2437033                                                                                                 |                                                                                                 aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                               >  1:2151534/1‑100 (MQ=255)  tCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGa                                                                                              >  2:2140415/1‑100 (MQ=255)   cGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAg                                                                                             <  2:492344/100‑1 (MQ=255)    gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                            <  1:528915/100‑1 (MQ=255)       cTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                         >  2:337899/1‑100 (MQ=255)        tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                        <  1:1371033/100‑1 (MQ=255)        tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                        <  1:1421069/100‑1 (MQ=255)         cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                       >  2:957600/1‑100 (MQ=255)         cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                       >  1:602352/1‑100 (MQ=255)          aCCCGGTATCCCGGCGACTATACGCAGCGGTGCGTTATCGCGGCGCACATCATACACCAAAGGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGAGCCa                                                                                      <  2:1676377/100‑1 (MQ=255)           cccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCa                                                                                      <  1:2173922/99‑1 (MQ=255)           cccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACACTATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCAt                                                                                     >  2:1965668/1‑100 (MQ=255)               gTATCCCGGCGACTTTACGCAGAGGTGCGTCATCGCGTCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAg                                                                                 <  2:1733659/100‑1 (MQ=255)                 aTCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCg                                                                               <  2:1042529/100‑1 (MQ=255)                 aTCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCg                                                                               >  2:2014841/1‑100 (MQ=255)                  tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGt                                                                              >  2:749115/1‑100 (MQ=255)                  tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGt                                                                              <  1:395598/100‑1 (MQ=255)                   cccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTa                                                                             >  2:2067014/1‑100 (MQ=255)                   cccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTa                                                                             >  1:761312/1‑100 (MQ=255)                    ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTa                                                                             >  2:1748588/1‑99 (MQ=255)                    ccGGCGACTTTACGCAGCGGTGCGACATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                            >  1:1805388/1‑100 (MQ=255)                      ggCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATg                                                                          <  2:792028/100‑1 (MQ=255)                      ggCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATg                                                                          <  1:1416654/100‑1 (MQ=255)                        cGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCg                                                                        <  1:2208180/100‑1 (MQ=255)                             ttACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCaa                                                                   >  1:821110/1‑100 (MQ=255)                               aCGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAc                                                                 <  1:2055045/100‑1 (MQ=255)                                   aGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTg                                                             <  1:1826007/100‑1 (MQ=255)                                    gCGGTGCGTCATCGCGGCGCACATCATACACCAACTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                            >  2:1159320/1‑100 (MQ=255)                                     cGGTGCGTCATCGCGGCGCACATCATACACCAAGTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgc                                                           <  2:2278844/100‑1 (MQ=255)                                     cGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgc                                                           <  1:81893/100‑1 (MQ=255)                                      ggTGCGTCATCGCGGCGCACATCTTACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                          <  1:378309/100‑1 (MQ=255)                                      ggTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                          >  1:2067419/1‑100 (MQ=255)                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                       <  1:1283573/100‑1 (MQ=255)                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                       >  1:1251783/1‑100 (MQ=255)                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                       >  1:2107944/1‑100 (MQ=255)                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                       >  2:2347522/1‑100 (MQ=255)                                          cGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcagg                                                      <  1:1377326/100‑1 (MQ=255)                                           gTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                                     >  2:690575/1‑100 (MQ=255)                                             cATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggcag                                                   <  1:1953900/100‑1 (MQ=255)                                              aTCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGt                                                  >  2:1920688/1‑100 (MQ=255)                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                              <  2:1940002/100‑1 (MQ=255)                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                              <  2:1889266/100‑1 (MQ=255)                                                    gcgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTa                                            <  2:1858179/100‑1 (MQ=255)                                                     cgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTa                                            <  2:510568/99‑1 (MQ=255)                                                     cgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTAt                                           >  1:130027/1‑100 (MQ=255)                                                      gcACATCGTACACCAGATGCAACATGCTTTTTTGCGTCGGTTATTCGAGCGGGCCATGCAGCGTAATGCGATCAAACTCTGGCGCAGGCAGTAATg                                              <  2:1533604/96‑1 (MQ=255)                                                      gcACATCATACACCAAATGCACCATCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATc                                          <  1:852708/100‑1 (MQ=255)                                                          atcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgc                                      <  1:716888/100‑1 (MQ=255)                                                             atACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCa                                   >  1:2015947/1‑100 (MQ=255)                                                              tACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCAATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCaa                                  >  2:1423200/1‑100 (MQ=255)                                                               acacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAAc                                 >  2:96893/1‑100 (MQ=255)                                                                cacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACg                                >  2:729323/1‑100 (MQ=255)                                                                cacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTACGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACg                                >  1:1260548/1‑100 (MQ=255)                                                                   cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAATCTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                             <  2:1364939/100‑1 (MQ=255)                                                                   cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                             <  2:611266/100‑1 (MQ=255)                                                                    aaaTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                             <  2:1671546/99‑1 (MQ=255)                                                                     aaTGCACCCGCCTTTTCAGCGTCGGTTGTTCGAACGGGCCATCCAGCGTAATGCGACCAAACTCTGCCGCACGCAGTAATGTATCCTCCCCAACTTGc                             >  2:2332638/1‑98 (MQ=255)                                                                     aaTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTg                           >  2:2518088/1‑100 (MQ=255)                                                                      aTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTg                           <  2:377849/99‑1 (MQ=255)                                                                      aTGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGt                          <  1:1126893/100‑1 (MQ=255)                                                                       tGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTg                           <  1:1093462/98‑1 (MQ=255)                                                                       tGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGtt                         >  1:917828/1‑100 (MQ=255)                                                                        gCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTc                        <  1:1398496/100‑1 (MQ=255)                                                                         cACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa                       <  2:130025/100‑1 (MQ=255)                                                                         cACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa                       >  2:1070586/1‑100 (MQ=255)                                                                            cAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGc                    <  1:791450/100‑1 (MQ=255)                                                                            cAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGc                    <  1:2494787/100‑1 (MQ=255)                                                                             aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                   <  1:2110033/100‑1 (MQ=255)                                                                             aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                   >  2:2487838/1‑100 (MQ=255)                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                  <  1:1459834/100‑1 (MQ=255)                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                  <  2:270320/100‑1 (MQ=255)                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                  >  1:802593/1‑100 (MQ=255)                                                                               ccTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCa                 <  1:62335/100‑1 (MQ=255)                                                                                cTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGATCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCaa                <  2:1927532/100‑1 (MQ=255)                                                                                cTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCaa                <  2:1633287/100‑1 (MQ=255)                                                                                 ttttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAAt               <  1:1992694/100‑1 (MQ=255)                                                                                   ttCAGGGTCGGTTGTTCGTGCGGGCGATGCGGCGTAATGCGACCAGAGCCTGGCGCAGGGAGTAATGTATCCAGCGCAACGTGCTGTTCa                       <  2:1721900/90‑1 (MQ=255)                                                                                   ttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATaa             >  1:895629/1‑100 (MQ=255)                                                                                   ttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCGGTTCATGCCCAATaa             >  1:2521142/1‑100 (MQ=255)                                                                                      aGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACt          >  2:366610/1‑100 (MQ=255)                                                                                      aGCGTCGGATGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAACGGACCGTGCTGTGGATGCCCAATAAACt          >  2:1669699/1‑100 (MQ=255)                                                                                       gCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAAc           <  2:2225656/98‑1 (MQ=255)                                                                                       gCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACtt         <  2:1803793/100‑1 (MQ=255)                                                                                       gCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACtt         >  1:1828406/1‑100 (MQ=255)                                                                                        cGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACtt         <  1:1765225/99‑1 (MQ=255)                                                                                        cGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTg        >  2:1518181/1‑100 (MQ=255)                                                                                           cGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCta     <  1:1001248/100‑1 (MQ=255)                                                                                           cGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCta     <  2:1828392/100‑1 (MQ=255)                                                                                            ggTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCtat    <  2:2015451/100‑1 (MQ=255)                                                                                            ggTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCtat    >  2:908722/1‑100 (MQ=255)                                                                                             gttgttCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCtata   <  1:2487853/100‑1 (MQ=255)                                                                                              ttGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTATaa  >  2:1816365/1‑100 (MQ=255)                                                                                                 |                                                                                                 ATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTATAA  >  NC_000913/2436841‑2437033

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