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.2% 579.7 / ‑1.2 163Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (66/94);  total (67/96)
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.94e-01

Only 100 of 163 total aligned reads displayed.

GTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCG  >  NC_000913/125291‑125481
                                                                                                |                                                                                              
gTGACGTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACTCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAgg                                                                                             <  1:4831385/100‑1 (MQ=255)
    cctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgat                                                                                         <  1:515148/100‑1 (MQ=255)
       cctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                      >  1:469099/1‑100 (MQ=255)
          tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                   <  1:4187115/100‑1 (MQ=255)
          tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                   >  1:835401/1‑100 (MQ=255)
          tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                   >  1:2378645/1‑100 (MQ=255)
           cACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAgg                                                                                             >  2:3783198/1‑89 (MQ=255)
           cACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                  <  1:454714/100‑1 (MQ=255)
             ccGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtg                                                                                           >  2:3898022/1‑89 (MQ=255)
             ccGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGc                                                                                >  1:1484614/1‑100 (MQ=255)
             ccGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGc                                                                                >  1:1490560/1‑100 (MQ=255)
              cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                          <  2:2741626/89‑1 (MQ=255)
               gAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTc                                                                              >  1:2681989/1‑100 (MQ=255)
                aGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                        <  2:1825791/89‑1 (MQ=255)
                aGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                        <  2:3708951/89‑1 (MQ=255)
                aGCTGGCGCGTGATGGTCAGGAATGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                        <  2:3296676/89‑1 (MQ=255)
                   tGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGc                                                                          <  1:4367035/100‑1 (MQ=255)
                     gcgcGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAg                                                                        <  1:3103058/100‑1 (MQ=255)
                      cgcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGt                                                                       <  1:4725847/100‑1 (MQ=255)
                          tGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTc                                                                              >  2:1608741/1‑89 (MQ=255)
                           gATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTcc                                                                             >  2:2740916/1‑89 (MQ=255)
                             tGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCt                                                                >  1:1723024/1‑100 (MQ=255)
                              ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGc                                                                          <  2:2669419/89‑1 (MQ=255)
                              ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTa                                                               <  1:4428450/100‑1 (MQ=255)
                              ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                          >  2:4099164/1‑89 (MQ=255)
                                 cAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCg                                                            >  1:910872/1‑100 (MQ=255)
                                   ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                          <  1:3332629/100‑1 (MQ=255)
                                   ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                          >  1:3512834/1‑100 (MQ=255)
                                    gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGc                                                                    <  2:2256650/89‑1 (MQ=255)
                                     aTTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACta                                                        <  1:3843565/100‑1 (MQ=255)
                                      ttGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAt                                                                  >  2:4686081/1‑89 (MQ=255)
                                       tgtgAACTCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATc                                                                 >  2:1679129/1‑89 (MQ=255)
                                       tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                      >  1:2203376/1‑100 (MQ=255)
                                       tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                      >  1:1290193/1‑100 (MQ=255)
                                        gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                     >  1:3823355/1‑100 (MQ=255)
                                        gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                     <  1:2509410/100‑1 (MQ=255)
                                        gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                     <  1:1293841/100‑1 (MQ=255)
                                          gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                              >  2:1378871/1‑89 (MQ=255)
                                          gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                   >  1:2082235/1‑100 (MQ=255)
                                          gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                   >  1:3607653/1‑100 (MQ=255)
                                           aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaa                                                  <  1:2443200/100‑1 (MQ=255)
                                             cGCTGGAACATGCTGCACCCGCTGGATACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCTTCAGTGGCATCTACCGACTATATGAAAc                                                <  1:19765/100‑1 (MQ=255)
                                              gCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                          <  2:4834631/89‑1 (MQ=255)
                                              gCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                          <  2:3825463/89‑1 (MQ=255)
                                               cTGGAACATGCTGCACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCAACt                                                         <  2:1546528/89‑1 (MQ=255)
                                               cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACt                                                         <  2:4047896/89‑1 (MQ=255)
                                                 ggAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtat                                                       >  2:3116601/1‑89 (MQ=255)
                                                 ggAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGtt                                            <  1:1577215/100‑1 (MQ=255)
                                                    aCATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGc                                         <  1:4623599/100‑1 (MQ=255)
                                                    aCATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGc                                         <  1:2652712/100‑1 (MQ=255)
                                                     cATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                   <  2:1917662/89‑1 (MQ=255)
                                                     cATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCt                                        <  1:2222214/100‑1 (MQ=255)
                                                      aTGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTg                                       <  1:1011957/100‑1 (MQ=255)
                                                      aTGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTg                                       <  1:4168899/100‑1 (MQ=255)
                                                       tgctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGa                                      >  1:3702134/1‑100 (MQ=255)
                                                        gctgctCTCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCCTCTACCGACTATATGAAAc                                                >  2:765877/1‑89 (MQ=255)
                                                        gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAc                                                >  2:1386825/1‑89 (MQ=255)
                                                        gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAg                                     >  1:631136/1‑100 (MQ=255)
                                                          tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCa                                   >  1:4719627/1‑100 (MQ=255)
                                                          ngcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCa                                   >  1:570757/2‑100 (MQ=255)
                                                             atccGCTGGAAACTCCTCGCGTTCCGTATATCACTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                           <  2:2511460/87‑1 (MQ=255)
                                                             aCCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                           <  2:2511562/89‑1 (MQ=255)
                                                              cccGCTGGGAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                          <  2:4832799/89‑1 (MQ=255)
                                                              cccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                          >  2:3029592/1‑89 (MQ=255)
                                                              cccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                          >  2:785733/1‑89 (MQ=255)
                                                              ccagCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGCTGAACGACGCGCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                          >  2:3556146/4‑89 (MQ=255)
                                                               ccGCTGGAAACTCCGCGCGTTCCGTATCTCGCTTAGGTGATGAACGACGCTCCGCCAGTGCCCTCTACCGCCTATATGAAACTGTTCGc                                         >  2:3672509/1‑89 (MQ=255)
                                                                cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                             <  1:1895657/100‑1 (MQ=255)
                                                                 gCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                            >  1:4107681/1‑100 (MQ=255)
                                                                 gCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                            <  1:2740916/100‑1 (MQ=255)
                                                                 gCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGAAGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTg                                       <  2:2207288/89‑1 (MQ=255)
                                                                  cTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGa                                      <  2:3977650/89‑1 (MQ=255)
                                                                   tGGAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCCTATATGAAACTGTTCGCTGAg                                     <  2:1168909/89‑1 (MQ=255)
                                                                    ggAAACTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                    <  2:3823355/89‑1 (MQ=255)
                                                                    ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                    <  2:3856207/89‑1 (MQ=255)
                                                                    ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                    >  2:958973/1‑89 (MQ=255)
                                                                    ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACt                         <  1:4035472/100‑1 (MQ=255)
                                                                      aaaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                       <  1:1477068/100‑1 (MQ=255)
                                                                       aaCTCCGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAgg                                 >  2:4433426/1‑89 (MQ=255)
                                                                       aaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTAc                      <  1:1022398/100‑1 (MQ=255)
                                                                        aCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACg                     <  1:3880589/100‑1 (MQ=255)
                                                                         cTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTc                               <  2:2217135/89‑1 (MQ=255)
                                                                         cTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTc                               <  2:4460248/89‑1 (MQ=255)
                                                                           ccGCGCGTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                  >  1:3131145/1‑100 (MQ=255)
                                                                           ccGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                  >  1:1710505/1‑100 (MQ=255)
                                                                            cgcgcgTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                            <  2:1291665/89‑1 (MQ=255)
                                                                            cgcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                            >  2:4157658/1‑89 (MQ=255)
                                                                             gcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTa                           >  2:4439202/1‑89 (MQ=255)
                                                                             gcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTa                           <  2:581924/89‑1 (MQ=255)
                                                                                 gTTCCGTATATCGCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                       >  2:4473490/1‑89 (MQ=255)
                                                                                  ntCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTga           >  1:3723955/2‑100 (MQ=255)
                                                                                     cGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTa                   <  2:3393006/89‑1 (MQ=255)
                                                                                        atatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTa     >  1:687871/1‑100 (MQ=255)
                                                                                         tatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTAc    >  1:4061316/1‑100 (MQ=255)
                                                                                          atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc              >  2:3754770/1‑89 (MQ=255)
                                                                                          atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc              <  2:1436881/89‑1 (MQ=255)
                                                                                          atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc              >  2:1120029/1‑89 (MQ=255)
                                                                                          atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGGACCGGCTGACGACTAcc   <  1:3518740/100‑1 (MQ=255)
                                                                                           tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg  <  1:937434/100‑1 (MQ=255)
                                                                                           tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg  >  1:4463193/1‑100 (MQ=255)
                                                                                                |                                                                                              
GTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCG  >  NC_000913/125291‑125481

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: