Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,075 G→A 5.1% D382N (GAC→AAC)  sdaA → L‑serine dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,0750GA 5.1% 393.9 / 8.7 117D382N (GAC→AAC) sdaAL‑serine dehydratase 1
Reads supporting (aligned to +/- strand):  ref base G (58/53);  new base A (3/3);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

Only 100 of 117 total aligned reads displayed.

AAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACG  >  NC_000913/1897938‑1898219
                                                                                                                                         |                                                                                                                                                
aaGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCAACCGGGtcgc                                                                                                                                        >  3:220477/1‑145 (MQ=255)
    tggttgCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTgcaggg                                                                                                                                    >  3:780031/1‑148 (MQ=255)
            aGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCa                                                                                                                            >  1:902277/1‑148 (MQ=255)
               gCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCa                                                                                                                            <  4:251675/145‑1 (MQ=255)
               gCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCa                                                                                                                            >  3:251675/1‑145 (MQ=255)
               gCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGt                                                                                                                         <  2:902277/148‑1 (MQ=255)
                  aaGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTCCAGGGCAGGTTCAGGTGcc                                                                                                                      >  4:775036/1‑148 (MQ=255)
                       ggTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGTGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCa                                                                                                                 <  1:177105/148‑1 (MQ=255)
                         tgtTGCCTGTTCAAGGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTCGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAtt                                                                                                               >  4:468459/1‑148 (MQ=255)
                                tGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTa                                                                                                        >  1:260398/1‑148 (MQ=255)
                                 gTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTa                                                                                                        <  4:228665/147‑1 (MQ=255)
                                  ttCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCg                                                                                                                     <  1:803368/133‑1 (MQ=255)
                                   tCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCAACCGGGTCgcaggttaa                                                                                                                                 >  1:119811/1‑116 (MQ=255)
                                   tCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCAACCGGGTCgcaggttaa                                                                                                                                 <  2:119811/120‑5 (MQ=255)
                                     aaTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTgca                                                                                                                                       <  1:790252/112‑1 (MQ=255)
                                     aaTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTgca                                                                                                                                       >  2:790252/1‑112 (MQ=255)
                                     aaTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGcc                                                                                                   <  3:468459/148‑1 (MQ=255)
                                      aTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCa                                                                                                  >  1:380035/1‑148 (MQ=255)
                                       tGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCg                                                                                                          >  1:613236/1‑139 (MQ=255)
                                       tGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCg                                                                                                          <  2:613236/139‑1 (MQ=255)
                                       tGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                                                                 >  1:436227/1‑148 (MQ=255)
                                        ggCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCAtt                                                                                                >  3:979036/1‑148 (MQ=255)
                                             ctgcGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctc                                                                                           <  4:780031/148‑1 (MQ=255)
                                                     cTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAgg                                                                                   <  2:436227/148‑1 (MQ=255)
                                                       tGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCa                                                                                                                            <  1:1056690/105‑1 (MQ=255)
                                                       tGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCa                                                                                                                            >  2:1056690/1‑105 (MQ=255)
                                                         cAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtg                                                                                         >  4:973678/1‑138 (MQ=255)
                                                         cAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtg                                                                                         <  3:973678/138‑1 (MQ=255)
                                                           gAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGTAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTa                                                                             <  1:1116977/148‑1 (MQ=255)
                                                           gAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTa                                                                             >  4:67803/1‑148 (MQ=255)
                                                                gctgGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACGACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTaa                                                                                                       >  3:810761/1‑117 (MQ=255)
                                                                gctgGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACGACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTaa                                                                                                       <  4:810761/117‑1 (MQ=255)
                                                                 ctgGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCAACCGGGTCgcaggttaa                                                                                                                                 >  2:69368/1‑86 (MQ=255)
                                                                 ctgGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCAACCGGGTCgcaggttaa                                                                                                                                 <  1:69368/90‑5 (MQ=255)
                                                                   gggCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGa                                                                                <  4:944911/137‑1 (MQ=255)
                                                                   gggCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGa                                                                                >  3:944911/1‑137 (MQ=255)
                                                                     gCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTg                                                                                                              <  1:50646/105‑1 (MQ=255)
                                                                     gCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTg                                                                                                              >  2:50646/1‑105 (MQ=255)
                                                                     gCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctc                                                                                           >  2:723607/1‑124 (MQ=255)
                                                                     gCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctc                                                                                           <  1:723607/124‑1 (MQ=255)
                                                                        gTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTg                                                                                                              >  1:66783/1‑102 (MQ=255)
                                                                        gTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTg                                                                                                              <  2:66783/102‑1 (MQ=255)
                                                                                gAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTgc                                                        <  3:67803/148‑1 (MQ=255)
                                                                                  aCAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCa                                                                                                                            >  2:398476/1‑78 (MQ=255)
                                                                                  aCAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATg                                                               >  2:819746/1‑139 (MQ=255)
                                                                                     ggTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGATCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGcgccgc                                                   <  3:775036/148‑1 (MQ=255)
                                                                                       tttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTa                                                                                                        >  4:857903/1‑93 (MQ=255)
                                                                                       tttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTa                                                                                                        <  3:857903/93‑1 (MQ=255)
                                                                                       tttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                 >  3:329219/1‑148 (MQ=255)
                                                                                        ttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGc                                                                                                    <  4:923032/96‑1 (MQ=255)
                                                                                        ttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGc                                                                                                    >  3:923032/1‑96 (MQ=255)
                                                                                        ttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCg                                                                                 <  2:420732/115‑1 (MQ=255)
                                                                                        ttGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCg                                                                                 >  1:420732/1‑115 (MQ=255)
                                                                                             tGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGt                                                                                                                    <  1:896563/75‑1 (MQ=255)
                                                                                             tGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGt                                                                                                                    >  2:896563/1‑75 (MQ=255)
                                                                                             tGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGc                                           <  2:260398/148‑1 (MQ=255)
                                                                                                   cggAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctct                                                                                          >  3:313837/1‑95 (MQ=255)
                                                                                                   cggAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctct                                                                                          <  4:313837/95‑1 (MQ=255)
                                                                                                   cggAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGAGCGGATGGCTCTGCGCCGCACCAGTGCACcgcg                                     >  4:751762/1‑148 (MQ=255)
                                                                                                     gAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTaa                                                                                                       >  2:210060/1‑80 (MQ=255)
                                                                                                     gAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCg                                                                                 >  1:1096915/1‑102 (MQ=255)
                                                                                                     gAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCg                                                                                 <  2:1096915/102‑1 (MQ=255)
                                                                                                     gAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCCTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgcg                                   <  3:751762/148‑1 (MQ=255)
                                                                                                        aTTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtct                                >  1:103799/1‑148 (MQ=255)
                                                                                                          tGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                              <  3:870110/148‑1 (MQ=255)
                                                                                                            gCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCGGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCt                            >  1:358868/1‑148 (MQ=255)
                                                                                                               tGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAt                                                                                                                <  1:32884/61‑1 (MQ=255)
                                                                                                               tGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAt                                                                                                                >  2:32884/1‑61 (MQ=255)
                                                                                                                 gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTg                                                         >  1:726857/1‑114 (MQ=255)
                                                                                                                 gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCCGATGGCTCTg                                                         <  2:726857/114‑1 (MQ=255)
                                                                                                                  aaCACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGa                                                                                      <  1:276386/84‑1 (MQ=255)
                                                                                                                  aaCACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGa                                                                                      >  2:276386/1‑84 (MQ=255)
                                                                                                                   acacAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTa                                                                                                        <  2:575626/65‑1 (MQ=255)
                                                                                                                   acacAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTa                                                                                                        >  1:575626/1‑65 (MQ=255)
                                                                                                                   acacAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                     <  3:686209/148‑1 (MQ=255)
                                                                                                                     acaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAgg                                                                                   >  4:884571/1‑84 (MQ=255)
                                                                                                                     acaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCg                                                                                 <  3:39109/86‑1 (MQ=255)
                                                                                                                     acaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCg                                                                                 >  4:39109/1‑86 (MQ=255)
                                                                                                                      caACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCa                                                                                                  >  3:432250/1‑68 (MQ=255)
                                                                                                                      caACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCa                                                                                                  <  4:432250/68‑1 (MQ=255)
                                                                                                                        aCCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGt                                                                                                         <  4:411002/59‑1 (MQ=255)
                                                                                                                        aCCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCgg                                                                  <  3:375500/98‑1 (MQ=255)
                                                                                                                        aCCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCgg                                                                  >  4:375500/1‑98 (MQ=255)
                                                                                                                          cTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCg              >  2:302620/1‑148 (MQ=255)
                                                                                                                           ttGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTaa                                                                                                       <  3:39905/58‑1 (MQ=255)
                                                                                                                           ttGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGa             >  4:628462/1‑148 (MQ=255)
                                                                                                                            tGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa            <  3:628462/148‑1 (MQ=255)
                                                                                                                               tttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGctc                                                           >  1:28097/1‑98 (MQ=255)
                                                                                                                               tttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGctc                                                           <  2:28097/98‑1 (MQ=255)
                                                                                                                                ttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCAt                >  4:910492/1‑140 (MQ=255)
                                                                                                                                ttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCAt                <  3:910492/140‑1 (MQ=255)
                                                                                                                                ttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa        >  1:513065/1‑148 (MQ=255)
                                                                                                                                ttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa        >  4:488371/1‑148 (MQ=255)
                                                                                                                                    ccTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATa                       <  1:942332/129‑1 (MQ=255)
                                                                                                                                    ccTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATa                       >  2:942332/1‑129 (MQ=255)
                                                                                                                                    ccTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATgca    <  3:488371/148‑3 (MQ=255)
                                                                                                                                    ccTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTa    <  4:329219/148‑1 (MQ=255)
                                                                                                                                      tGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa            <  4:729577/138‑1 (MQ=255)
                                                                                                                                      tGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACg  >  2:740506/1‑148 (MQ=255)
                                                                                                                                      tGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACg  >  4:110669/1‑148 (MQ=255)
                                                                                                                                         |                                                                                                                                                
AAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACG  >  NC_000913/1897938‑1898219

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