Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,649 A→T 61.0% 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,6490AT61.0% 127.5 / 176.8 146I124F (ATC→TTC) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (46/43);  total (75/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 147 total aligned reads displayed.

ACTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATTCCTA  >  NC_000913/4214551‑4214748
                                                                                                  |                                                                                                   
 cTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATc                                                                                                   <  2:1699034/100‑1 (MQ=255)
  tttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCa                                                                                                  <  2:2851568/100‑1 (MQ=255)
  tttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCa                                                                                                  >  2:2969918/1‑100 (MQ=255)
   ttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCaa                                                                                                 <  2:722535/100‑1 (MQ=255)
    tGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCaaa                                                                                                >  1:1799086/1‑100 (MQ=255)
     gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAAc                                                                                               >  2:354665/1‑100 (MQ=255)
     gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAAc                                                                                               <  2:1930607/100‑1 (MQ=255)
     gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAAc                                                                                               <  2:1720711/100‑1 (MQ=255)
     gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAAc                                                                                               >  2:395096/1‑100 (MQ=255)
        tacatTTAGGATCAGAACTTTGACGGTTTGATTTGATCGTTTGCTCCTCTGGGCCGGGTGGAGTTTAATGATGTCGCTTACTTGCCGCAGTTCAAAc                                                                                               <  1:3117546/94‑1 (MQ=255)
        gATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAgg                                                                                            <  2:1269069/100‑1 (MQ=255)
         atatTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGt                                                                                           <  1:361340/100‑1 (MQ=255)
         atatTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGt                                                                                           <  1:2687949/100‑1 (MQ=255)
           atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                         >  1:2898376/1‑100 (MQ=255)
            ttCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCt                                                                                        >  1:3091072/1‑100 (MQ=255)
            ttCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCt                                                                                        >  2:1471499/1‑100 (MQ=255)
            ttCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCt                                                                                        >  2:11082/1‑100 (MQ=255)
              cAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTgg                                                                                      <  2:3085486/100‑1 (MQ=255)
                ggATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAg                                                                                    <  1:1275243/100‑1 (MQ=255)
                 gATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGt                                                                                   >  2:2183419/1‑100 (MQ=255)
                   tCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTgg                                                                                 <  2:245966/100‑1 (MQ=255)
                   tCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTgg                                                                                 >  1:2094697/1‑100 (MQ=255)
                    cAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGt                                                                                <  2:2563773/100‑1 (MQ=255)
                     aGAACTTTGACGGTTTGATTGTAACTGGTGTGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGCTGCTGGCGTGGTc                                                                               <  1:2137380/100‑1 (MQ=255)
                     aGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTc                                                                               <  1:1274471/100‑1 (MQ=255)
                     aGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTc                                                                               <  2:908952/100‑1 (MQ=255)
                     aGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTc                                                                               >  1:744580/1‑100 (MQ=255)
                       aaCTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGa                                                                             <  1:3317984/100‑1 (MQ=255)
                        aCTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGaa                                                                            >  2:1314975/1‑100 (MQ=255)
                         cTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGaaa                                                                           >  2:3335677/1‑100 (MQ=255)
                          tttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAg                                                                          <  1:2145261/100‑1 (MQ=255)
                           ttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGa                                                                         >  1:3364583/1‑100 (MQ=255)
                           ttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGa                                                                         >  1:1488409/1‑100 (MQ=255)
                             gACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATc                                                                       >  2:1202232/1‑100 (MQ=255)
                             gACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATc                                                                       <  1:1581301/100‑1 (MQ=255)
                                ggTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACg                                                                    <  1:1035028/100‑1 (MQ=255)
                                ggTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACg                                                                    <  1:2969895/100‑1 (MQ=255)
                                ggTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACg                                                                    <  2:2477005/100‑1 (MQ=255)
                                   ttGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCa                                                                 >  2:2480293/1‑100 (MQ=255)
                                      aTTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCt                                                              <  2:2950885/100‑1 (MQ=255)
                                       ttGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTc                                                             <  2:183876/100‑1 (MQ=255)
                                       ttGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTc                                                             >  1:2557081/1‑100 (MQ=255)
                                       ttGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTc                                                             >  1:2367690/1‑100 (MQ=255)
                                        tGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCg                                                            >  1:722160/1‑100 (MQ=255)
                                        tGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCg                                                            <  1:2480274/100‑1 (MQ=255)
                                         gTAACTGGTGCGCCGCTTGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGa                                                           >  1:1164739/1‑100 (MQ=255)
                                         gTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGa                                                           <  1:1398389/100‑1 (MQ=255)
                                          tAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGAc                                                          >  1:146217/1‑100 (MQ=255)
                                          tAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGAc                                                          >  2:1906748/1‑100 (MQ=255)
                                           aaCTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACg                                                         <  2:628702/100‑1 (MQ=255)
                                               ggTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGt                                                     <  1:2698271/100‑1 (MQ=255)
                                               ggTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCAAGTCACCTCGACGCTGt                                                     >  1:2939667/1‑100 (MQ=255)
                                               ggTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCAAGTCACCTCGACGCTGt                                                     >  1:2069733/1‑100 (MQ=255)
                                               ggTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGt                                                     <  1:395093/100‑1 (MQ=255)
                                                gTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGtt                                                    <  2:1488423/100‑1 (MQ=255)
                                                gTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGtt                                                    <  2:2094714/100‑1 (MQ=255)
                                                 tGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGttt                                                   >  2:125294/1‑100 (MQ=255)
                                                  gcgcCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTg                                                  <  1:3090904/100‑1 (MQ=255)
                                                    gccgcTGGGCCTGGTGGAGTTTAATGATGTCGCGTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGCAAGATCACGTCACCTCGACGCTGTTTGTc                                                >  1:2080504/1‑100 (MQ=255)
                                                        cTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTctgct                                            >  1:1375627/1‑100 (MQ=255)
                                                         tGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTctgctg                                           <  1:354662/100‑1 (MQ=255)
                                                          gggCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTgg                                          >  2:269285/1‑100 (MQ=255)
                                                            gCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGc                                        <  2:1239103/100‑1 (MQ=255)
                                                            gCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGc                                        <  2:118016/100‑1 (MQ=255)
                                                             ccTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCg                                       <  2:1030373/100‑1 (MQ=255)
                                                                 gtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTAc                                   <  1:890983/100‑1 (MQ=255)
                                                                  tggAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACa                                  >  2:441277/1‑100 (MQ=255)
                                                                   ggAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAg                                 >  2:3100388/1‑100 (MQ=255)
                                                                   ggAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAg                                 <  1:1652733/100‑1 (MQ=255)
                                                                     aGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGc                               <  2:1799102/100‑1 (MQ=255)
                                                                      gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                              >  2:1999729/1‑100 (MQ=255)
                                                                      gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                              >  2:265206/1‑100 (MQ=255)
                                                                        ttAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgc                            >  2:360576/1‑100 (MQ=255)
                                                                        ttAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgc                            >  1:895323/1‑100 (MQ=255)
                                                                         tAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgcg                           >  2:1860824/1‑100 (MQ=255)
                                                                         tAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgcg                           <  1:2971346/100‑1 (MQ=255)
                                                                          aaTGATGTCGCTGACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCgcgc                          <  2:3161493/100‑1 (MQ=255)
                                                                           atgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAGCTCGACGCTGTTTGTCTGCTGTGCGGTGCAGGCCGTGCt                         >  2:693759/1‑100 (MQ=255)
                                                                           atgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCgcgc                          >  1:325394/1‑99 (MQ=255)
                                                                           atgatgTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCt                         <  2:653047/100‑1 (MQ=255)
                                                                            tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTc                        <  2:851719/100‑1 (MQ=255)
                                                                            tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTc                        >  2:1402182/1‑100 (MQ=255)
                                                                               tgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAat                     >  2:2301266/1‑100 (MQ=255)
                                                                                gTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAata                    <  2:1375639/100‑1 (MQ=255)
                                                                                   gCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATcc                 >  1:257115/1‑100 (MQ=255)
                                                                                   gCTTACTGGCCGCAGATCAAACAGGTGTTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATcc                 <  1:1956906/100‑1 (MQ=255)
                                                                                       aCTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTa             >  1:2860198/1‑100 (MQ=255)
                                                                                        cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc            <  1:914328/100‑1 (MQ=255)
                                                                                        cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc            >  2:614001/1‑100 (MQ=255)
                                                                                          ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACgg          >  1:1185513/1‑100 (MQ=255)
                                                                                           gCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGc         >  2:1810945/1‑100 (MQ=255)
                                                                                            ccGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCa        >  2:3244721/1‑100 (MQ=255)
                                                                                            ccGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCa        <  2:3155942/100‑1 (MQ=255)
                                                                                            ccGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCATCTACGGc         >  1:3187752/1‑99 (MQ=255)
                                                                                             cGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGTCCGCGCTCAATATCCTCTACGGCAt       >  2:587204/1‑100 (MQ=255)
                                                                                             cGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCAt       <  2:3036767/100‑1 (MQ=255)
                                                                                             cGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCAt       >  1:2805368/1‑100 (MQ=255)
                                                                                             cGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCAt       >  1:1935059/1‑100 (MQ=255)
                                                                                              gCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCAtt      <  1:3195747/100‑1 (MQ=255)
                                                                                                 gATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATTCCt   >  1:1857518/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ACTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATTCCTA  >  NC_000913/4214551‑4214748

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