Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,214,649 A→T 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,6490AT97.8% 317.8 / ‑0.8 90I124F (ATC→TTC) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (48/40);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATTC  >  NC_000913/4214552‑4214745
                                                                                                 |                                                                                                
cTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTc                                                                                                <  1:1160178/100‑1 (MQ=255)
 tttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTTGCCGCAGTTCa                                                                                               >  1:2114431/1‑100 (MQ=255)
 tttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCa                                                                                               >  1:2111564/1‑100 (MQ=255)
 tttGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATcg                                                                                               >  1:2112914/1‑99 (MQ=255)
   tGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCaaa                                                                                             >  2:128823/1‑100 (MQ=255)
   tGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCaaa                                                                                             >  2:1883408/1‑100 (MQ=255)
    gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAAc                                                                                            <  1:470260/100‑1 (MQ=255)
    gAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAAc                                                                                            >  2:1118193/1‑100 (MQ=255)
      aGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAg                                                                                          >  2:1217807/1‑100 (MQ=255)
         tatTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTg                                                                                       >  2:143507/1‑100 (MQ=255)
          atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                      >  2:1542080/1‑100 (MQ=255)
          atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                      >  2:306885/1‑100 (MQ=255)
          atTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGc                                                                                      >  2:2699496/1‑100 (MQ=255)
             cAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTgg                                                                                   <  1:1019624/100‑1 (MQ=255)
             cAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTgg                                                                                   >  1:1960193/1‑100 (MQ=255)
                gATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGt                                                                                >  1:2563987/1‑100 (MQ=255)
                 aTCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTg                                                                               >  1:1115231/1‑100 (MQ=255)
                 aTCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTg                                                                               <  2:1795496/100‑1 (MQ=255)
                  tCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTgg                                                                              <  1:2052516/100‑1 (MQ=255)
                   cAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGt                                                                             >  2:2444887/1‑100 (MQ=255)
                   cAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGATGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGt                                                                             >  2:723792/1‑100 (MQ=255)
                    aGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTc                                                                            >  2:249906/1‑100 (MQ=255)
                    aGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTc                                                                            >  2:2613126/1‑100 (MQ=255)
                     gAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCg                                                                           <  1:843109/100‑1 (MQ=255)
                       aCTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGaa                                                                         <  2:839019/100‑1 (MQ=255)
                        cTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGaaa                                                                        >  2:1052008/1‑100 (MQ=255)
                        cTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGaaa                                                                        >  1:1809674/1‑100 (MQ=255)
                         tttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAg                                                                       <  2:507929/100‑1 (MQ=255)
                         tttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAg                                                                       <  1:123330/100‑1 (MQ=255)
                         tttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAg                                                                       <  1:1563422/100‑1 (MQ=255)
                         tttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAg                                                                       <  2:312915/100‑1 (MQ=255)
                          ttGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGa                                                                      <  2:732046/100‑1 (MQ=255)
                           tGACGGTTTGATAGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGAt                                                                     >  1:1815018/1‑100 (MQ=255)
                            gACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATc                                                                    <  2:2473369/100‑1 (MQ=255)
                              cGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCAc                                                                  >  1:559585/1‑100 (MQ=255)
                              cGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCAc                                                                  >  1:2433319/1‑100 (MQ=255)
                               ggTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACg                                                                 <  2:931645/100‑1 (MQ=255)
                                 tttGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTc                                                               <  1:306888/100‑1 (MQ=255)
                                 tttGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTc                                                               <  1:1474297/100‑1 (MQ=255)
                                  ttGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCa                                                              <  2:2084089/100‑1 (MQ=255)
                                  ttGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCa                                                              <  2:1005635/100‑1 (MQ=255)
                                  ttGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCa                                                              <  1:2523920/100‑1 (MQ=255)
                                    gATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAcc                                                            <  2:2147886/100‑1 (MQ=255)
                                    gATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCAcc                                                            <  1:2608978/100‑1 (MQ=255)
                                      ttGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTc                                                          >  2:2506189/1‑100 (MQ=255)
                                          aaCTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACg                                                      >  1:1740398/1‑100 (MQ=255)
                                          aaCTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACg                                                      >  2:1731832/1‑100 (MQ=255)
                                               gTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGtt                                                 <  2:2563968/100‑1 (MQ=255)
                                                 gcgcCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTg                                               <  1:18835/100‑1 (MQ=255)
                                                  cgccgcTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGt                                              >  2:2376577/1‑100 (MQ=255)
                                                   gccgcTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTc                                             >  2:791166/1‑100 (MQ=255)
                                                   gccgcTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTc                                             <  1:1542092/100‑1 (MQ=255)
                                                    ccgcTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTGCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTct                                            <  1:1041466/100‑1 (MQ=255)
                                                           gCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGc                                     <  2:175499/100‑1 (MQ=255)
                                                           gCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGc                                     <  2:1012136/100‑1 (MQ=255)
                                                               ggtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTa                                 >  1:2256378/1‑100 (MQ=255)
                                                               ggtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTa                                 <  2:2256360/100‑1 (MQ=255)
                                                                gtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACATGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTAc                                >  2:2497151/1‑100 (MQ=255)
                                                                gtggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTAc                                <  1:858899/100‑1 (MQ=255)
                                                                 tggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACa                               <  1:1625971/100‑1 (MQ=255)
                                                                  ggAGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAg                              >  2:1159595/1‑100 (MQ=255)
                                                                    aGTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGc                            <  2:1136009/100‑1 (MQ=255)
                                                                     gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                           >  2:616514/1‑100 (MQ=255)
                                                                     gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                           >  2:2445597/1‑100 (MQ=255)
                                                                     gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGcc                           >  1:1323827/1‑100 (MQ=255)
                                                                     gTTTAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGc                            >  1:2248610/1‑99 (MQ=255)
                                                                       ttAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcgc                         <  1:1851617/100‑1 (MQ=255)
                                                                       ttAATGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCcg                          <  2:1849478/99‑1 (MQ=255)
                                                                         aaTGATGTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCgcgc                       <  1:2613145/100‑1 (MQ=255)
                                                                           tgatgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTc                     <  1:2699516/100‑1 (MQ=255)
                                                                             atgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCaa                   <  1:461319/100‑1 (MQ=255)
                                                                             atgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCaa                   <  1:310270/100‑1 (MQ=255)
                                                                             atgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCaa                   <  1:249909/100‑1 (MQ=255)
                                                                             atgTCGCTTACTGGCCGCAGTTCAAACAGGCGCTGGAGTGGTCGAAAGATCACGTCACCTCGAAGCTGTTTGTCTGCTCGGCGGTACAGGCCGCGCTc                     >  2:713370/1‑98 (MQ=255)
                                                                              tgTCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAat                  >  2:525875/1‑100 (MQ=255)
                                                                                tCGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAatat                >  2:2647434/1‑100 (MQ=255)
                                                                                 cGCTTACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATc               >  2:2342696/1‑100 (MQ=255)
                                                                                     tACTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCctct           <  1:128825/100‑1 (MQ=255)
                                                                                      aCTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCTATTTCCTCTa          >  2:1843639/1‑100 (MQ=255)
                                                                                       cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc         >  1:2319385/1‑100 (MQ=255)
                                                                                       cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc         >  1:1331001/1‑100 (MQ=255)
                                                                                       cTGGCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTAc         >  1:483528/1‑100 (MQ=255)
                                                                                         ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCACTATCCTCTACgg       >  1:2541279/1‑100 (MQ=255)
                                                                                         ggCCGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACgg       <  1:2444906/100‑1 (MQ=255)
                                                                                          gCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCAGCGCTCAATATCCTCTACGGc      >  2:1979370/1‑100 (MQ=255)
                                                                                           ccGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCa     >  2:2631782/1‑100 (MQ=255)
                                                                                           ccGCAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCa     >  1:862564/1‑100 (MQ=255)
                                                                                           ccGCAGTACAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCa     >  1:1454917/1‑100 (MQ=255)
                                                                                            cGCAGTTCAAACAGGTGCTGGAGCGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCAt    >  2:2380876/1‑100 (MQ=255)
                                                                                              cAGTTCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATTc  >  1:2332355/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTTTGAAGATATTCAGGATCAGAACTTTGACGGTTTGATTGTAACTGGTGCGCCGCTGGGCCTGGTGGAGTTTAATGATGTCGCTTACTGGCCGCAGATCAAACAGGTGCTGGAGTGGTCGAAAGATCACGTCACCTCGACGCTGTTTGTCTGCTGGGCGGTACAGGCCGCGCTCAATATCCTCTACGGCATTC  >  NC_000913/4214552‑4214745

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: