Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,372 T→C 16.0% 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,3720TC16.0% 399.9 / 64.4 174intergenic (+33/+785)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base T (58/88);  new base C (16/12);  total (74/100)
Fisher's exact test for biased strand distribution p-value = 9.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 177 total aligned reads displayed.

ACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACC  >  NC_000913/2653273‑2653470
                                                                                                   |                                                                                                  
aCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAttt                                                                                                    >  1:2615190/1‑100 (MQ=255)
 cTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg                                                                                                   <  2:1642314/100‑1 (MQ=255)
   gtgggagggCGCGGGGAGTGAATGGGGAGGGGGGGCAGATTTACCGGTTGAGCGAGTGAAATAAGTATTTGACAGGCAATAAAAAACCGCCGAATTTg                                                                                                   <  2:4015405/91‑1 (MQ=255)
   gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                                 <  1:719007/100‑1 (MQ=255)
     acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                               <  1:601120/100‑1 (MQ=255)
     acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                               <  1:2410474/100‑1 (MQ=255)
       gacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtt                                                                                             >  1:2440867/1‑100 (MQ=255)
       gacgGAGCATGTAGTGAATGGGGGGCGCGGTCAGATTTACCGGTTGAGCCATTGATATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtt                                                                                             <  2:3307884/100‑1 (MQ=37)
         cgGCTCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                                <  2:3072390/95‑1 (MQ=255)
         cgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttt                                                                                           <  2:2239600/100‑1 (MQ=255)
          ggCGCATGGAGTGAATCGGGCGCGCCGGCAGATTTACCTGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                          <  2:89166/100‑1 (MQ=38)
           gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                         <  1:500382/100‑1 (MQ=255)
            cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTa                                                                                        <  2:684701/100‑1 (MQ=255)
            cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                        <  1:745564/100‑1 (MQ=255)
             gcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTCGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                       <  1:2435887/100‑1 (MQ=38)
             gcATGCAGTGAATGCGGCGCGCGAGCACATTTACCAGTTGAGCCAGTCAAATACGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                       <  2:3128243/100‑1 (MQ=25)
              cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                      >  1:1327246/1‑100 (MQ=255)
              cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAtt                                                                                      <  2:157140/100‑1 (MQ=255)
               aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                     >  1:2336909/1‑100 (MQ=255)
                tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                    >  1:326328/1‑100 (MQ=255)
                tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGc                                                                                    <  2:1063155/100‑1 (MQ=255)
                 ggAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                                   <  1:3032142/100‑1 (MQ=255)
                 ggAGTGAATGGGCCGCGCGTGCAGATCTACCGGTTAAGCCAGTGAAATGAGTATTTTTCATGCAATAAAAAACCGCCAAATTCTGCGGTTTTTTATTg                                                                                     <  2:2098126/98‑1 (MQ=21)
                    gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGCATTTGGCGGTTTTTTATTGCTAga                                                                                >  1:2481934/1‑99 (MQ=255)
                    gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGt                                                                                >  1:1265320/1‑100 (MQ=255)
                    gTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGt                                                                                >  2:4084455/1‑100 (MQ=255)
                      gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                              >  2:3107959/1‑100 (MQ=255)
                      gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                              <  2:3156520/100‑1 (MQ=255)
                      gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                              >  2:630929/1‑100 (MQ=255)
                      gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                              >  1:3941248/1‑100 (MQ=255)
                      gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGcctgtaa                                                                              <  1:2034140/100‑7 (MQ=25)
                        aTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg                                                                            >  1:3388607/1‑100 (MQ=255)
                           gggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTc                                                                         >  1:2588887/1‑100 (MQ=255)
                              cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                      <  2:249695/100‑1 (MQ=255)
                                 gcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                                   <  1:101815/100‑1 (MQ=255)
                                 gcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                                   <  1:2672265/100‑1 (MQ=255)
                                     gCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc                                                               <  1:2853588/100‑1 (MQ=255)
                                      cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                              <  1:895889/100‑1 (MQ=255)
                                       aGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                             >  2:382498/1‑100 (MQ=255)
                                        gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                            >  2:3447630/1‑100 (MQ=255)
                                         aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagc                                                           <  1:1916561/100‑1 (MQ=255)
                                          tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGGTCGCGGCCTTTCcagca                                                          >  1:3280032/1‑100 (MQ=255)
                                           ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                         <  2:1337855/100‑1 (MQ=255)
                                           ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                         <  2:2691298/100‑1 (MQ=255)
                                           ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                         <  1:3702322/100‑1 (MQ=255)
                                           ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                         >  2:1518922/1‑100 (MQ=255)
                                            tACCGGTTGAGCCATTGAAATAAGTATTTTCCAGGTAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                                   <  2:2988202/89‑1 (MQ=255)
                                            tACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAgg                                                        <  1:1742815/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)
                                                ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGATTTTTATTGCTAGTCTGGTTCGCGGCCTGTCCAGCAGGttta                                                    <  1:2813796/100‑3 (MQ=255)
                                                    gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTGTGGTTCGCGGCCTTTCCGGCAGGTTGACtcg                                                <  1:2818717/100‑3 (MQ=255)
                                                      gCCAGTGAAATAAGTATTTTCCAGGCAATAAAAAACCGCAGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                              >  2:2554275/1‑100 (MQ=255)
                                                      gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                              >  2:490011/1‑100 (MQ=255)
                                                      cttccgatcttAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                              <  2:2619304/90‑1 (MQ=255)
                                                      cgatctGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                              <  2:3491331/95‑1 (MQ=255)
                                                        cAGTGCAATAAGTATTTGACAGGCAATAAATTACCGCCGAATTTGGCGGTTTTTTATTTCTAGTCTGGTTCCCCGCCTTTCAGCAGGTTGACTTGTGTTa                                           <  2:3934711/100‑1 (MQ=255)
                                                        cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                            <  1:3601554/100‑1 (MQ=255)
                                                         aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGCGTTa                                           >  2:2327322/1‑100 (MQ=255)
                                                          gTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc                                          <  1:827553/100‑1 (MQ=255)
                                                          gTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc                                          <  1:1868595/100‑1 (MQ=255)
                                                            gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtgagatcg                                        >  1:3491357/1‑94 (MQ=255)
                                                            gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                        >  2:2253856/1‑100 (MQ=255)
                                                             aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                       <  2:1150025/100‑1 (MQ=255)
                                                              aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGa                                      >  2:1327990/1‑100 (MQ=255)
                                                              aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGa                                      <  1:3226188/100‑1 (MQ=255)
                                                                tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTTGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                                    <  2:356374/100‑1 (MQ=255)
                                                                tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                                    <  1:3773913/100‑1 (MQ=255)
                                                                 aaGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCa                                   <  1:3529286/100‑1 (MQ=255)
                                                                  aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCaa                                  >  1:3696599/1‑100 (MQ=255)
                                                                    tATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACg                                <  1:2549528/100‑1 (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)
                                                                             ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                       <  1:3482008/100‑1 (MQ=255)
                                                                              gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTTGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                      <  2:2909211/100‑1 (MQ=255)
                                                                               cAATAAAAAACCGCCAAATTCGGCGGTTTTATAGTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGGAGGTGCTTca                     >  2:2734076/1‑100 (MQ=25)
                                                                                 aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca                   >  1:3394672/1‑100 (MQ=255)
                                                                                 aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca                   >  1:2247955/1‑100 (MQ=255)
                                                                                  tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTGTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg                  <  2:1983348/100‑1 (MQ=255)
                                                                                    aaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa                >  2:550415/1‑100 (MQ=255)
                                                                                    aaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa                <  2:3186992/100‑1 (MQ=255)
                                                                                    aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGTCTTGTGGTACATGTGCAAGGCAGGTGCTTCACAg                  >  2:4072743/1‑98 (MQ=21)
                                                                                    aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCCGGTGCTTCGCAg                  >  2:643641/1‑98 (MQ=25)
                                                                                      aaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaa              >  1:2222484/1‑100 (MQ=255)
                                                                                        aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc            <  1:1495275/100‑1 (MQ=255)
                                                                                         ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGCCTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc            <  2:3717540/99‑1 (MQ=255)
                                                                                          cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaa          >  2:1624066/1‑100 (MQ=255)
                                                                                           gccgAATTTGGCGGTTTTTTATCGCTAGGCTAGTTCGCGGCCTTTTAAGCAGGTTGACTTGTTTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAAt         <  1:2601659/100‑1 (MQ=255)
                                                                                            ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa        <  1:3826679/100‑1 (MQ=255)
                                                                                            ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa        <  1:735060/100‑1 (MQ=255)
                                                                                            ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa        >  1:3216442/1‑100 (MQ=255)
                                                                                              gAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGGAGGTGCTTCACAGCAAAACAATACt      >  1:310179/1‑100 (MQ=255)
                                                                                                  ttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc  <  2:3789166/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: