Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,369,793 G→A 100% E35E (GAG→GAA pspE → thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,369,7930GA98.0% 305.6 / ‑2.5 99E35E (GAG→GAApspEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (47/50);  total (48/51)
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.04e-01

AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTG  >  NC_000913/1369696‑1369888
                                                                                                 |                                                                                               
aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACa                                                                                               <  1:1590959/100‑1 (MQ=255)
         tACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCATGTTCCAGAGCAGTATCAGCAAGAACACGTTCAggg                                                                                      >  1:624568/1‑100 (MQ=255)
            ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGc                                                                                   <  2:1681601/100‑1 (MQ=255)
            ttGCTCTGGCACTGGTGTTTTCAATGCCCGTTTTCGCAGCTGAACACTTGATCGATGTTGGTGTTCCAGAGCAGTATCAGTAAGAACACGTTCAGGGGGc                                                                                   <  2:599847/100‑1 (MQ=255)
              gCTCTGGCACTGGTGTTTTTTTTGTTTTTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTGCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGcc                                                                                  >  2:965225/1‑99 (MQ=255)
              gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCa                                                                                 <  1:1517653/100‑1 (MQ=255)
               ctctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCAt                                                                                <  2:2867714/100‑1 (MQ=255)
               ctctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACAAGTTCAGGGGGCCAt                                                                                >  2:2876703/1‑100 (MQ=255)
                tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATc                                                                               <  1:1389233/100‑1 (MQ=255)
                 ctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCa                                                                              <  2:2375259/100‑1 (MQ=255)
                  tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCaa                                                                             >  1:1447953/1‑100 (MQ=255)
                   ggCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAat                                                                            >  1:2116164/1‑100 (MQ=255)
                    gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGATCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAata                                                                           <  2:511578/100‑1 (MQ=255)
                     cACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAatat                                                                          <  1:1750644/100‑1 (MQ=255)
                      aCTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCATTATCAGCAAGAACACTTTCAGGGGGCCATCAATAtt                                                                         <  2:245336/100‑1 (MQ=255)
                         ggTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTccc                                                                      >  2:2190164/1‑100 (MQ=255)
                          gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCCGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc                                                                     <  2:1647209/100‑1 (MQ=255)
                          gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc                                                                     <  2:654663/100‑1 (MQ=255)
                            gtTTTCACTGCCCGTTTTCGCCGCTGAACACAGGATCGATGTTCGTGTTCAAGAACAGTATCAGCAAGAACCCGTTCAGGGGGCCAAAAATATTcccc                                                                     >  2:1038971/1‑98 (MQ=255)
                               ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                <  1:1297387/100‑1 (MQ=255)
                               ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                <  1:824928/100‑1 (MQ=255)
                               ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                <  2:997005/100‑1 (MQ=255)
                                  aCTGCCCGTTTTCGCCGCTGAACACTGGATAGATGTTCGTGTTCCAGAGCAGTATCAACAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaagga                                                              >  2:630298/1‑99 (MQ=255)
                                   cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                            >  2:1882966/1‑100 (MQ=255)
                                   cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                            >  1:2224331/1‑100 (MQ=255)
                                    tGCCAGTTTTCGACGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATAAGCAAGAACACGTTAAGGGCGCCATTAATATTCCCATgaagaaa                                                             >  2:356916/1‑98 (MQ=255)
                                     gcgcGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                          <  2:2148055/97‑1 (MQ=255)
                                     gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                            >  2:1812839/1‑98 (MQ=255)
                                     gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                          <  2:2206955/100‑1 (MQ=255)
                                      cccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGa                                                         <  1:436780/100‑1 (MQ=255)
                                         gTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAag                                                      <  2:2441554/100‑1 (MQ=255)
                                           tttCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAagag                                                    >  1:1701295/1‑100 (MQ=255)
                                            ttCGCCGCTGAACACTGGTGAGATGTTCGTGTTCCAGAGTAGTATCAGTAAGAACACGTTCAGGGGGCCATCAATTTTTCCCTGAAGGAAGTGAAAGAgc                                                   <  2:55622/100‑1 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                 <  1:2949508/100‑1 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                 >  2:1196746/1‑100 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTACAGAGCAGTATCAGCAAGAACACGTTCAGGGGCCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                 >  2:1546306/1‑100 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                <  2:711218/100‑1 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                >  1:1729566/1‑100 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTACAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                >  2:1885556/1‑100 (MQ=255)
                                                ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                               <  1:2634092/100‑1 (MQ=255)
                                                ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTACAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                               >  1:869309/1‑100 (MQ=255)
                                                 cgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAtt                                              <  1:2990206/100‑1 (MQ=255)
                                                 cgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                               <  1:3025265/99‑1 (MQ=255)
                                                 cgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTAAAGGAAGACAACAATATTCAACTGAAGGAAGTGAAAGAGCGAAtt                                              >  1:848789/1‑100 (MQ=255)
                                                  gcTGAACACTGGATCGGTGTTGGTGTTCCAGAGTAGTATCAGCAAGAACACGTTCAGGGTGCCATCAATATTCCTCTGAAGGAAGTGAAAGAGCGCATTg                                             <  2:127037/100‑1 (MQ=255)
                                                  gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg                                             >  1:130830/1‑100 (MQ=255)
                                                  gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg                                             <  1:915827/100‑1 (MQ=255)
                                                   cTGAACACTGGATCGATGTTCGTGTTCCAGATCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGc                                            <  1:1555997/100‑1 (MQ=255)
                                                     gAACACTGGCTCCATATTAGTGTTAAACAGCAGTATCAGCAACAACAGGTTGAGGGTCCAATCACTATTCCCCTCAAGGAAGTGAAagaca                                                   >  2:2884658/1‑89 (MQ=255)
                                                       acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATCGCCa                                          >  2:1579913/1‑98 (MQ=255)
                                                       acacTGGATCGATGTTCGTGTACCAGAGCAGTATCAGCAAGAACACGTTCAGGGTGCCATCAATATTCCCATGAAGGAAGTAAAAGAGCGCATTCCCAcc                                        >  2:1446353/1‑100 (MQ=255)
                                                         acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGTGGGCCATCAATATTCCCCTGAAGGAAGTTAAAGAGCGCATTGCCAccgc                                      <  2:1267697/100‑1 (MQ=255)
                                                         acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgc                                      >  1:2240347/1‑100 (MQ=255)
                                                         acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgc                                      <  2:2219383/100‑1 (MQ=255)
                                                           tGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGTGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgccg                                    <  2:2866185/100‑1 (MQ=255)
                                                           tGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgccg                                    >  2:1069175/1‑100 (MQ=255)
                                                            ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt                                   <  2:130829/100‑1 (MQ=255)
                                                            ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATAGCCAccgccg                                    >  2:3006336/1‑99 (MQ=255)
                                                              aTCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTc                                 >  1:2148155/1‑100 (MQ=255)
                                                                cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTc                                 >  2:877542/1‑98 (MQ=255)
                                                                 gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                              <  2:2224308/100‑1 (MQ=255)
                                                                 gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                              <  2:2337826/100‑1 (MQ=255)
                                                                 gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAAGAAATGAAAGAGCGCATTGCCACCGCCGTTCCgg                              >  1:1911991/1‑100 (MQ=255)
                                                                  aTGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGa                             <  1:2085302/100‑1 (MQ=255)
                                                                   tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                            <  1:1882988/100‑1 (MQ=255)
                                                                     ttCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaa                          <  1:520316/100‑1 (MQ=255)
                                                                      tCGTGTTCCAGAGCAGTATCAGCAAGAACACGTACAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGCTCCGGATaaa                         >  2:326974/1‑100 (MQ=255)
                                                                       cGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaaa                        >  1:2003788/1‑100 (MQ=255)
                                                                       cGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaaa                        >  2:1137417/1‑100 (MQ=255)
                                                                         tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCTATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                      <  1:1347815/100‑1 (MQ=255)
                                                                         tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTTAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                      >  1:419193/1‑100 (MQ=255)
                                                                         tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                      >  1:1102321/1‑100 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                     <  1:832820/100‑1 (MQ=255)
                                                                            tCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac                   <  2:3029371/100‑1 (MQ=255)
                                                                             ccAGAGCAGTATCAGCAAGAACACGTTCAGAGGGACATCAATATTCCCATGAAGGAAATAAAAGAGCGCATTGCCACCGCCGTTCCGGACAAAAAAGaca                  >  2:2482729/1‑100 (MQ=255)
                                                                                gagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg               <  1:1717953/100‑1 (MQ=255)
                                                                                gagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg               <  1:339023/100‑1 (MQ=255)
                                                                                gagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg               >  2:4637/1‑100 (MQ=255)
                                                                                gagCAGTATCAGAAAGAATACGTTCAGGGGGCCATCAATATATCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg               <  1:2944411/100‑1 (MQ=255)
                                                                                    aGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa           >  1:2854882/1‑100 (MQ=255)
                                                                                    aGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa           <  1:143351/100‑1 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg         <  2:2404993/100‑1 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg         <  1:2488166/100‑1 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg         <  2:624564/100‑1 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCATGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg         <  1:2190187/100‑1 (MQ=255)
                                                                                       ntCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgt        >  1:83248/2‑100 (MQ=255)
                                                                                       aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCACCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa          >  2:1154885/1‑98 (MQ=255)
                                                                                        tCAGCAAGAACACGTTCAGGGGGCCATCAATCTTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCAGTTCCGGATAAAAACGACACCGTGAAAgt        >  1:727224/1‑99 (MQ=255)
                                                                                        tCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtg       <  1:1685055/100‑1 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa     >  1:1459695/1‑100 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa     >  2:630081/1‑100 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa     >  1:2574800/1‑100 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGAGCACCACCAATAAACCCCTGAAGGAAGCGAAAGAGCGCATTGCAACCGCAGTTCCGGATAAAAACGACAACGTGAAAGAGTa     >  2:1433419/1‑100 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt    >  1:1391847/1‑100 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt    >  1:1360916/1‑100 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt    <  2:1205631/100‑1 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCAGGAAGGAAGTGAAAGAGCGCATTGCCACCGCGGTTCCGGATAAAAACGACACCGTGAAAGTGTa     >  2:1199245/1‑99 (MQ=255)
                                                                                            cAAGAACACGTTCAGGGGGCCATCAATATTACCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTACGGATAAAAACGACACCGTGAAAGTGTAtt   >  2:1831096/1‑100 (MQ=255)
                                                                                             aaGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTg  <  1:1541978/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTG  >  NC_000913/1369696‑1369888

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