Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,646,678 G→A 100% 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,6780GA89.6% 347.1 / 25.5 134G127D (GGC→GAC) gorglutathione oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (5/9);  new base A (69/50);  total (74/60)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 135 total aligned reads displayed.

TTCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATG  >  NC_000913/3646582‑3646773
                                                                                                |                                                                                               
ttCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAAcgacga                                                                                              <  2:915210/100‑1 (MQ=255)
 tCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTAGAGGTAAACGGCGaa                                                                                             >  2:2868778/1‑100 (MQ=255)
   aTACTTCCTATGAAAACGTGCTTGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAAc                                                                                           <  2:397295/100‑1 (MQ=255)
    tACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAAcc                                                                                          <  1:1971366/100‑1 (MQ=255)
          cTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACg                                                                                    >  2:350008/1‑100 (MQ=255)
          cTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACg                                                                                    <  2:168733/100‑1 (MQ=255)
           tATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACgg                                                                                   >  1:1017900/1‑100 (MQ=255)
           tATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACgg                                                                                   >  1:1374933/1‑100 (MQ=255)
           tATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACgg                                                                                   >  2:1360312/1‑100 (MQ=255)
            aTGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGc                                                                                  <  2:1858740/100‑1 (MQ=255)
            aTGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGc                                                                                  >  2:14355/1‑100 (MQ=255)
             tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGcc                                                                                 <  1:2129757/100‑1 (MQ=255)
             tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGcc                                                                                 >  1:1017344/1‑100 (MQ=255)
             tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGcc                                                                                 >  2:1389547/1‑100 (MQ=255)
             tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGcc                                                                                 >  1:1975403/1‑100 (MQ=255)
             tGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGcc                                                                                 >  1:1017305/1‑100 (MQ=255)
               acaaCGTGCTCGGAAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGa                                                                               >  2:1963933/3‑100 (MQ=255)
               aaaaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCg                                                                                >  1:862938/1‑99 (MQ=255)
                 aaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                             <  1:2105007/100‑1 (MQ=255)
                 aaCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                             <  2:2788077/100‑1 (MQ=255)
                 aaCGTGCTCGGTAAAAATAACGTTGACGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATc                                                                             <  1:2344127/100‑1 (MQ=255)
                  aCGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCa                                                                            <  2:1133431/100‑1 (MQ=255)
                    gTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCata                                                                          >  1:1666547/1‑100 (MQ=255)
                     tGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCatat                                                                         <  2:2059729/100‑1 (MQ=255)
                       cTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTc                                                                       <  2:2740162/100‑1 (MQ=255)
                       cTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTc                                                                       >  2:2187791/1‑100 (MQ=255)
                        tCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCt                                                                      >  1:1280863/1‑100 (MQ=255)
                          ggTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGa                                                                    >  2:611416/1‑100 (MQ=255)
                           gTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGAt                                                                   >  2:2897300/1‑100 (MQ=255)
                            tAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGAt                                                                   >  1:2446959/1‑99 (MQ=255)
                            tAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATc                                                                  >  1:1834192/1‑100 (MQ=255)
                             aaaaaTAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCg                                                                 <  2:170586/100‑1 (MQ=255)
                                  tAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCcaca                                                            >  2:2325094/1‑100 (MQ=255)
                                   ccgatctaTGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAg                                                           <  1:2237275/93‑1 (MQ=255)
                                     cGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGc                                                         >  2:2084835/1‑100 (MQ=255)
                                     cGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGc                                                         >  2:605230/1‑100 (MQ=255)
                                      gTTGATGTAATCAAAGGCTTTGGCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGATACCATCACGGCCGATCATATTTTGATCGCCACAg                                                           <  1:1392620/97‑1 (MQ=255)
                                       ttGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCTAAAACGCTGGAGGTAAACGACGAAACCATAACGGCCTATCATATTCTGATCGCcaca                                                            >  1:342345/1‑95 (MQ=255)
                                        tGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgt                                                      <  1:1350606/100‑1 (MQ=255)
                                           tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgt                                                   >  2:461702/1‑100 (MQ=255)
                                           tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgt                                                   >  2:601041/1‑100 (MQ=255)
                                           tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtcgt                                                   >  1:2228482/1‑100 (MQ=255)
                                           tGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGgtc                                                     >  2:2791294/1‑98 (MQ=255)
                                             tAATCAAAGGCTTTGCCCGCGTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAGCCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTcc                                                 <  1:2313884/100‑1 (MQ=255)
                                              aaTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCg                                                >  2:2807479/1‑100 (MQ=255)
                                               aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCg                                                >  2:772856/1‑99 (MQ=255)
                                               aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                               >  1:43359/1‑100 (MQ=255)
                                               aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                               >  1:44128/1‑100 (MQ=255)
                                               aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                               >  1:998831/1‑100 (MQ=255)
                                               aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                               <  2:2686646/100‑1 (MQ=255)
                                               aTCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGa                                               >  1:802037/1‑100 (MQ=255)
                                                tCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAg                                              <  1:928787/100‑1 (MQ=255)
                                                 cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             <  2:1630856/100‑1 (MQ=255)
                                                 cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGACACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             >  1:714577/1‑100 (MQ=255)
                                                 cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             <  1:867617/100‑1 (MQ=255)
                                                 cAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGc                                             >  1:152359/1‑100 (MQ=255)
                                                          tGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGat                                    >  1:2187430/1‑100 (MQ=255)
                                                          tGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGat                                    >  1:2847324/1‑100 (MQ=255)
                                                              cGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTc                                >  1:2393314/1‑100 (MQ=255)
                                                                cTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCg                              <  1:371267/100‑1 (MQ=255)
                                                                cTTCGTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCg                              >  2:455801/1‑100 (MQ=255)
                                                                  tcattGCTGCCACAACGCTGGAGTTAAACGTCGAAACCATCACGTCTGGTCATATTCTGATCTCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCgg                             <  1:1901861/96‑1 (MQ=255)
                                                                    gTTGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCg                          <  1:264659/100‑1 (MQ=255)
                                                                     ttGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCg                          <  2:1280868/99‑1 (MQ=255)
                                                                      tGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTg                        <  2:2516019/100‑1 (MQ=255)
                                                                      tGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTg                        <  2:2680594/100‑1 (MQ=255)
                                                                      tGATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTg                        >  1:1338748/1‑100 (MQ=255)
                                                                       ggtGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATGTTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTgg                       <  1:272439/98‑1 (MQ=255)
                                                                       gATGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTgg                       >  1:2720368/1‑100 (MQ=255)
                                                                        aTGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGa                      <  2:395329/100‑1 (MQ=255)
                                                                         tGCCAAAAGGCTGGAGGTAACCGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCCGATATTCCGGGCGTGGaa                     <  1:699330/100‑1 (MQ=255)
                                                                         tGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCCTAGTCTGATCGCCACAGGCGGTCGGCCGAGCCACCCGGCTATTCCGGGTGTgg                       >  1:1652607/1‑98 (MQ=255)
                                                                         tGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGaa                     >  2:1401875/1‑100 (MQ=255)
                                                                         tGCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGaa                     >  2:881088/1‑100 (MQ=255)
                                                                          gCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAAt                    <  2:1503572/100‑1 (MQ=255)
                                                                          gCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGaa                     <  2:1315516/99‑1 (MQ=255)
                                                                          gCCAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAAt                    >  2:472064/1‑100 (MQ=255)
                                                                           ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                   >  2:2095155/1‑100 (MQ=255)
                                                                           ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                   >  1:649432/1‑100 (MQ=255)
                                                                           ccAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATa                   >  2:695293/1‑100 (MQ=255)
                                                                            naaaaCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATAc                  >  1:915087/2‑100 (MQ=255)
                                                                            cAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATAc                  >  1:1639049/1‑100 (MQ=255)
                                                                            cAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATAc                  <  2:1720753/100‑1 (MQ=255)
                                                                            cAAAACGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATAc                  >  2:89215/1‑100 (MQ=255)
                                                                              aaaCGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACgg                <  1:908175/100‑1 (MQ=255)
                                                                               aaCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGt               <  2:2187443/100‑1 (MQ=255)
                                                                                aCGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTa              <  2:2187699/100‑1 (MQ=255)
                                                                                 cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAt             <  1:1195935/100‑1 (MQ=255)
                                                                                 cGCTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGCCCGAGCCACCCGGATATTCCGGGCTTGGAATACGGTAt             <  1:1797617/100‑1 (MQ=255)
                                                                                  gCTGGATGTAAACTACGAAACCATCTCAGCCGAGCATATTCTGATCGCCACAGACGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTAtt            <  1:973226/100‑1 (MQ=255)
                                                                                   cTGGAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTg           <  2:1541665/100‑1 (MQ=255)
                                                                                     ggAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGAt         <  1:1651383/100‑1 (MQ=255)
                                                                                     ggAGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGAt         >  2:2859393/1‑100 (MQ=255)
                                                                                       aGGTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTc       <  2:2272063/100‑1 (MQ=255)
                                                                                        ggTAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCt      <  1:723711/100‑1 (MQ=255)
                                                                                          tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa    >  2:2793975/1‑100 (MQ=255)
                                                                                          tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa    >  2:2793955/1‑100 (MQ=255)
                                                                                          tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa    <  1:2793957/100‑1 (MQ=255)
                                                                                          tAAACGACGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGa    <  1:2793937/100‑1 (MQ=255)
                                                                                            aaCGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATg  >  1:1543904/1‑100 (MQ=255)
                                                                                                |                                                                                               
TTCATACTTCCTATGAAAACGTGCTCGGTAAAAATAACGTTGATGTAATCAAAGGCTTTGCCCGCTTCGTTGATGCCAAAACGCTGGAGGTAAACGGCGAAACCATCACGGCCGATCATATTCTGATCGCCACAGGCGGTCGTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATG  >  NC_000913/3646582‑3646773

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: