Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,646,678 G→A 57.3% 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,6780GA57.3% 37.3 / 90.8 84G127D (GGC→GAC) gorglutathione oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (18/18);  new base A (23/25);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

ATTCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGG  >  NC_000913/3646581‑3646774
                                                                                                 |                                                                                                
aTTCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAacgacg                                                                                                <  2:1978946/100‑1 (MQ=255)
   cATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAGTCAAAGGCTTTGACCGCTTCGTTGATGACCAAACGCTGGAGGTAAacgacg                                                                                                >  2:1876472/1‑97 (MQ=255)
   cATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGaaa                                                                                             >  1:1982036/1‑100 (MQ=255)
   cATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGaaa                                                                                             <  1:1190472/100‑1 (MQ=255)
      aCTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCa                                                                                          >  2:918998/1‑100 (MQ=255)
       cTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCAt                                                                                         <  1:1854484/100‑1 (MQ=255)
        ttCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATc                                                                                        <  1:1045765/100‑1 (MQ=255)
          ccTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAATGGCGAAACCATCAc                                                                                      >  2:617781/1‑100 (MQ=255)
             aTGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAATGCTGGTGGTTATCGGCGAAACCATCACGGc                                                                                   <  2:542655/100‑1 (MQ=255)
             aTGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCCCGGc                                                                                   <  2:1668781/100‑1 (MQ=255)
             aTGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGc                                                                                   >  1:2245937/1‑100 (MQ=255)
              tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACg                                                                                     >  2:1921473/1‑97 (MQ=255)
              tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGcc                                                                                  <  1:1849257/100‑1 (MQ=255)
                aaaaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGa                                                                                >  1:702215/1‑100 (MQ=255)
                 aaaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGAt                                                                               <  1:1756590/100‑1 (MQ=255)
                  aaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                              <  2:282167/100‑1 (MQ=255)
                  aaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                              <  1:1689088/100‑1 (MQ=255)
                      tGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCatat                                                                          <  1:1359056/100‑1 (MQ=255)
                       gCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGCAACCATCACGGCCGATCATAtt                                                                         >  2:319796/1‑100 (MQ=255)
                        cTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTc                                                                        <  1:918999/100‑1 (MQ=255)
                         tCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCt                                                                       >  1:1430674/1‑100 (MQ=255)
                          cGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTg                                                                      <  1:1771885/100‑1 (MQ=255)
                          cGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGTCCGATCATATTCTg                                                                      >  1:886222/1‑100 (MQ=255)
                             tAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATc                                                                   >  2:1722102/1‑100 (MQ=255)
                              aaaaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCg                                                                  >  1:588345/1‑100 (MQ=255)
                                aaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGcc                                                                <  2:582646/100‑1 (MQ=255)
                                aaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGcc                                                                <  2:1679975/100‑1 (MQ=255)
                                aaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCAATCATATTCTGATCGcc                                                                <  2:1126009/100‑1 (MQ=255)
                                  aTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAATCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCcac                                                              <  2:571648/100‑1 (MQ=255)
                                   tAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCACAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTAATCGCcaca                                                             >  1:691688/1‑100 (MQ=255)
                                   tAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCcaca                                                             <  1:1018683/100‑1 (MQ=255)
                                   tAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCcaca                                                             <  2:853030/100‑1 (MQ=255)
                                    aaCGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAg                                                            >  2:1114430/1‑100 (MQ=255)
                                    aaCGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAg                                                            <  2:1256441/100‑1 (MQ=255)
                                       gTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCg                                                         <  2:365419/100‑1 (MQ=255)
                                        ttGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCGAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCgg                                                        >  2:1192951/1‑100 (MQ=255)
                                        ttGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCgg                                                        <  1:1829066/100‑1 (MQ=255)
                                         tGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCg                                                         >  1:756985/1‑98 (MQ=255)
                                         tGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgt                                                       <  1:1287448/100‑1 (MQ=255)
                                          gATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtc                                                      >  1:1174325/1‑100 (MQ=255)
                                          gATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtc                                                      >  1:1383679/1‑100 (MQ=255)
                                           aTGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcg                                                     >  1:259227/1‑100 (MQ=255)
                                            tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgt                                                    >  2:1051542/1‑100 (MQ=255)
                                            tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgt                                                    >  2:1194313/1‑100 (MQ=255)
                                               aaTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCg                                                 <  2:1255255/100‑1 (MQ=255)
                                                aTCAAAGTCTTTGCCCGCTTCGTTGATGCCAAAAAGCAGGAGGTACACGGCGAGACCATCACAACCGATCAGATTATGATCGCCACAAGCGGTCGTCCGa                                                >  2:1855935/1‑100 (MQ=255)
                                                aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                                >  1:1781644/1‑100 (MQ=255)
                                                  cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                              >  1:147005/1‑100 (MQ=255)
                                                   aaaGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGcc                                             >  1:1850959/1‑100 (MQ=255)
                                                         tttGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGACCGAGCCACCCgg                                       <  2:1383673/100‑1 (MQ=255)
                                                           tGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGat                                     >  2:353947/1‑100 (MQ=255)
                                                            gCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGata                                    >  1:31782/1‑100 (MQ=255)
                                                              ccGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATAtt                                  >  1:49282/1‑100 (MQ=255)
                                                               cGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTc                                 >  1:1168041/1‑100 (MQ=255)
                                                               cGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTc                                 <  1:73595/100‑1 (MQ=255)
                                                                gCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGaatatcc                                <  2:247646/100‑7 (MQ=255)
                                                                  ttCGTTGATGTCCAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCAACAGGCGTTCGTCTGAGCCACCCGGATATTCCgg                              <  2:937952/100‑1 (MQ=255)
                                                                  ttCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCgg                              <  1:421891/100‑1 (MQ=255)
                                                                   tCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCggg                             <  2:134994/100‑1 (MQ=255)
                                                                   tCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCggg                             <  1:1362645/100‑1 (MQ=255)
                                                                    cGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGc                            >  2:1108980/1‑100 (MQ=255)
                                                                     gTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCg                           <  1:639596/100‑1 (MQ=255)
                                                                        gATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTgg                        >  2:747540/1‑100 (MQ=255)
                                                                        gATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTgg                        <  2:1659393/100‑1 (MQ=255)
                                                                          tGCCAAATCTTTGGAGGTAATTTGCGTATTCTTCTCGTCTGATCATATTCTGATCGCCACAGGCGGTCGTCCGGGCCACCCGGATATTCCGGGCGTgg                        <  1:115859/98‑1 (MQ=255)
                                                                           gCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAAt                     <  1:441037/100‑1 (MQ=255)
                                                                           gCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCAACCCGGATATTCCGGGCGTGGAAt                     >  2:211204/1‑100 (MQ=255)
                                                                           gCCAAAACGCTGGAGGTAAAAGGAGAAACCAGCACGGCCGATCAGATTCTGATCGCCACAGGCGGGCCTCCGAGCCACCCGGATATTTCGGGCGAGGa                       >  2:1940999/1‑98 (MQ=255)
                                                                            ccAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                    >  1:1321523/1‑100 (MQ=255)
                                                                            ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                    >  1:1424778/1‑100 (MQ=255)
                                                                            ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                    >  1:942793/1‑100 (MQ=255)
                                                                             cAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCATCCGAGCCACCCGGATATTCCGGGCGTGGACTAc                   >  1:303095/1‑100 (MQ=255)
                                                                               aaaCGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACgg                 <  1:200266/100‑1 (MQ=255)
                                                                                aaCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGt                <  1:2035984/100‑1 (MQ=255)
                                                                                    cTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTg            <  1:1135652/100‑1 (MQ=255)
                                                                                      ggAGGTAAACGGTGAAACCATCACTGACGTTCATATTCTGATCGCCACAGGCTGTCGTCCGAGCCACCTGGATATTCCGGGCGTGGAATACGGTATTGAt          <  2:486820/100‑1 (MQ=255)
                                                                                      ggAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCAGGCGTGGAATACGGTATTGAt          >  2:29307/1‑100 (MQ=255)
                                                                                        aGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGTCACTCGGATATTCCGGGCGTGGAATACGGTATTGAtt         <  2:1113314/99‑1 (MQ=255)
                                                                                        aGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTc        <  1:211204/100‑1 (MQ=255)
                                                                                        aGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTc        <  1:1997848/100‑1 (MQ=255)
                                                                                         ggTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCt       >  1:928854/1‑100 (MQ=255)
                                                                                           tAAACTACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGAt          <  1:604571/95‑1 (MQ=255)
                                                                                           tAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGAATTGATTc        >  1:1728335/1‑97 (MQ=255)
                                                                                           tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa     >  1:2034709/1‑100 (MQ=255)
                                                                                            aaaCGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGAt    <  2:702215/100‑1 (MQ=255)
                                                                                            aaaCGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGAt    >  2:1047436/1‑100 (MQ=255)
                                                                                              acgacgAAACCATGACGGCCGAGCATAGTGTGATCGCCACAGGCGGTTGTGCGAGACACCCGGATATTACGGGCGTGGAATACGGTCTTGATGCTGATgg  <  2:1320088/100‑1 (MQ=255)
                                                                                              acgacgAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATgg  >  2:633781/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ATTCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGG  >  NC_000913/3646581‑3646774

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