Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,372 T→C 15.1% intergenic (+33/+785) 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,3720TC15.1% 268.7 / 37.1 119intergenic (+33/+785)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base T (39/61);  new base C (10/8);  total (50/69)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 126 total aligned reads displayed.

ACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACC  >  NC_000913/2653273‑2653470
                                                                                                   |                                                                                                  
aCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAttt                                                                                                    <  2:1855644/100‑1 (MQ=255)
 cTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg                                                                                                   >  2:2334083/1‑100 (MQ=255)
  tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGTAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                                  <  1:1704596/100‑1 (MQ=255)
   gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                                 <  1:2934210/100‑1 (MQ=255)
   gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                                 <  2:613603/100‑1 (MQ=255)
   gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                                 <  1:173331/100‑1 (MQ=255)
     acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGTTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                               <  1:1987060/100‑1 (MQ=255)
     acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                               <  1:562202/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)
          ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                          <  2:705601/100‑1 (MQ=255)
          ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                          <  2:44568/100‑1 (MQ=255)
          ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                          <  2:349366/100‑1 (MQ=255)
           gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                         <  2:1669211/100‑1 (MQ=255)
            cgcATGGAGTGAATGGGGCGCGCGGTCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                          <  1:1244004/98‑1 (MQ=255)
            cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACTGTTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                        <  1:143903/100‑1 (MQ=38)
            cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTa                                                                                        >  1:369716/1‑100 (MQ=255)
             gcATGGAGTGAGTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                       <  2:1561375/100‑1 (MQ=38)
             gcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCTATTTAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAt                                                                                       <  1:1003820/100‑1 (MQ=255)
               aTGGAGTGAATGGGGCGCGCGGGCAGATTTACTGTTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                     <  1:398519/100‑1 (MQ=255)
               aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                     <  1:2910100/100‑1 (MQ=255)
                 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)
                           gggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTc                                                                         >  2:2830740/1‑100 (MQ=255)
                              cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                      >  1:2416826/1‑100 (MQ=255)
                               gcgcgGGCAGATTTACCGGTTGAGTCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCCGGTTCGCgg                                                                     <  1:2450245/100‑1 (MQ=255)
                               gcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCgg                                                                     <  1:2173725/100‑1 (MQ=255)
                                cgcgGGCAGATTTACCGTTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                                                                    <  2:727448/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                                                              <  2:1190567/100‑1 (MQ=255)
                                        gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                            <  1:916749/100‑1 (MQ=255)
                                        gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGACAATAAATAATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                            >  1:798863/1‑100 (MQ=255)
                                             aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                       <  1:2621877/100‑1 (MQ=255)
                                             aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                       >  2:792973/1‑100 (MQ=255)
                                                ttttGAGTCAGTTAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                       <  1:555630/95‑1 (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)
                                                      gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                              <  2:2660013/100‑1 (MQ=255)
                                                      gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                              >  1:1162706/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)
                                                           tGAAATAAGTATTTTACAGGCAATAAAAAACCACCGAATATGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                         <  2:2127675/100‑1 (MQ=255)
                                                            gAAATAAGTATTTTACAGGCAATAAAAAACTGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                        <  1:360982/100‑1 (MQ=255)
                                                            gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                        <  2:3006684/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)
                                                                    tATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACg                                <  1:789398/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)
                                                                              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)
                                                                                  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)
                                                                                     aaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTACAGAAAGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaa               >  1:1091200/1‑100 (MQ=38)
                                                                                        aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc            >  1:1970420/1‑100 (MQ=255)
                                                                                         ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc            <  1:1845535/99‑1 (MQ=255)
                                                                                         ccgccgAATTTGGCGGTTTTTTATTGCTAGTATGGTTCGCGGCCATTCCTGCTGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCTCAGGaaatc            <  2:2112753/99‑3 (MQ=255)
                                                                                          cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaa          >  2:2214746/1‑100 (MQ=255)
                                                                                          cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCATTTGCAGCGGGTTGACTTGTGTTAGATGAGAATCGCTGGTGATTCACAGCCAAACa           >  2:154206/1‑99 (MQ=255)
                                                                                            ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa        >  2:2740988/1‑100 (MQ=255)
                                                                                             cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc       >  1:1848872/1‑100 (MQ=255)
                                                                                             cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAAATGGGCAACGCAGGTGCTTCACAGCAAAACAATAc       >  2:2171872/1‑100 (MQ=255)
                                                                                              aaaTTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACt      >  2:251554/2‑100 (MQ=255)
                                                                                                aTTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCa    <  2:1091425/100‑1 (MQ=255)
                                                                                                  ttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc  <  2:2988237/100‑1 (MQ=255)
                                                                                                  ttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc  <  2:355346/100‑1 (MQ=255)
                                                                                                  ttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc  <  2:2438559/100‑1 (MQ=255)
                                                                                                  ttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc  <  1:293687/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACC  >  NC_000913/2653273‑2653470

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: