Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,105 G→C 15.1% W204C (TGG→TGC) ‡ sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,1050GC15.1% 537.4 / 78.8 192W204C (TGG→TGCsseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (83/80);  new base C (13/16);  total (96/96)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 198 total aligned reads displayed.

TATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGC  >  NC_000913/2653006‑2653204
                                                                                                   |                                                                                                   
tattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg                                                                                                     <  1:2337472/100‑1 (MQ=255)
tattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg                                                                                                     >  2:2276207/1‑100 (MQ=255)
 attGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGa                                                                                                    <  2:4184821/100‑1 (MQ=255)
    gATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACgg                                                                                                 >  1:3700576/1‑100 (MQ=255)
    gATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACgg                                                                                                 <  1:2726999/100‑1 (MQ=255)
     aTGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGa                                                                                                <  2:3899966/100‑1 (MQ=255)
        cccgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGaa                                                                                               >  1:2576196/1‑98 (MQ=255)
        cccgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACt                                                                                             >  2:1660446/1‑100 (MQ=255)
          cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTgg                                                                                           <  1:3626724/100‑1 (MQ=255)
          cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTgg                                                                                           <  1:2064341/100‑1 (MQ=255)
          cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcatt                                                                                           <  1:1611686/100‑4 (MQ=255)
           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattc                                                                                          <  1:475646/100‑5 (MQ=255)
           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattc                                                                                          <  2:1040923/100‑5 (MQ=255)
           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattc                                                                                          <  2:3517175/100‑5 (MQ=255)
             ccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcag                                                                                        <  2:3751599/100‑7 (MQ=255)
              cgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTTGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcagt                                                                                       <  2:3425695/100‑8 (MQ=255)
              cgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcg                                                                                       <  1:3213183/100‑1 (MQ=255)
              cgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcagt                                                                                       <  1:4100246/100‑8 (MQ=255)
              cgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATGCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcg                                                                                       <  2:2586730/100‑1 (MQ=255)
               ggCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcgc                                                                                      <  1:2360953/100‑1 (MQ=255)
               ggCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcgc                                                                                      <  2:2926563/100‑1 (MQ=255)
               ggCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcagtg                                                                                      <  1:2299829/100‑9 (MQ=255)
                gCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGcgcg                                                                                     <  1:246982/100‑1 (MQ=255)
                 cTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGa                                                                                    >  1:1076202/1‑100 (MQ=255)
                    cACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAgg                                                                                 >  2:1322161/1‑100 (MQ=255)
                     aCGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGc                                                                                <  2:2128836/100‑1 (MQ=255)
                      cGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATGCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCg                                                                               <  2:174015/100‑1 (MQ=255)
                        ttttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACTTCGCGGCAATATTCCCGGAGCCCTGAATGTTCCGTGGACGGGACTGGTGCGCGAAGGc                                                                                >  2:1780122/1‑97 (MQ=255)
                        ttttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGaa                                                                             <  2:3164694/100‑1 (MQ=255)
                         tttAACGCAGAAGTTGATGAACCTCGCTCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGaan                                                                            <  2:653072/100‑2 (MQ=255)
                         tttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAAc                                                                            <  1:1020033/100‑1 (MQ=255)
                         tttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAAc                                                                            <  1:376182/100‑1 (MQ=255)
                          ttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt                                                                           >  2:2373583/1‑100 (MQ=255)
                           tAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTa                                                                          >  1:707832/1‑100 (MQ=255)
                               gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGCGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                      >  2:2996847/1‑100 (MQ=255)
                               gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGCAGGCGAGCTaaaa                                                                       >  2:2601413/1‑99 (MQ=255)
                               gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                      <  1:1995148/100‑1 (MQ=255)
                               gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAAATGGTGCGCGAAGGCGAACTaaaaa                                                                      >  2:2081573/1‑100 (MQ=255)
                                  gAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGa                                                                   <  1:3185255/100‑1 (MQ=255)
                                    aGTTGATGAACCTCGCCCAGGTTTACGTCGCGTACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAcc                                                                 <  2:2943905/100‑1 (MQ=255)
                                        gatgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATg                                                             <  2:1081879/100‑1 (MQ=255)
                                        gatgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATg                                                             >  1:348442/1‑100 (MQ=255)
                                         attaaCCTCGCCCAGGTTTACGGCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATg                                                             <  2:2492433/96‑1 (MQ=255)
                                         atgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGa                                                            >  2:3747946/1‑100 (MQ=255)
                                         atgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGa                                                            >  2:450376/1‑100 (MQ=255)
                                                tCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                     >  1:4158163/1‑100 (MQ=255)
                                                tCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                     <  2:2408812/100‑1 (MQ=255)
                                                  gCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGc                                                   <  1:648847/100‑1 (MQ=255)
                                                   cccAGGTTTACGTCGCGGACATATTCCCGGAGCCCTGCATGTTTCGTGGACGGAACTGGTGCGGGAAGGCGAAATAAAAAGGACCGATGAACTGGATGc                                                   >  2:821529/1‑99 (MQ=255)
                                                   cccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCg                                                  >  2:364334/1‑100 (MQ=255)
                                                   cccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAAATGGATGCg                                                  >  2:729713/1‑100 (MQ=255)
                                                     cAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGat                                                >  2:2669347/1‑100 (MQ=255)
                                                     cAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGat                                                <  1:1146570/100‑1 (MQ=255)
                                                       ggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGGACTAAAAACGACCGATGAACTGGATGCGatat                                              >  2:2904494/1‑100 (MQ=255)
                                                       ggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGatat                                              >  2:3802880/1‑100 (MQ=255)
                                                         tttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttt                                            <  1:1679907/100‑1 (MQ=255)
                                                         tttACGTCGCGGACATATGCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtt                                             <  2:3591487/99‑1 (MQ=255)
                                                          ttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttt                                           <  2:3132361/100‑1 (MQ=255)
                                                           tACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGCTAttttt                                          >  2:632529/1‑100 (MQ=255)
                                                           tACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                          <  2:926052/100‑1 (MQ=255)
                                                           tACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                          >  1:352393/1‑100 (MQ=255)
                                                            aCGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                         <  1:1432262/100‑1 (MQ=255)
                                                            aCGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                         >  1:871876/1‑100 (MQ=255)
                                                             cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                        <  1:3086643/100‑1 (MQ=255)
                                                             cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                        <  1:3277650/100‑1 (MQ=255)
                                                             cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                        <  1:358277/100‑1 (MQ=255)
                                                              gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGatagatcggaa                                       >  1:2197914/1‑92 (MQ=255)
                                                              gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTgg                                       <  2:121328/100‑1 (MQ=255)
                                                                cgcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGCAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTc                                     >  2:3032197/1‑100 (MQ=255)
                                                                  cgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTcgc                                   >  2:3808839/1‑100 (MQ=255)
                                                                     aCATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggc                                <  1:1362559/100‑1 (MQ=255)
                                                                        tatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACTGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTc                             >  1:1893567/1‑100 (MQ=255)
                                                                         nttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                            >  2:3268110/2‑100 (MQ=255)
                                                                         atTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                            >  2:1619723/1‑100 (MQ=255)
                                                                         atTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                            >  1:3292173/1‑100 (MQ=255)
                                                                          ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAg                           <  1:2927629/100‑1 (MQ=255)
                                                                             ccGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTa                        <  1:4061608/100‑1 (MQ=255)
                                                                                 aGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAc                    <  1:1326648/100‑1 (MQ=255)
                                                                                 aGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAc                    <  1:264529/100‑1 (MQ=255)
                                                                                 aGCACTGAATGTTCCGGGGACGGAACTGGTGGGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAc                    <  2:1934951/100‑1 (MQ=255)
                                                                                   gcgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaa                  >  1:393063/9‑100 (MQ=255)
                                                                                    cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                 >  1:2053191/8‑100 (MQ=255)
                                                                                    cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                 >  1:3006592/8‑100 (MQ=255)
                                                                                    cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                 >  1:1274177/8‑100 (MQ=255)
                                                                                     gcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAc                >  1:3807007/7‑100 (MQ=255)
                                                                                      tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAcc               <  2:1459698/100‑1 (MQ=255)
                                                                                      caccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAcc               >  1:3693290/6‑100 (MQ=255)
                                                                                      caccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAcc               >  1:1016560/6‑100 (MQ=255)
                                                                                       gAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCa              <  2:455932/100‑1 (MQ=255)
                                                                                       gAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTCCGACAAACCaa             >  1:386814/1‑100 (MQ=255)
                                                                                       gAATGTTCCGTGGACGGAACTGGTGCGCGAAGGAGAACTAAAAACGACCGATGAACTGGATGCGATATTTGTTGGTCGCGGCGTCAGCTACGACAAACCa              >  2:2519516/1‑100 (MQ=255)
                                                                                       accagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCa              >  2:769013/5‑100 (MQ=255)
                                                                                           gTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTTTTTTTTTCTTCAGCTACGACAAACCAATTa          >  2:1030922/1‑100 (MQ=255)
                                                                                           gTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTa          >  1:2888888/1‑100 (MQ=255)
                                                                                             tCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc        >  2:3914909/1‑100 (MQ=255)
                                                                                             tCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc        >  1:1699340/1‑100 (MQ=255)
                                                                                              ccGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCg       <  2:2175978/100‑1 (MQ=255)
                                                                                               cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGt      <  1:3419156/100‑1 (MQ=255)
                                                                                               cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGt      >  1:3141142/1‑100 (MQ=255)
                                                                                                   gACGGAACTGGTGCGCTAAGGCGAACTACAAACGAACGATGCACTGACTCCGATATTCTTTGGCAGCGGAGTCAGCTACGACAAACCAATTATCGTCatc  <  2:2442888/100‑3 (MQ=255)
                                                                                                   |                                                                                                   
TATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGC  >  NC_000913/2653006‑2653204

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: