Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 15.6% 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.6% 261.3 / 32.5 113intergenic (+28/+790)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (35/59);  new base A (9/9);  total (44/69)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

Only 100 of 122 total aligned reads displayed.

GTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAC  >  NC_000913/2653268‑2653465
                                                                                                   |                                                                                                  
     aCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAttt                                                                                               <  2:1855644/100‑1 (MQ=255)
     aCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAAttc                                                                                               <  2:1284140/100‑2 (MQ=255)
      cTGTACGACGGCGCATGGATTAAATGTAGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAATTATTTTACAGGCAATAAAAAACCGCCGAAttt                                                                                               <  2:649932/99‑1 (MQ=255)
        gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                            <  1:2934210/100‑1 (MQ=255)
        gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                            <  1:173331/100‑1 (MQ=255)
          acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGTTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          <  1:1987060/100‑1 (MQ=255)
           cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGt                                                                                         >  2:241258/1‑100 (MQ=255)
           cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGt                                                                                         <  1:418814/100‑1 (MQ=255)
            gacgGCTCATGTAGTGAATGGGGCGCGCGTTCAGATTTACTGGTTTAACCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          <  1:1635847/98‑1 (MQ=255)
             acgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttt                                                                                       <  1:2320794/100‑1 (MQ=255)
              cgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttt                                                                                      <  1:1454675/100‑1 (MQ=255)
               ggCGCATGGATTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                     <  1:2016557/100‑1 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                     <  2:705601/100‑1 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                     <  2:349366/100‑1 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                     <  2:44568/100‑1 (MQ=255)
                gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                    <  2:1669211/100‑1 (MQ=255)
                 cgcATGGAGTGAATGGGGCGCGCGGTCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                     <  1:1244004/98‑1 (MQ=255)
                 cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTa                                                                                   >  1:369716/1‑100 (MQ=255)
                  gcATGGAGTGAGTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                  <  2:1561375/100‑1 (MQ=38)
                  gcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCTATTTAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAt                                                                                  <  1:1003820/100‑1 (MQ=255)
                   cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                 >  1:671591/1‑100 (MQ=255)
                    aTGGAGTGAATGGGGCGCGCGGGCAGATTTACTGTTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                <  1:398519/100‑1 (MQ=255)
                    aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                <  1:2910100/100‑1 (MQ=255)
                    aTGGAGTGAATGCGGCGCTCGGGCAGATTTACCGGTTGAGGCAGTGAAATAAGTATTTTACAGGCAATAAAAATCCGCTAAATTCGGCGGTTTTTTAt                                                                                  <  2:2568432/98‑1 (MQ=37)
                      ggAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGcc                                                                              <  1:1407541/100‑2 (MQ=255)
                       gAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTa                                                                             >  2:1591768/1‑100 (MQ=255)
                         gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGt                                                                           >  1:2745788/1‑100 (MQ=255)
                          tGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTc                                                                          <  1:2577596/100‑1 (MQ=255)
                           gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                         <  2:1519927/100‑1 (MQ=255)
                           gAATGGGACGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                         >  1:1583120/1‑100 (MQ=255)
                                gggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTc                                                                    >  2:2830740/1‑100 (MQ=255)
                                   cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                 >  1:2416826/1‑100 (MQ=255)
                                    gcgcgGGCAGATTTACCGGTTGAGTCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCCGGTTCGCgg                                                                <  1:2450245/100‑1 (MQ=255)
                                     cgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                                                               <  1:2212023/100‑1 (MQ=255)
                                     cgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                                                               <  2:1431788/100‑1 (MQ=255)
                                      gcgGGCAGATTTTCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                                                               <  1:1429194/99‑1 (MQ=255)
                                       cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGGATTTTACAGGCAATAAAAAACCGCCGAATTGGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                             >  1:2127698/1‑100 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         <  1:852712/100‑1 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         <  2:2961896/100‑1 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         <  2:1190567/100‑1 (MQ=255)
                                             gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       <  1:916749/100‑1 (MQ=255)
                                             gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGACAATAAATAATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       >  1:798863/1‑100 (MQ=255)
                                              aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagc                                                      <  2:369713/100‑1 (MQ=255)
                                                  aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                  <  1:2621877/100‑1 (MQ=255)
                                                  aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                  <  2:2464292/100‑1 (MQ=255)
                                                  aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                  >  2:792973/1‑100 (MQ=255)
                                                   ccGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGtt                                                 >  1:1159055/1‑100 (MQ=255)
                                                     ttttGAGTCAGTTAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                  <  1:555630/95‑1 (MQ=255)
                                                     ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                               >  1:355349/1‑100 (MQ=255)
                                                      gTTGAGTCAGTGAAATAAGTTTTTTATAGTTAATAAAAATCCGCTGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAc                                              <  1:2638737/100‑1 (MQ=25)
                                                        tGAGTCAGTGAAATAAGTATTTTATAGGCAATAAAAGTTCGTCTAATTTGGCGGTTGTTTATTGCTAGTCGGGTTCGCGGCCTTTCCAGCAGGTTGACtt                                            <  1:676873/100‑1 (MQ=25)
                                                        tGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACtt                                            >  2:570811/1‑100 (MQ=255)
                                                          aGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                          >  2:293684/1‑100 (MQ=255)
                                                           gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                          >  1:624779/1‑99 (MQ=255)
                                                             cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTCGTATGGTTCGCCGCCTTTACAGCAGGTTGACTTGTGtt                                       >  2:2809977/1‑100 (MQ=38)
                                                              aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGATAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTa                                      <  1:1127924/100‑1 (MQ=255)
                                                                tgatattattATTTTACAGGTAATAAAAAACTGTCGAATTTGGTGGTTTTTTATTGCTAGTCTGGTTCGCGGACTTTCCAGCAGGTTGACTTGTGTTACa                                    <  1:97456/91‑1 (MQ=37)
                                                                tGAAATAAGTATTTTACAGGCAATAAAAAACCGTCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                    <  1:2772482/100‑1 (MQ=255)
                                                                tGAAATAAGTATTTTACAGGCAATAAAAAACCACCGAATATGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                    <  2:2127675/100‑1 (MQ=255)
                                                                 gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                   <  2:1465554/100‑1 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                  <  2:253809/100‑1 (MQ=255)
                                                                    aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAATCTGATTAGCGGACTTTCCAGCAGGTTGACTTGTGTTACATGa                                 >  1:1768362/1‑99 (MQ=38)
                                                                     tAAGTATTTTACAGGCAATAAAAATCTGCCGTTTTTTGTTGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAg                                <  1:65239/99‑1 (MQ=25)
                                                                      aaGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCa                              <  1:2313629/100‑1 (MQ=255)
                                                                       aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCaa                             >  1:2168144/1‑100 (MQ=255)
                                                                        gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGATAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAAc                            >  2:2825370/1‑100 (MQ=255)
                                                                           ttttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa                         <  1:1829276/100‑1 (MQ=255)
                                                                            tttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAg                        <  2:919864/100‑1 (MQ=255)
                                                                             ttACAGGTAATAAAAAACCGTCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGGGGCCTTTCTAGCAGGTTGAGTTGTGTTACATGAGCAACGCAgg                       <  2:232112/100‑1 (MQ=38)
                                                                             ttACAGGCAATAAAAAACAGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGAAGGTTGACTTGTGTTACATGAGCAACGAAgg                       >  2:2502667/1‑100 (MQ=255)
                                                                                  ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                  <  2:624775/100‑1 (MQ=255)
                                                                                  ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                  >  2:775321/1‑100 (MQ=255)
                                                                                   gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                 <  2:2067259/100‑1 (MQ=255)
                                                                                   gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                 >  2:1275079/1‑100 (MQ=255)
                                                                                    cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                >  2:276862/1‑100 (MQ=255)
                                                                                    cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                 <  1:2334110/99‑1 (MQ=255)
                                                                                    cAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                >  1:2449642/1‑100 (MQ=255)
                                                                                     aaTAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcac               <  1:1572895/100‑1 (MQ=255)
                                                                                     aaTAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAACAGGTTGAATAGTGATACATGAGCAACGCAGGTGCTTca                >  1:1138974/1‑99 (MQ=37)
                                                                                     aaTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcac               >  2:716980/1‑100 (MQ=255)
                                                                                      aTAACAAACCGCCGAATTTAGCGAATTTTTATTGCTAGTGTGGTTCGCGGCCTTTCCAGCACGTTGACTTATTTTACATAAACAATGCATGTGCtt                  >  2:884875/1‑96 (MQ=19)
                                                                                      aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAGTTGTGTTACATCATCAACGGAAGGGCTTtaca              >  2:699555/1‑100 (MQ=18)
                                                                                      aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca              >  1:1415359/1‑100 (MQ=255)
                                                                                       tgaaaatCCGCCGCATTTGGCGTTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcac               <  2:2914278/96‑1 (MQ=38)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGATTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             >  1:333929/1‑100 (MQ=255)
                                                                                       tAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAATTGTATTACAGAAGCAAGGCAGGTGGttaa                >  2:2954977/1‑95 (MQ=18)
                                                                                       tAAAAAACCGCCAAATTCGAAGGAATTTTATTGCCTGTCTGGTTCACGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGAAACGCAGGTGCTTCACAg             >  2:1291233/1‑100 (MQ=25)
                                                                                        aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc            >  1:240342/1‑100 (MQ=255)
                                                                                         aaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa           >  1:2183598/1‑100 (MQ=255)
                                                                                         aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAAGAGCAACGCAGGCGCTTCACAg             >  2:1947994/1‑98 (MQ=25)
                                                                                         aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAAGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa           >  1:1173130/1‑100 (MQ=37)
                                                                                          aaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACACAGGTGCTTCACAGCaa          >  1:551030/1‑100 (MQ=255)
                                                                                          aaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTACAGAAAGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaa          >  1:1091200/1‑100 (MQ=38)
                                                                                              ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc       <  1:1845535/99‑1 (MQ=255)
                                                                                               cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaa     >  2:2214746/1‑100 (MQ=255)
                                                                                               cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCATTTGCAGCGGGTTGACTTGTGTTAGATGAGAATCGCTGGTGATTCACAGCCAAACa      >  2:154206/1‑99 (MQ=255)
                                                                                               cgccgAATTTGGCGGTTTGTTATTGATAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTATGTTACATAAGCAACGCAGGTGCTTCACAGCAAAACaa     >  2:2081930/1‑100 (MQ=38)
                                                                                                 ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa   >  2:2740988/1‑100 (MQ=255)
                                                                                                  cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc  >  1:1848872/1‑100 (MQ=255)
                                                                                                  cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAAATGGGCAACGCAGGTGCTTCACAGCAAAACAATAc  >  2:2171872/1‑100 (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 < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: