Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,387 C→T Q791* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,3870CT98.9% 628.2 / ‑2.7 176Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (92/82);  total (94/82)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 176 total aligned reads displayed.

GACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCG  >  NC_000913/125293‑125483
                                                                                              |                                                                                                
gACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtg                                                                                             <  1:1210155/100‑1 (MQ=255)
     cctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                        >  1:2424543/1‑100 (MQ=255)
     cctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                        >  1:2710244/1‑100 (MQ=255)
      ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                       <  1:3612155/100‑1 (MQ=255)
      ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                       >  1:2475198/1‑100 (MQ=255)
      ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                       >  1:1670480/1‑100 (MQ=255)
        tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAg                                                                                                >  2:4005748/1‑89 (MQ=255)
        tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                     <  1:3930404/100‑1 (MQ=255)
         cACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                    >  1:1063232/1‑100 (MQ=255)
         cACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                    <  1:2010054/100‑1 (MQ=255)
          aCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGt                                                                                              <  2:2194107/89‑1 (MQ=255)
          aCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgacg                                                                                   >  1:3980731/1‑100 (MQ=255)
           ccGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGc                                                                                  <  1:722906/100‑1 (MQ=255)
            cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                            <  2:710277/89‑1 (MQ=255)
              aGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                          <  2:2804627/89‑1 (MQ=255)
              aGCTGGCGCGTGATCGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                          <  2:2101939/89‑1 (MQ=255)
                cTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                             <  1:4657912/100‑1 (MQ=255)
                 tGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                            >  1:4407840/1‑100 (MQ=255)
                 tGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                            <  1:5027268/100‑1 (MQ=255)
                  ggCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCa                                                                           >  1:4989317/1‑100 (MQ=255)
                  ggCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCa                                                                           >  1:4847479/1‑100 (MQ=255)
                    cgcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                    <  2:171591/89‑1 (MQ=255)
                    cgcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                    >  2:2913372/1‑89 (MQ=255)
                    cgcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGt                                                                         >  1:1427488/1‑100 (MQ=255)
                      cgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                       <  1:1161975/100‑1 (MQ=255)
                       gTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGc                                                                      <  1:1016014/100‑1 (MQ=255)
                        tGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCa                                                                     <  1:3413298/100‑1 (MQ=255)
                         gATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTcc                                                                               <  2:871263/89‑1 (MQ=255)
                           tGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCt                                                                  >  1:4682072/1‑100 (MQ=255)
                            ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                            <  2:4670360/89‑1 (MQ=255)
                            ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                            <  2:1670480/89‑1 (MQ=255)
                            ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                            <  2:3790866/89‑1 (MQ=255)
                            ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTa                                                                 >  1:223864/1‑100 (MQ=255)
                             gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                                <  1:4056170/100‑1 (MQ=255)
                              tCAGGATTGTGAACGCTGGCACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAcc                                                               >  1:2968953/1‑100 (MQ=255)
                              tCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAg                                                                          >  2:2112043/1‑89 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                       >  2:1617074/1‑89 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                       >  2:1795380/1‑89 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                       <  2:1325837/89‑1 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                       >  2:2740328/1‑89 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                       <  2:4592312/89‑1 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                            <  1:96596/100‑1 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                            >  1:3014347/1‑100 (MQ=255)
                                 ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                            >  1:2382780/1‑100 (MQ=255)
                                   aTTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACta                                                          <  1:2364069/100‑1 (MQ=255)
                                     tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                        <  1:4352068/100‑1 (MQ=255)
                                      gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCt                                                                  >  2:4104798/1‑89 (MQ=255)
                                      gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                       >  1:4147054/1‑100 (MQ=255)
                                      gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                       <  1:621210/100‑1 (MQ=255)
                                       tgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTa                                                                 >  2:3977610/1‑89 (MQ=255)
                                       tgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATg                                                      >  1:1484619/1‑100 (MQ=255)
                                        gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                                >  2:3579710/1‑89 (MQ=255)
                                        gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                     >  1:1647755/1‑100 (MQ=255)
                                         aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaa                                                    <  1:972332/100‑1 (MQ=255)
                                          aCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCGGGCAGTGGCATCTAcca                                                              <  2:2412291/89‑2 (MQ=255)
                                           cGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAc                                                  >  1:1733325/1‑100 (MQ=255)
                                           cGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAc                                                  <  1:3252405/100‑1 (MQ=255)
                                           cGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAc                                                  >  1:255082/1‑100 (MQ=255)
                                             cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTg                                                >  1:521250/1‑100 (MQ=255)
                                             cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTg                                                >  1:2143006/1‑100 (MQ=255)
                                               ggAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtat                                                         >  2:4815874/1‑89 (MQ=255)
                                                gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                        <  2:1880541/89‑1 (MQ=255)
                                                 aaCATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                       <  2:1989236/89‑1 (MQ=255)
                                                  aCATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGc                                           >  1:9520/1‑100 (MQ=255)
                                                    aTGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTACATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaa                                                    <  2:1235342/89‑1 (MQ=255)
                                                      gctgcACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGCGATAAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAc                                                  <  2:4796245/89‑1 (MQ=255)
                                                      gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAc                                                  >  2:2660007/1‑89 (MQ=255)
                                                        tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTg                                                <  2:2615271/89‑1 (MQ=255)
                                                        tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCa                                     >  1:4217251/1‑100 (MQ=255)
                                                         gcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGt                                               >  2:532765/1‑89 (MQ=255)
                                                           aCCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                             >  2:1553989/1‑89 (MQ=255)
                                                            cccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTc                                 >  1:1830428/1‑100 (MQ=255)
                                                             ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGc                                           <  2:3936845/89‑1 (MQ=255)
                                                              cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                               >  1:922233/1‑100 (MQ=255)
                                                                  ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                      <  2:2143006/89‑1 (MQ=255)
                                                                  ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACt                           >  1:3367794/1‑100 (MQ=255)
                                                                  ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACt                           >  1:4631447/1‑100 (MQ=255)
                                                                  ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACt                           <  1:4097608/100‑1 (MQ=255)
                                                                  ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACt                           >  1:2604159/1‑100 (MQ=255)
                                                                   gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACtt                          <  1:2050680/100‑1 (MQ=255)
                                                                      aCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACg                       <  1:610520/100‑1 (MQ=255)
                                                                         ccGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                               >  2:3758275/1‑89 (MQ=255)
                                                                         ccGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                               >  2:3077165/1‑89 (MQ=255)
                                                                          cgcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCGGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc                   >  1:3685015/1‑100 (MQ=255)
                                                                          cgcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc                   >  1:4181430/1‑100 (MQ=255)
                                                                               gTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                         <  2:4914961/89‑1 (MQ=255)
                                                                                  ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGt                      >  2:2440829/1‑89 (MQ=255)
                                                                                  ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacg           >  1:1356794/1‑100 (MQ=255)
                                                                                   cGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTa                     <  2:4670653/89‑1 (MQ=255)
                                                                                    gTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                    <  2:2258301/89‑1 (MQ=255)
                                                                                      atatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTa       >  1:4240709/1‑100 (MQ=255)
                                                                                       tatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCgg                 >  2:1366641/1‑89 (MQ=255)
                                                                                        atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc                <  2:3685015/89‑1 (MQ=255)
                                                                                        atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTAcc     <  1:481029/100‑1 (MQ=255)
                                                                                         tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTagat    >  1:4599358/1‑97 (MQ=255)
                                                                                         tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg    <  1:3108592/100‑1 (MQ=255)
                                                                                         tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg    <  1:1370338/100‑1 (MQ=255)
                                                                                          cGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCTACTATATGAACCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTg              <  2:1832652/89‑1 (MQ=255)
                                                                                           gCTTAGGTGATGAAGGACGGGCGGGCAGTGGCATCTACCGTCTATATGAAAGTGTTCGCTCAGCAGGTCCGTACTTCCGTACCGGCTga             >  2:4954317/1‑89 (MQ=255)
                                                                                           gCTTAGGTGATGAACGACGCTCCGGCAGTGGCAGCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTga             >  2:3858412/1‑89 (MQ=255)
                                                                                              |                                                                                                
GACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCG  >  NC_000913/125293‑125483

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: