Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 138,888 T→C T780A (ACG→GCG)  gcd ← glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913138,8880TC86.7% 733.9 / 94.5 302T780A (ACG→GCG) gcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (16/24);  new base C (133/128);  total (149/153)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 303 total aligned reads displayed.

CATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTGGT  >  NC_000913/138796‑138979
                                                                                            |                                                                                           
      caacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCaaa                                                                                       <  2:5221655/93‑1 (MQ=255)
        acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATg                                                                                     <  1:5359526/93‑1 (MQ=255)
         cgcgCACGTTATCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGa                                                                                    <  2:5913760/93‑1 (MQ=255)
          gcgcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGaa                                                                                   <  2:6705318/93‑1 (MQ=255)
           cgcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAc                                                                                  <  2:4761629/93‑1 (MQ=255)
           cgcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAAc                                                                                  <  1:7268609/93‑1 (MQ=255)
            gcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAAcc                                                                                 <  2:1448900/93‑1 (MQ=255)
            gcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAAcc                                                                                 <  1:6106767/93‑1 (MQ=255)
            gcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAAcc                                                                                 <  2:1737369/93‑1 (MQ=255)
                gTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGa                                                                             >  1:751165/1‑93 (MQ=255)
                gTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGa                                                                             >  1:5750923/1‑93 (MQ=255)
                 ttGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGAc                                                                            >  2:6014548/1‑93 (MQ=255)
                 ttGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGAc                                                                            >  1:905882/1‑93 (MQ=255)
                     cAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCCCCAAATGAACCGTGACCGcc                                                                        <  1:1158183/93‑1 (MQ=255)
                      aGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCt                                                                       >  1:6995395/1‑93 (MQ=255)
                      aGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCt                                                                       <  1:4812263/93‑1 (MQ=255)
                      aGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCt                                                                       >  2:7304434/1‑93 (MQ=255)
                      aGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCt                                                                       >  1:1288939/1‑93 (MQ=255)
                        cAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGc                                                                     >  2:6208037/1‑93 (MQ=255)
                         aaTCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCg                                                                    >  1:3854932/1‑93 (MQ=255)
                             tGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGgaga                                                                <  1:5982862/93‑1 (MQ=255)
                             tGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGgaga                                                                <  2:4449001/93‑1 (MQ=255)
                             tGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGgaga                                                                >  2:5693058/1‑93 (MQ=255)
                              gAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGAt                                                               <  1:5837636/93‑1 (MQ=255)
                              gAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCCCAAATGAACCGTGACCGCCTGCGGAGATc                                                              <  2:458799/93‑1 (MQ=255)
                               aGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATc                                                              <  2:2230450/93‑1 (MQ=255)
                                gCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATca                                                             <  2:3922798/93‑1 (MQ=255)
                                gCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATca                                                             <  1:6764302/93‑1 (MQ=255)
                                 cAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcac                                                            <  2:1849987/93‑1 (MQ=255)
                                   aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                          >  2:2560207/1‑93 (MQ=255)
                                   aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                          >  1:3210870/1‑93 (MQ=255)
                                   aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                          >  2:6990131/1‑93 (MQ=255)
                                   aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                          <  1:3433619/93‑1 (MQ=255)
                                    gTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcgccac                                                         <  2:508058/93‑1 (MQ=255)
                                    gTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcaccac                                                         >  2:3705946/1‑93 (MQ=255)
                                     tCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACcaca                                                        >  1:7321656/1‑93 (MQ=255)
                                     tCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACcaca                                                        >  2:4089628/1‑93 (MQ=255)
                                      cTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACAt                                                       >  2:6051851/1‑93 (MQ=255)
                                       ttACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATa                                                      >  1:5948932/1‑93 (MQ=255)
                                       ttACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATa                                                      >  1:4998730/1‑93 (MQ=255)
                                          cTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTg                                                   >  2:2870960/1‑93 (MQ=255)
                                           ttCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGCTCACCACATACTGc                                                  >  2:7207692/1‑93 (MQ=255)
                                           ttCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc                                                  >  2:4564791/1‑93 (MQ=255)
                                           ttCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc                                                  >  1:4263068/1‑93 (MQ=255)
                                             cacaTCACCCGGCAGCGCATAAGCCACAATATAGTCCCCCATCTTCGCACCAAATGAACCGTGCCCGCCTGCGGAGATCACCACATACTGCtt                                                <  1:6606031/93‑1 (MQ=255)
                                               catcatCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTAc                                              >  1:1582701/1‑93 (MQ=255)
                                                atcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTAcc                                             >  2:607653/1‑93 (MQ=255)
                                                atcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTAcc                                             >  2:3039606/1‑93 (MQ=255)
                                                 tcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCa                                            <  2:826084/93‑1 (MQ=255)
                                                 tcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCa                                            <  1:1911908/93‑1 (MQ=255)
                                                 tcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCa                                            >  1:3365612/1‑93 (MQ=255)
                                                  catcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCAt                                           <  1:1036978/93‑1 (MQ=255)
                                                   atcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCAtt                                          >  2:1024339/1‑93 (MQ=255)
                                                    tcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTc                                         <  2:1279092/93‑1 (MQ=255)
                                                    tcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTc                                         >  1:2284675/1‑93 (MQ=255)
                                                    tcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTc                                         >  1:6406639/1‑93 (MQ=255)
                                                     ccGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCa                                        <  1:5239313/93‑1 (MQ=255)
                                                     ccGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCa                                        <  2:4969107/93‑1 (MQ=255)
                                                      cGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCAc                                       >  1:311620/1‑93 (MQ=255)
                                                      cGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCAc                                       >  1:1359801/1‑93 (MQ=255)
                                                       ggCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACt                                      <  2:5589456/93‑1 (MQ=255)
                                                        gCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACtt                                     >  2:7509934/1‑93 (MQ=255)
                                                        gCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACtt                                     >  1:2172442/1‑93 (MQ=255)
                                                         cAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTc                                    >  2:6455351/1‑93 (MQ=255)
                                                          aGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCa                                   <  1:1422771/93‑1 (MQ=255)
                                                           gcgcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCAt                                  >  2:1884210/1‑93 (MQ=255)
                                                           gcgcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCAt                                  <  2:3808408/93‑1 (MQ=255)
                                                            cgcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATa                                 <  2:3127971/93‑1 (MQ=255)
                                                              cATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAgg                               >  2:4414069/1‑93 (MQ=255)
                                                              cATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAgg                               <  2:5456021/93‑1 (MQ=255)
                                                              cATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAgg                               >  1:3558296/1‑93 (MQ=255)
                                                              cATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAgg                               >  2:425844/1‑93 (MQ=255)
                                                               aTAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGt                              >  2:3479066/1‑93 (MQ=255)
                                                                tAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTc                             <  1:5077709/93‑1 (MQ=255)
                                                                tAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTc                             <  2:7022159/93‑1 (MQ=255)
                                                                 aaGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCa                            <  2:3338696/93‑1 (MQ=255)
                                                                  aGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCAt                           >  1:1740523/1‑93 (MQ=255)
                                                                   gCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCAtt                          <  2:6391358/93‑1 (MQ=255)
                                                                    ccACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTg                         <  2:6883107/93‑1 (MQ=255)
                                                                    ccACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTg                         >  2:2937273/1‑93 (MQ=255)
                                                                    ccACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTg                         <  1:7449047/93‑1 (MQ=255)
                                                                     cacaATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTgg                        >  2:1199590/1‑93 (MQ=255)
                                                                     cacaATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTgg                        >  1:2734325/1‑93 (MQ=255)
                                                                       caATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCg                      >  1:4416289/1‑93 (MQ=255)
                                                                         atatAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTa                    >  2:467483/1‑93 (MQ=255)
                                                                          tataGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAg                   >  2:5305475/1‑93 (MQ=255)
                                                                           ataGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGc                  >  2:5537008/1‑93 (MQ=255)
                                                                             aGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCt                >  1:1807320/1‑93 (MQ=255)
                                                                             aGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCt                <  1:1107501/93‑1 (MQ=255)
                                                                             aGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCt                <  1:6920228/93‑1 (MQ=255)
                                                                               tCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGa              >  1:412682/1‑93 (MQ=255)
                                                                                cGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGac             <  2:5812274/93‑1 (MQ=255)
                                                                                 gCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGacc            >  2:1972313/1‑93 (MQ=255)
                                                                                    cATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGaccacc         >  2:1645392/1‑93 (MQ=255)
                                                                                     aTCTTCGTACCAAATGAACCGTGACCGTCTGCGGAGATCACCACATACTGCTTAGCATTCACTTCATAGGTCATTGGCGTAGCCTGACCAccc        >  1:5662244/1‑93 (MQ=255)
                                                                                     aTCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCAccc        <  2:1123254/93‑1 (MQ=255)
                                                                                       cTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGc      <  1:2459764/93‑1 (MQ=255)
                                                                                         tCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTg    <  1:5210163/93‑1 (MQ=255)
                                                                                         tCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTg    <  2:1307321/93‑1 (MQ=255)
                                                                                           gTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTGGt  >  2:5977314/1‑93 (MQ=255)
                                                                                            |                                                                                           
CATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTGGT  >  NC_000913/138796‑138979

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