Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 15.3% intergenic (+28/+790) sseA → / ← sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,3670GA15.3% 286.3 / 37.9 118intergenic (+28/+790)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (44/56);  new base A (8/10);  total (52/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

Only 100 of 125 total aligned reads displayed.

GTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAC  >  NC_000913/2653268‑2653465
                                                                                                   |                                                                                                  
gTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTAAAAAAAAAAAAAAAccgccg                                                                                                    <  1:1636329/100‑1 (MQ=255)
 tGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGa                                                                                                   >  2:2464036/1‑100 (MQ=255)
  gAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGaa                                                                                                  <  2:1662594/100‑1 (MQ=255)
       tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                             >  2:2585607/1‑100 (MQ=255)
       tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                             <  1:2454805/100‑1 (MQ=255)
       tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg                                                                                              <  1:1441806/99‑1 (MQ=255)
        ccgatctggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                            <  1:1136501/93‑1 (MQ=255)
         tACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                           <  2:2505541/100‑1 (MQ=255)
         tACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                           <  2:1597077/100‑1 (MQ=255)
          acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          <  1:1997969/100‑1 (MQ=255)
           cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGt                                                                                         >  2:1060160/1‑100 (MQ=255)
           caacgGCGCATGGGGTGAATGGGGCGCACGAGCAGATTTACCAGTTGAGCCCGTAAAATAAGCTTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGt                                                                                         >  2:1748618/3‑100 (MQ=255)
            gacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtt                                                                                        <  2:1337373/100‑1 (MQ=255)
              cgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtttt                                                                                      <  2:1011893/100‑1 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGcagatcga                                                                                     >  2:1136505/1‑93 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                     >  1:58773/1‑100 (MQ=255)
                gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                    >  2:2385620/1‑100 (MQ=255)
                gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                    <  2:2527173/100‑1 (MQ=255)
                 cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                   <  1:1232110/100‑1 (MQ=255)
                  gcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                  <  2:2248422/100‑1 (MQ=255)
                   cATGGAGTGAATGGGGCGCGCGGGCAGATTTGCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                 <  1:1745038/100‑1 (MQ=255)
                   cATGGAGTGAATGGGGCGCACGGGCAGATTTACCGGTTGAGCCAGGGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAt                                                                                  <  1:2305206/99‑1 (MQ=255)
                    aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                <  2:1439320/100‑1 (MQ=255)
                    aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTg                                                                                <  2:330531/100‑1 (MQ=255)
                    aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTg                                                                                <  2:1319824/100‑1 (MQ=255)
                     tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                               <  2:714330/100‑1 (MQ=255)
                     tGGAGTGAATGGGGCGCGCGGGAAGATTTACCGGTTGAGCGAGTTAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                               <  2:537191/100‑1 (MQ=255)
                      ggAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                              <  1:2050821/100‑1 (MQ=255)
                        aGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAg                                                                            <  1:1839989/100‑1 (MQ=255)
                           gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                         >  2:1935985/1‑100 (MQ=255)
                               ggggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGtt                                                                     >  1:246900/1‑100 (MQ=255)
                                gggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTc                                                                    >  2:2017640/1‑100 (MQ=255)
                                   cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                 <  1:1935975/100‑1 (MQ=255)
                                   cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                 <  2:1979893/100‑1 (MQ=255)
                                      gcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                              <  1:294824/100‑1 (MQ=255)
                                       cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                             >  1:1177765/1‑100 (MQ=255)
                                         ggCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCttt                                                           >  1:653874/1‑100 (MQ=255)
                                         ggCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCttt                                                           >  1:1349836/1‑100 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         >  1:2018446/1‑100 (MQ=255)
                                            aGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                        <  2:1472176/100‑1 (MQ=255)
                                             gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       >  1:1291953/1‑100 (MQ=255)
                                             gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       >  2:767370/1‑100 (MQ=255)
                                              aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       <  1:530066/99‑1 (MQ=255)
                                               tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagca                                                     >  2:246062/1‑100 (MQ=255)
                                                  aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                  >  2:1268023/1‑100 (MQ=255)
                                                   ccGGTTGAGCCCGTAAAATAAGCTTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTGGTCTGGTTCGCGGCCTTTCCAGCAGGtt                                                 <  2:2001193/100‑1 (MQ=25)
                                                    cGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTGTGGTTCGCGGCCTTTGCAGCAgg                                                   >  1:2158235/1‑97 (MQ=255)
                                                     ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                               >  2:114196/1‑100 (MQ=255)
                                                       ttGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACt                                             <  2:2085336/100‑1 (MQ=255)
                                                         gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtg                                           <  1:1160503/100‑1 (MQ=255)
                                                           gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                         >  2:2397177/1‑100 (MQ=255)
                                                            ccAGGGAAATAGGTATGTGCCAGGCCATAAAGAACTGCCGAATTTGGCGGTGTTTTATTGCGAGTTTGGTTCGCGGCCTTTCCAGCAGGTTGACGtgtg                                         <  1:2192038/99‑1 (MQ=21)
                                                             cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                       >  1:1237193/1‑100 (MQ=255)
                                                             cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                       <  2:1397131/100‑1 (MQ=255)
                                                                 gAAATAATTATTTTACAGGCAATAATCTATCGCCGGATTTGGCGGTTTTTTATTGCTAGTCTTGTTCGCGGCCTTTCCAGCAGGTTGACGTGTGTTa                                      <  1:1355319/97‑1 (MQ=25)
                                                                 gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGGATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                    <  1:1143607/99‑1 (MQ=255)
                                                                 gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                   >  1:2332811/1‑100 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                  <  2:579087/100‑1 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                  <  2:1177769/100‑1 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                  <  2:124729/100‑1 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                  <  1:2522234/100‑1 (MQ=255)
                                                                    aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAg                                >  2:897440/1‑100 (MQ=255)
                                                                    aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAg                                <  1:722394/100‑1 (MQ=255)
                                                                    aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAg                                <  2:507580/100‑1 (MQ=255)
                                                                     tAAGTATTTTACAGGCAATAAATATCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                               <  2:528083/100‑1 (MQ=255)
                                                                     tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                               >  1:466769/1‑100 (MQ=255)
                                                                      aaGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTATCGTTCGCGGCCTTTCGAGCATGTTTCCTTGTGTTACATGAGc                               >  1:1017441/1‑99 (MQ=37)
                                                                        gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAAc                            >  2:654030/1‑100 (MQ=255)
                                                                          aTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGc                          >  1:684315/1‑100 (MQ=255)
                                                                           ccagactagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa                         >  1:224191/9‑100 (MQ=25)
                                                                             ttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                       <  2:1163896/100‑1 (MQ=255)
                                                                             ttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                       <  1:1961559/100‑1 (MQ=255)
                                                                               aCAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTg                     <  2:1877286/100‑1 (MQ=255)
                                                                                ctagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGc                    >  2:1117394/4‑100 (MQ=37)
                                                                                cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGc                    <  1:1177375/100‑1 (MQ=255)
                                                                                cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTAACATGAGCAACGCGGGTGc                    >  2:2199993/1‑100 (MQ=255)
                                                                                 ccgatctAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                   <  1:54043/94‑1 (MQ=255)
                                                                                 aGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                   <  1:2490977/100‑1 (MQ=255)
                                                                                 aGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                   <  2:192805/100‑1 (MQ=255)
                                                                                  ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                  >  2:663912/1‑100 (MQ=255)
                                                                                   gCAATAAAAAACCGCCGCATTTGGCGGTTTTTTATGGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                 <  1:321186/100‑1 (MQ=255)
                                                                                   gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACGTGTGTTACATGAGCAACGCAGGTGCTTc                 >  1:169468/1‑100 (MQ=255)
                                                                                    cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                >  2:2035879/1‑100 (MQ=255)
                                                                                     aaTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCAGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcc                >  1:2040430/1‑98 (MQ=37)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             <  1:675942/100‑1 (MQ=255)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             <  1:1916948/100‑1 (MQ=255)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             >  2:2588442/1‑100 (MQ=255)
                                                                                        aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc            >  2:891843/1‑100 (MQ=255)
                                                                                        aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc            <  1:509230/100‑1 (MQ=255)
                                                                                        aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCCGGATCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtagatcgg            >  2:54044/1‑93 (MQ=37)
                                                                                        aaaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAAGATGACATGTGCGACATGAGCGGCGCAGGTGCGGCACAg             >  1:1838200/1‑99 (MQ=18)
                                                                                         aaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGTCTTATGTTAGTTGAGCAACGCAGGTGCTTCGCAg             >  1:2288290/1‑98 (MQ=37)
                                                                                         aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa           >  2:793611/1‑100 (MQ=255)
                                                                                         aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa           >  2:237077/1‑100 (MQ=255)
                                                                                            aaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaaa        <  1:2470386/100‑1 (MQ=255)
                                                                                              ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc       <  1:1906007/99‑1 (MQ=255)
                                                                                               cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaa     <  2:2445178/100‑1 (MQ=255)
                                                                                                gccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACCTGAGCAACGCAGGTGCTTCACAGCAAAACAAt    >  2:900047/1‑100 (MQ=255)
                                                                                                gccgAAGTTGGCGGATTTTGATGGCTAGTAAGGTCCGCGGCCTTACCAGCAGGTTGGCTTGTGTTACATGAGCAACGAAGGTACTTCACAGCAAAACa      <  1:1781598/98‑1 (MQ=255)
                                                                                                  cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTGCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGAAAAACAATAc  <  1:121564/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAC  >  NC_000913/2653268‑2653465

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: