Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,369,793 G→A 20.9% 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,7930GA20.9% 221.9 / 48.5 105E35E (GAG→GAA) pspEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base G (46/37);  new base A (14/8);  total (60/45)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

TAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAA  >  NC_000913/1369694‑1369891
                                                                                                   |                                                                                                  
tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAagag                                                                                                    >  2:3276090/1‑100 (MQ=255)
tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAagag                                                                                                    <  1:2118966/100‑1 (MQ=255)
tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTATAGCCGCTGAACACTGGATCGATGTTCGTGTGCCAGAGCAGTATCAGCAagag                                                                                                    >  2:11218/1‑100 (MQ=255)
 aaaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGc                                                                                                   <  2:1547293/100‑1 (MQ=255)
  aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCa                                                                                                  >  1:547028/1‑100 (MQ=255)
  aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCa                                                                                                  <  2:1893456/100‑1 (MQ=255)
       ggCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc                                                                                             <  1:104460/100‑1 (MQ=255)
       ggCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTc                                                                                             >  1:110877/1‑100 (MQ=255)
       ggCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTc                                                                                             >  1:382431/1‑100 (MQ=255)
        gCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCa                                                                                            <  1:3188800/100‑1 (MQ=255)
        gCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCa                                                                                            <  1:2432003/100‑1 (MQ=255)
         cTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAg                                                                                           >  1:2621827/1‑100 (MQ=255)
            aCTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAgggg                                                                                        <  1:3335088/100‑1 (MQ=255)
             cTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAggggg                                                                                       <  2:3001220/100‑1 (MQ=255)
             cTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAggggg                                                                                       >  1:2954549/1‑100 (MQ=255)
             cTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAggggg                                                                                       <  1:241232/100‑1 (MQ=255)
                 ctctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCAt                                                                                   >  2:1455603/1‑100 (MQ=255)
                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATc                                                                                  <  1:1001475/100‑1 (MQ=255)
                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATc                                                                                  <  1:199262/100‑1 (MQ=255)
                   ctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGCACACTGGATCGATGTTCGTTTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATc                                                                                  <  2:2250170/99‑1 (MQ=255)
                   ctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCa                                                                                 >  1:3507773/1‑100 (MQ=255)
                    tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCaa                                                                                <  2:1626629/100‑1 (MQ=255)
                       cACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAatat                                                                             <  1:794406/100‑1 (MQ=255)
                       cACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAatat                                                                             <  1:768349/100‑1 (MQ=255)
                         cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTc                                                                           >  1:21722/1‑100 (MQ=255)
                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTcc                                                                          <  2:2551595/100‑1 (MQ=255)
                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTcc                                                                          <  2:2339787/100‑1 (MQ=255)
                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcc                                                                          <  1:2002867/100‑1 (MQ=255)
                           ggTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTccc                                                                         >  2:2852023/1‑100 (MQ=255)
                           ggTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTccc                                                                         >  2:202738/1‑100 (MQ=255)
                           ggTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTcc                                                                          >  1:4110976/1‑99 (MQ=255)
                           ggTCTTTCCACTGCCCGTTTTCGCCGCTGCACACTGGATCGATGTTCGTGCTCCAGAGCAGTATCCGCAAGAACACGTCCAGGGGGCCATCAATATTccc                                                                         >  2:916651/1‑100 (MQ=255)
                            gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTcccc                                                                        <  2:1383507/100‑1 (MQ=255)
                              gtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTg                                                                      <  1:2722250/100‑1 (MQ=255)
                              aacaccagtGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTg                                                                      <  1:663308/92‑1 (MQ=255)
                                tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                    >  1:2342554/1‑100 (MQ=255)
                                tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGCTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                    >  1:2382035/1‑100 (MQ=255)
                                 ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                   <  2:4130781/100‑1 (MQ=255)
                                 ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                   <  1:1603471/100‑1 (MQ=255)
                                     cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                               <  1:2996817/100‑1 (MQ=255)
                                      tGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                               >  2:3062799/1‑99 (MQ=255)
                                      tGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGt                                                              >  2:1122532/1‑100 (MQ=255)
                                       gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                             <  1:1455613/100‑1 (MQ=255)
                                       gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                             >  1:2095423/1‑100 (MQ=255)
                                       gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                             >  1:68586/1‑100 (MQ=255)
                                       gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                             >  1:2725288/1‑100 (MQ=255)
                                       gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg                                                             <  1:3499730/100‑1 (MQ=255)
                                        cccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGa                                                            <  1:303395/100‑1 (MQ=255)
                                        cccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGa                                                            <  1:3294010/100‑1 (MQ=255)
                                                cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                    >  2:1299790/1‑100 (MQ=255)
                                                cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGTGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcg                                                     >  2:1808260/1‑99 (MQ=255)
                                                    gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg                                                >  1:1733080/1‑100 (MQ=255)
                                                    gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg                                                <  1:3575011/100‑1 (MQ=255)
                                                    gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAAGATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                                  >  2:3873779/1‑98 (MQ=255)
                                                    gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg                                                >  1:1666740/1‑100 (MQ=255)
                                                    gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg                                                <  2:3727944/100‑1 (MQ=255)
                                                       gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                             >  2:3566979/1‑100 (MQ=255)
                                                        aaCACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAc                                            >  1:1436752/1‑100 (MQ=255)
                                                         acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc                                           <  2:744820/100‑1 (MQ=255)
                                                            ctcgATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc                                        >  2:3379547/4‑100 (MQ=255)
                                                               gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt                                     >  2:1868345/1‑100 (MQ=255)
                                                                 tCGATGTTCGTGTGCCAGAGCAGTATCAGCATGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTcc                                   >  1:2780219/1‑100 (MQ=255)
                                                                  cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg                                  >  2:889652/1‑100 (MQ=255)
                                                                  cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg                                  >  1:2088182/1‑100 (MQ=255)
                                                                   gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                                 <  2:1232140/100‑1 (MQ=255)
                                                                   gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                                 <  1:606189/100‑1 (MQ=255)
                                                                   gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                                 <  1:2038121/100‑1 (MQ=255)
                                                                   gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCgg                                 <  1:3506092/100‑1 (MQ=255)
                                                                    aTGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGa                                <  2:6074/100‑1 (MQ=255)
                                                                     tGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                               >  1:3614733/1‑100 (MQ=255)
                                                                     tGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                               <  1:741781/100‑1 (MQ=255)
                                                                     tGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                               <  2:2805091/100‑1 (MQ=255)
                                                                      gtgcGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCCATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                               <  2:3839493/96‑1 (MQ=255)
                                                                      gTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa                              <  2:3540276/100‑1 (MQ=255)
                                                                      gTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa                              <  1:971765/100‑1 (MQ=255)
                                                                      gTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa                              >  2:935160/1‑100 (MQ=255)
                                                                        tCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaa                            >  1:1118723/1‑100 (MQ=255)
                                                                          gtgtTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaaaa                          >  1:1990794/1‑100 (MQ=255)
                                                                           tgtTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                         <  1:2805995/100‑1 (MQ=255)
                                                                            gtTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                        >  2:128968/1‑100 (MQ=255)
                                                                             ttCCAGAGCAGTATCAGCACGAGCACGTTCAGGGGGCCATCAATATTACCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTACGGATAAAAACGa                       >  2:2901142/1‑100 (MQ=255)
                                                                              tCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac                      >  1:3399916/1‑100 (MQ=255)
                                                                               ccAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                     >  2:445227/1‑100 (MQ=255)
                                                                               ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                     >  1:2369580/1‑100 (MQ=255)
                                                                                cAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac                    >  1:123627/1‑100 (MQ=255)
                                                                                cAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGCGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac                    >  2:2893428/1‑100 (MQ=255)
                                                                                 agagCAGTATCAGCAAGAGCACGTTAAGGGGGCCAGCAATATTCCCCTGAAGGAAGTGCAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac                    >  2:2695646/1‑99 (MQ=255)
                                                                                    gCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg                >  1:3155405/1‑100 (MQ=255)
                                                                                     cAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGa               >  1:702522/1‑100 (MQ=255)
                                                                                      aGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGGGAAAGAGCGCATTGTCACCGCCGTTCCGGATGAAAACGACACCGTGaa              >  2:851903/1‑100 (MQ=255)
                                                                                       gTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa             >  1:2932261/1‑100 (MQ=255)
                                                                                       gTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa             >  2:861890/1‑100 (MQ=255)
                                                                                       gTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa             >  1:525986/1‑100 (MQ=255)
                                                                                         aTCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgt           >  2:2961709/1‑100 (MQ=255)
                                                                                           cAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtgt         >  1:2139781/1‑100 (MQ=255)
                                                                                            aGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtgt         >  2:1526779/1‑99 (MQ=255)
                                                                                            aGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGAGTa        >  2:3535512/1‑100 (MQ=255)
                                                                                               aaGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGGCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTg     >  1:2902162/1‑100 (MQ=255)
                                                                                                agagCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc    >  1:2478102/1‑100 (MQ=255)
                                                                                                  agCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCa   >  2:1808980/1‑99 (MQ=255)
                                                                                                   |                                                                                                  
TAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAA  >  NC_000913/1369694‑1369891

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