Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 16.7% 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,3670GA16.7% 408.2 / 56.8 166intergenic (+28/+790)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (56/82);  new base A (16/12);  total (72/94)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

TGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACT  >  NC_000913/2653269‑2653466
                                                                                                  |                                                                                                   
tGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGa                                                                                                    >  1:3484557/1‑100 (MQ=255)
tGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGa                                                                                                    >  1:2305187/1‑100 (MQ=255)
 gAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGaa                                                                                                   >  2:473151/1‑100 (MQ=255)
    aCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAttt                                                                                                >  1:2615190/1‑100 (MQ=255)
     cTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg                                                                                               <  2:1642314/100‑1 (MQ=255)
       gtgggagggCGCGGGGAGTGAATGGGGAGGGGGGGCAGATTTACCGGTTGAGCGAGTGAAATAAGTATTTGACAGGCAATAAAAAACCGCCGAATTTg                                                                                               <  2:4015405/91‑1 (MQ=255)
       gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                             <  2:620220/100‑1 (MQ=255)
         acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                           >  1:3374026/1‑100 (MQ=255)
           gacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtt                                                                                         >  2:216056/1‑100 (MQ=255)
           gacgGAGCATGTAGTGAATGGGGGGCGCGGTCAGATTTACCGGTTGAGCCATTGATATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtt                                                                                         <  2:3307884/100‑1 (MQ=37)
            acgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttt                                                                                        >  1:3270692/1‑100 (MQ=255)
             cgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttt                                                                                       <  1:721920/100‑1 (MQ=255)
              ggCGCATGGAGTGAATGGGGCGCGCGTGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                      <  2:656375/100‑1 (MQ=255)
              ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                      >  1:625737/1‑100 (MQ=255)
              ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                      <  1:4001120/100‑1 (MQ=255)
              ggCGCATGGAGTGAATCGGGCGCGCCGGCAGATTTACCTGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                      <  2:89166/100‑1 (MQ=38)
               gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                     <  1:500382/100‑1 (MQ=255)
                cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                    <  1:745564/100‑1 (MQ=255)
                 gcATGCAGTGAATGCGGCGCGCGAGCACATTTACCAGTTGAGCCAGTCAAATACGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                   <  2:3128243/100‑1 (MQ=25)
                  naTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                  >  2:2936374/2‑100 (MQ=255)
                  cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                  >  1:1327246/1‑100 (MQ=255)
                  cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAtt                                                                                  <  2:157140/100‑1 (MQ=255)
                   aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACGGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTg                                                                                 <  2:2013274/100‑1 (MQ=255)
                   aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                 <  1:3241203/100‑1 (MQ=255)
                   aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTg                                                                                 <  1:2575347/100‑1 (MQ=255)
                    tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGc                                                                                <  2:1063155/100‑1 (MQ=255)
                     ggAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGcc                                                                               <  1:3203384/100‑2 (MQ=255)
                      gAGTGAATGGGGGGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGGATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                 <  2:494035/97‑1 (MQ=255)
                        gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGCATTTGGCGGTTTTTTATTGCTAga                                                                            >  1:2481934/1‑99 (MQ=255)
                        gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGt                                                                            >  2:4084455/1‑100 (MQ=255)
                        gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGt                                                                            >  1:1265320/1‑100 (MQ=255)
                        gTGAATGGGGCGCGCGGGCAAATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGt                                                                            >  2:2549513/1‑100 (MQ=255)
                          gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                          >  2:630929/1‑100 (MQ=255)
                          gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                          >  2:3107959/1‑100 (MQ=255)
                          gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                          >  1:3501657/1‑100 (MQ=255)
                          gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGcctgtaa                                                                          <  1:2034140/100‑7 (MQ=25)
                           aaTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTg                                                                         <  1:3798895/100‑1 (MQ=255)
                            aTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg                                                                        <  2:3185281/100‑1 (MQ=255)
                            aTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg                                                                        <  1:2349307/100‑1 (MQ=255)
                                  cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                  <  2:249695/100‑1 (MQ=255)
                                     gcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                               <  1:101815/100‑1 (MQ=255)
                                          cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                          >  2:3006934/1‑100 (MQ=255)
                                           aGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                         >  2:382498/1‑100 (MQ=255)
                                            gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                        >  2:1401494/1‑100 (MQ=255)
                                             aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagc                                                       >  1:3567697/1‑100 (MQ=255)
                                               ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                     <  1:2509351/100‑1 (MQ=255)
                                                tACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAgg                                                    <  1:1742815/100‑1 (MQ=255)
                                                 aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                   <  1:3490212/100‑1 (MQ=255)
                                                   cGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg                                                 <  2:1668789/100‑1 (MQ=255)
                                                   cGGTTGAGCCAGTGAAAGAAGTATTGTACAGGCAATAAAAAACGGCCGTATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCCGGtt                                                  <  2:3424795/99‑1 (MQ=37)
                                                    ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                                >  1:597384/1‑100 (MQ=255)
                                                    ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                                >  1:1566041/1‑100 (MQ=255)
                                                          cttccgatcttAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                          <  2:2619304/90‑1 (MQ=255)
                                                            cAGTGCAATAAGTATTTGACAGGCAATAAATTACCGCCGAATTTGGCGGTTTTTTATTTCTAGTCTGGTTCCCCGCCTTTCAGCAGGTTGACTTGTGTTa                                       <  2:3934711/100‑1 (MQ=255)
                                                             aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTa                                       <  1:2438715/100‑1 (MQ=255)
                                                             aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGCGTTa                                       >  2:2327322/1‑100 (MQ=255)
                                                              gTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc                                      <  2:3563615/100‑1 (MQ=255)
                                                              gTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc                                      <  1:1868595/100‑1 (MQ=255)
                                                              gTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc                                      <  1:827553/100‑1 (MQ=255)
                                                                gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                    >  2:2253856/1‑100 (MQ=255)
                                                                 aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                   <  2:1150025/100‑1 (MQ=255)
                                                                  aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGa                                  <  1:3044356/100‑1 (MQ=255)
                                                                  aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGa                                  >  2:1327990/1‑100 (MQ=255)
                                                                  aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGa                                  <  1:3226188/100‑1 (MQ=255)
                                                                    tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtgagatcggaag                                >  1:2619319/1‑90 (MQ=255)
                                                                      aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCaa                              >  1:3696973/1‑100 (MQ=255)
                                                                      aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCaa                              >  1:3696599/1‑100 (MQ=255)
                                                                        tATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTCACATGAGCAACg                            >  2:74657/1‑100 (MQ=255)
                                                                         aTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGc                           >  1:1679242/1‑100 (MQ=255)
                                                                          ttttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa                          <  1:2317343/100‑1 (MQ=255)
                                                                            ttACAGGCAATAAAAGACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGc                           <  1:4121149/97‑1 (MQ=255)
                                                                            ttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgn                        <  1:2936391/100‑2 (MQ=255)
                                                                            ttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                        <  1:1366957/100‑1 (MQ=255)
                                                                              aCAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTg                      <  2:2792444/100‑1 (MQ=255)
                                                                              aCAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTg                      <  2:3664663/100‑1 (MQ=255)
                                                                              aCAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTg                      <  2:1539193/100‑1 (MQ=255)
                                                                               cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCTTTTCCAGCAGCTTGACTTGTGTTACATGCGCAACGCAGGtac                     <  1:2818103/100‑3 (MQ=38)
                                                                                  gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                  <  2:2644289/100‑1 (MQ=255)
                                                                                   cAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                 >  1:1947647/1‑100 (MQ=255)
                                                                                   cAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                 >  1:2837497/1‑100 (MQ=255)
                                                                                   cAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                 >  2:4052404/1‑100 (MQ=255)
                                                                                   cAATAAAAAACCGCCAAATTCGGCGGTTTTATAGTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGGAGGTGCTTca                 >  2:2734076/1‑100 (MQ=25)
                                                                                     aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca               >  1:2247955/1‑100 (MQ=255)
                                                                                     aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca               >  1:3394672/1‑100 (MQ=255)
                                                                                      tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg              <  1:2598308/100‑1 (MQ=255)
                                                                                      tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTGTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg              <  2:1983348/100‑1 (MQ=255)
                                                                                      tAAAAAACCGCCAAATTCGGCGGGTTTTTATTGATAGTATGGTTCGCGGCCTTTCAAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTg                      >  2:3296216/1‑92 (MQ=37)
                                                                                       aaaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc             >  1:3738437/1‑100 (MQ=255)
                                                                                       aaaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc             >  1:991880/1‑100 (MQ=255)
                                                                                       aaaaaaCCGCCAAATTCGGCGGTTTTTTATTCCTAGTCTGGTGCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc             >  1:2811743/1‑100 (MQ=37)
                                                                                        aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa            >  1:2433877/1‑100 (MQ=255)
                                                                                         aaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaa           >  2:3082109/1‑100 (MQ=255)
                                                                                          aaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaa          >  1:2222484/1‑100 (MQ=255)
                                                                                            aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc        <  1:1495275/100‑1 (MQ=255)
                                                                                             ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGCCTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc        <  2:3717540/99‑1 (MQ=255)
                                                                                             ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACCTGCGCAACGCAGGTGCTTCACAGCAAAACa       >  2:2704405/1‑100 (MQ=255)
                                                                                                ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa    <  1:735060/100‑1 (MQ=255)
                                                                                                ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa    <  1:3826679/100‑1 (MQ=255)
                                                                                                ccgAATTTGGCGGTGTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGAAACGCAGGGGCTTCACAGCAAAGCAAt     >  2:3484439/1‑99 (MQ=38)
                                                                                                 cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc   >  1:3025083/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACT  >  NC_000913/2653269‑2653466

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: