Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,649 A→T 100% I124F (ATC→TTC)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,6490AT88.5% 305.8 / 25.5 104I124F (ATC→TTC) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base T (46/46);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 104 total aligned reads displayed.

ACTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATT  >  NC_000913/4214551‑4214744
                                                                                                  |                                                                                               
aCTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGAt                                                                                                <  1:2476576/100‑1 (MQ=255)
 cTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGAt                                                                                                <  1:113465/99‑1 (MQ=255)
  tttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCa                                                                                              >  1:1282129/1‑100 (MQ=255)
  tttGAAGATATTCAGGATCAGAACTTTGACAGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCa                                                                                              <  2:1974722/100‑1 (MQ=255)
   ttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCaa                                                                                             >  2:9236/1‑100 (MQ=255)
     gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAAc                                                                                           <  1:1003952/100‑1 (MQ=255)
     gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAAc                                                                                           >  1:1523588/1‑100 (MQ=255)
        gATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAgg                                                                                        >  1:11239/1‑100 (MQ=255)
           atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                     >  1:1353157/1‑100 (MQ=255)
           atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                     <  1:648319/100‑1 (MQ=255)
           atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                     <  2:1353150/100‑1 (MQ=255)
           atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                     <  1:2254372/100‑1 (MQ=255)
            ttCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCt                                                                                    <  1:586513/100‑1 (MQ=255)
              cAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCt                                                                                    >  2:2229603/1‑98 (MQ=255)
              cAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTgg                                                                                  >  2:91935/1‑100 (MQ=255)
                ggATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAg                                                                                <  1:302735/100‑1 (MQ=255)
                ggATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAg                                                                                <  2:1359892/100‑1 (MQ=255)
                ggATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAg                                                                                <  1:2308930/100‑1 (MQ=255)
                 gATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGt                                                                               >  2:543434/1‑100 (MQ=255)
                 gATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACACGTTCTGGa                                                                                 >  2:1463272/1‑98 (MQ=255)
                   tCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTgg                                                                             >  2:2077232/1‑100 (MQ=255)
                    cAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGt                                                                            <  2:2331505/100‑1 (MQ=255)
                    cAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGt                                                                            >  1:2126944/1‑100 (MQ=255)
                    cAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGt                                                                            >  2:1544989/1‑100 (MQ=255)
                     aGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTc                                                                           >  1:1930333/1‑100 (MQ=255)
                      gAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCg                                                                          >  2:323654/1‑100 (MQ=255)
                       aaCTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGa                                                                         <  1:172862/100‑1 (MQ=255)
                        aCTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGGTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGaa                                                                        >  1:488721/1‑100 (MQ=255)
                         cTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGaaa                                                                       <  2:2571040/100‑1 (MQ=255)
                          tttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAg                                                                      <  1:1547573/100‑1 (MQ=255)
                           ttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGa                                                                     <  2:1161148/100‑1 (MQ=255)
                             gACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTGGAAAGATc                                                                   >  1:2012610/1‑100 (MQ=255)
                              aCGGCTTGATTGTAACTGGTGCGCCGCTGAGCCTTCTGAAGGTTAATGATTTAACTTACCGGACGCAGTTCAAACAGCTGCTGGAGTAGTCGtaaaacca                                                                  <  2:670280/100‑6 (MQ=255)
                               cGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCAc                                                                 >  1:2490358/1‑100 (MQ=255)
                               cGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCAc                                                                 <  2:206695/100‑1 (MQ=255)
                                ggTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACg                                                                <  2:2004210/100‑1 (MQ=255)
                                 gTTTGATTGTAACTGGTGCGCCGCTGTGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACg                                                                <  2:1105499/99‑1 (MQ=255)
                                    tGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAc                                                            <  1:2574176/100‑1 (MQ=255)
                                     gATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAcc                                                           >  1:408759/1‑100 (MQ=255)
                                     gATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAcc                                                           >  2:553560/1‑100 (MQ=255)
                                     gATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAcc                                                           >  1:510356/1‑100 (MQ=255)
                                       ttGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTc                                                         >  2:756824/1‑100 (MQ=255)
                                       ttGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTc                                                         <  2:763347/100‑1 (MQ=255)
                                        tGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCg                                                        <  1:1691012/100‑1 (MQ=255)
                                        tGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCg                                                        >  1:1837448/1‑100 (MQ=255)
                                         gTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGa                                                       >  2:757404/1‑100 (MQ=255)
                                         gTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGa                                                       >  1:518480/1‑100 (MQ=255)
                                           aaCTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACg                                                     <  1:2354173/100‑1 (MQ=255)
                                              tGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTg                                                  <  2:1507366/100‑1 (MQ=255)
                                                  gcgcCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTg                                              >  2:2336998/1‑100 (MQ=255)
                                                  gcgcCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTg                                              <  1:1532848/100‑1 (MQ=255)
                                                   cgccgcTGGGCCTGGTGGAGTTTCCTGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTg                                              <  2:902309/99‑1 (MQ=255)
                                                       gcTGGGCCTGTTGGAGTTTAAGGATGTCGCTTTCTGGCGGCAATTCAAACAGGTGCTGGAGTGGTCGAACGATCACGTCACCGCGACGCTGTTTGTct                                           <  2:257869/98‑1 (MQ=255)
                                                       gcTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTctgc                                         >  1:17091/1‑100 (MQ=255)
                                                       gcTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTctgc                                         <  1:2553611/100‑1 (MQ=255)
                                                        cTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGAGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTct                                           >  2:61053/1‑97 (MQ=255)
                                                          gggCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTgg                                      >  1:734035/1‑100 (MQ=255)
                                                            gCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGc                                    <  2:496112/100‑1 (MQ=255)
                                                             ccTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCg                                   <  1:960632/100‑1 (MQ=255)
                                                              cTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGc                                    <  1:2377197/98‑1 (MQ=255)
                                                              cTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCgg                                  <  1:1849894/100‑1 (MQ=255)
                                                              cTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCgg                                  <  1:2479149/100‑1 (MQ=255)
                                                               tcgtgcAGCTTCATGATATCGCTTACTTGCCGGAATTCAACCATGTTCTGGAGTGGTCGAACGATCACGTCATCTCGACGCTGTGTGTCTGCTGGGc                                    <  2:2320493/95‑1 (MQ=255)
                                                                ggtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTa                                >  2:1328767/1‑100 (MQ=255)
                                                                 gtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTAc                               >  1:2525951/1‑100 (MQ=255)
                                                                   ggAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTACTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAg                             <  2:1824918/100‑1 (MQ=255)
                                                                    gggTTTAATGATGTCGCGTACGGGCCGCAGTTCAAACAGGTGCAGGAGTGGTCGAAAGATGACGTCACCTCGGCGCTGTTTGTctgct                                        <  2:862200/86‑1 (MQ=255)
                                                                     aGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGc                           >  2:527648/1‑100 (MQ=255)
                                                                      gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                          >  1:2039719/1‑100 (MQ=255)
                                                                      gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                          >  1:2309374/1‑100 (MQ=255)
                                                                        ntAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgc                        >  1:1734284/2‑100 (MQ=255)
                                                                         tAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgcg                       >  2:2102452/1‑100 (MQ=255)
                                                                         tAATGATGGCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCGGTTTGTCTGCTGGGCGGTACAGGCcgc                        >  2:917385/1‑99 (MQ=255)
                                                                          aaTGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCgcgc                      <  2:1407227/100‑1 (MQ=255)
                                                                           atgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCt                     <  2:543005/100‑1 (MQ=255)
                                                                            tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTc                    >  1:38578/1‑100 (MQ=255)
                                                                            tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTc                    >  2:2408083/1‑100 (MQ=255)
                                                                            tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTc                    <  1:1145826/100‑1 (MQ=255)
                                                                            tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGATGGGCGGTACAGGCCGCGCTc                    >  1:1126897/1‑100 (MQ=255)
                                                                             gatgTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCa                   <  2:2082049/100‑1 (MQ=255)
                                                                              atgtcggacaCTGGCCGCCGTTCAAACAGGTGCTGGAGTGGTCGACAGATCACGCCCCCTCGACGCTGTTTGTCTGATGGGCGGTACAGGCCGCGCGCac                  <  2:1516663/91‑2 (MQ=255)
                                                                              atgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCaa                  >  1:172413/1‑100 (MQ=255)
                                                                                gTCTCTTACTGGCCGCACTTCAAACAGGTGCTGGAGTTGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAata                <  2:2118160/100‑1 (MQ=255)
                                                                                  cGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATc              <  1:774829/100‑1 (MQ=255)
                                                                                  cGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATc              >  1:574522/1‑100 (MQ=255)
                                                                                  cGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTAGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGAGCGGTACAGGCCGCGCTCAatat               >  1:614483/1‑99 (MQ=255)
                                                                                   gCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCTAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATcc             >  2:1598088/1‑100 (MQ=255)
                                                                                   gCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATcc             <  2:1046013/100‑1 (MQ=255)
                                                                                   gCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATcc             <  1:73403/100‑1 (MQ=255)
                                                                                      tACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCctct          <  1:211276/100‑1 (MQ=255)
                                                                                        cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc        >  1:1125400/1‑100 (MQ=255)
                                                                                        cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc        <  2:64160/100‑1 (MQ=255)
                                                                                        cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc        <  1:1223174/100‑1 (MQ=255)
                                                                                        cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTATAc        >  2:328761/1‑100 (MQ=255)
                                                                                          ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACgg      >  2:697542/1‑100 (MQ=255)
                                                                                          ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACgg      >  2:659030/1‑100 (MQ=255)
                                                                                          ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACgg      >  2:1056005/1‑100 (MQ=255)
                                                                                          ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACgg      >  1:281722/1‑100 (MQ=255)
                                                                                            ccGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCa    <  2:17091/100‑1 (MQ=255)
                                                                                              gCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCAtt  <  1:432894/100‑1 (MQ=255)
                                                                                                  |                                                                                               
ACTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATT  >  NC_000913/4214551‑4214744

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: