Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,646,678 G→A G127D (GGC→GAC)  gor → glutathione oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,646,6780GA100.0% 362.5 / NA 100G127D (GGC→GAC) gorglutathione oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (49/51);  total (49/51)

CATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATG  >  NC_000913/3646584‑3646773
                                                                                              |                                                                                               
cATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGaaa                                                                                            <  1:1457817/100‑1 (MQ=255)
cATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGaaa                                                                                            >  1:1500231/1‑100 (MQ=255)
 aTACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAAc                                                                                           >  2:966935/1‑100 (MQ=255)
 aTACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAAc                                                                                           <  1:777732/100‑1 (MQ=255)
  tACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGTGGTAAACGACGAAAcc                                                                                          >  2:1612698/1‑100 (MQ=255)
  tACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAAcc                                                                                          >  2:1071948/1‑100 (MQ=255)
  tACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAAcc                                                                                          <  1:775906/100‑1 (MQ=255)
   aCTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCa                                                                                         <  2:1222890/100‑1 (MQ=255)
   aCTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCa                                                                                         >  1:1586191/1‑100 (MQ=255)
       ccTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCAc                                                                                     >  1:1770258/1‑100 (MQ=255)
        cTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACg                                                                                    <  1:1315030/100‑1 (MQ=255)
        cTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACg                                                                                    >  2:808430/1‑100 (MQ=255)
        cTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACg                                                                                    >  1:2652786/1‑100 (MQ=255)
         tATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACgg                                                                                   >  1:661934/1‑100 (MQ=255)
         tATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACgg                                                                                   >  2:642741/1‑100 (MQ=255)
         tATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACgg                                                                                   >  1:247054/1‑100 (MQ=255)
          aTGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGc                                                                                  <  1:2289725/100‑1 (MQ=255)
            gAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCg                                                                                <  2:618950/100‑1 (MQ=255)
             aaaaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGa                                                                               >  2:1340212/1‑100 (MQ=255)
             aaaaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGa                                                                               >  2:1074408/1‑100 (MQ=255)
              aaaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGAt                                                                              <  1:2633757/100‑1 (MQ=255)
               aaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                             <  2:606179/100‑1 (MQ=255)
               aaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCACGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                             <  2:327785/100‑1 (MQ=255)
                aCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCa                                                                            <  1:937637/100‑1 (MQ=255)
                aCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCa                                                                            <  1:479442/100‑1 (MQ=255)
                aCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCa                                                                            <  1:1331694/100‑1 (MQ=255)
                aCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCa                                                                            <  2:2298592/100‑1 (MQ=255)
                   tGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCatat                                                                         <  2:729802/100‑1 (MQ=255)
                     cTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTc                                                                       <  1:725767/100‑1 (MQ=255)
                     cTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTc                                                                       <  1:808436/100‑1 (MQ=255)
                      tCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCt                                                                      <  1:265350/100‑1 (MQ=255)
                        ggTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGa                                                                    >  2:1522702/1‑100 (MQ=255)
                         gTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGAt                                                                   >  1:2120111/1‑100 (MQ=255)
                         gTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGAt                                                                   >  1:1329851/1‑100 (MQ=255)
                           aaaaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCg                                                                 <  1:334040/100‑1 (MQ=255)
                           aaaaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCg                                                                 <  1:54642/100‑1 (MQ=255)
                           aaaaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCg                                                                 <  1:966941/100‑1 (MQ=255)
                            aaaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGc                                                                >  1:2605748/1‑100 (MQ=255)
                               aTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCcac                                                             <  1:773374/100‑1 (MQ=255)
                               aTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCcac                                                             >  2:606999/1‑100 (MQ=255)
                               aTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCcac                                                             <  2:620076/100‑1 (MQ=255)
                                tAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCcaca                                                            <  2:1329843/100‑1 (MQ=255)
                                       gATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtc                                                     >  2:1858759/1‑100 (MQ=255)
                                       gATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtc                                                     >  2:762030/1‑100 (MQ=255)
                                         tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgt                                                   <  2:997539/100‑1 (MQ=255)
                                          gTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgtc                                                  <  1:500358/100‑1 (MQ=255)
                                          gTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgtc                                                  >  1:1552042/1‑100 (MQ=255)
                                          gTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgtc                                                  <  2:2397888/100‑1 (MQ=255)
                                          gTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgtc                                                  <  2:1500220/100‑1 (MQ=255)
                                           tAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTcc                                                 >  2:1490304/1‑100 (MQ=255)
                                            aaTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCg                                                <  1:1040875/100‑1 (MQ=255)
                                             aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                               <  2:901271/100‑1 (MQ=255)
                                               cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             >  2:203360/1‑100 (MQ=255)
                                               cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             <  1:681474/100‑1 (MQ=255)
                                               cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             <  1:1915009/100‑1 (MQ=255)
                                                aaaGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGcc                                            >  2:260856/1‑100 (MQ=255)
                                                aaaGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGcc                                            <  1:145781/100‑1 (MQ=255)
                                                    gCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCAccc                                        <  1:1441210/100‑1 (MQ=255)
                                                     cTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCg                                       >  1:2597233/1‑100 (MQ=255)
                                                     cTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCg                                       >  1:1321475/1‑100 (MQ=255)
                                                      tttGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCgg                                      >  2:2147489/1‑100 (MQ=255)
                                                      tttGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCgg                                      >  1:729627/1‑100 (MQ=255)
                                                       ttGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGa                                     >  1:1520108/1‑100 (MQ=255)
                                                        tGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGat                                    <  2:2652766/100‑1 (MQ=255)
                                                         gCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGata                                   >  1:2308808/1‑100 (MQ=255)
                                                            cGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTc                                <  1:1851013/100‑1 (MQ=255)
                                                             gCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTcc                               >  1:1663390/1‑100 (MQ=255)
                                                             gCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTcc                               >  2:1945790/1‑100 (MQ=255)
                                                               ttCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCgg                             <  1:1426417/100‑1 (MQ=255)
                                                                tCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCgg                             >  2:1893284/1‑99 (MQ=255)
                                                                 cGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGc                           <  1:1612711/100‑1 (MQ=255)
                                                                   ttGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGt                         >  1:682289/1‑100 (MQ=255)
                                                                    tGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTg                        <  2:793401/100‑1 (MQ=255)
                                                                    tGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTg                        <  2:2308790/100‑1 (MQ=255)
                                                                     gATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACATGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTgg                       <  1:1945806/100‑1 (MQ=255)
                                                                     gATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTgg                       <  1:1893299/100‑1 (MQ=255)
                                                                      aTGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGa                      >  2:2506493/1‑100 (MQ=255)
                                                                      aTGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGa                      >  1:1757595/1‑100 (MQ=255)
                                                                       tGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGaa                     >  1:598471/1‑100 (MQ=255)
                                                                       tGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGaa                     >  1:2052374/1‑100 (MQ=255)
                                                                       tGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGaa                     >  2:30224/1‑100 (MQ=255)
                                                                         ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                   >  1:942454/1‑100 (MQ=255)
                                                                         ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                   >  2:756140/1‑100 (MQ=255)
                                                                           aaaaCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACg                 <  2:2052357/100‑1 (MQ=255)
                                                                             aaCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGt               >  2:503838/1‑100 (MQ=255)
                                                                              aCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTa              <  1:1749703/100‑1 (MQ=255)
                                                                               cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAt             <  1:1071955/100‑1 (MQ=255)
                                                                               cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAt             <  1:1340221/100‑1 (MQ=255)
                                                                               cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAt             <  1:2208130/100‑1 (MQ=255)
                                                                               cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAt             <  1:2339541/100‑1 (MQ=255)
                                                                               cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAt             <  1:2599080/100‑1 (MQ=255)
                                                                                gCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAtt            <  1:203363/100‑1 (MQ=255)
                                                                                  tGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGa          <  1:2445869/100‑1 (MQ=255)
                                                                                   ggAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGAt         >  2:1253762/1‑100 (MQ=255)
                                                                                     aGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTc       >  1:1451400/1‑100 (MQ=255)
                                                                                        tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa    >  1:2316647/1‑100 (MQ=255)
                                                                                        tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa    >  2:849906/1‑100 (MQ=255)
                                                                                         aaaCGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGAt   >  2:2492639/1‑100 (MQ=255)
                                                                                         aaaCGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGAt   >  2:2180592/1‑100 (MQ=255)
                                                                                          aaCGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATg  >  2:2183303/1‑100 (MQ=255)
                                                                                              |                                                                                               
CATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATG  >  NC_000913/3646584‑3646773

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