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

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

Only 100 of 197 total aligned reads displayed.

TTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAG  >  NC_000913/2653005‑2653203
                                                                                                   |                                                                                                   
 tattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg                                                                                                    <  1:2337472/100‑1 (MQ=255)
     gATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACgg                                                                                                <  1:2726999/100‑1 (MQ=255)
        gcccgcccGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGa                                                                                               >  2:3806046/1‑98 (MQ=255)
         cccgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGaa                                                                                              >  1:2576196/1‑98 (MQ=255)
          ccgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcat                                                                                           <  2:3058485/100‑3 (MQ=255)
           cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTgg                                                                                          <  2:3802198/100‑1 (MQ=255)
           cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTgg                                                                                          <  1:2064341/100‑1 (MQ=255)
           cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcatt                                                                                          <  1:1611686/100‑4 (MQ=255)
             cccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTg                                                                                        <  1:3652187/100‑1 (MQ=255)
              ccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcag                                                                                       <  1:3600193/100‑7 (MQ=255)
              ccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattca                                                                                        <  2:3900557/99‑6 (MQ=255)
               cgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcg                                                                                      <  1:3213183/100‑1 (MQ=255)
               cgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcagt                                                                                      <  2:1760415/100‑8 (MQ=255)
                ggCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcgc                                                                                     <  1:2360953/100‑1 (MQ=255)
                ggCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcagtg                                                                                     <  1:2299829/100‑9 (MQ=255)
                 gCTGCACGTTTTAACGCAGAAGTTGATGATCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcag                                                                                       <  2:1953470/97‑7 (MQ=255)
                 gCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGcgcg                                                                                    <  1:246982/100‑1 (MQ=255)
                  cTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGa                                                                                   >  1:1076202/1‑100 (MQ=255)
                    gCACGTTTTAACGCAGAAGTGGACGAGCCTCGGCCAGGTTTACGTCGCGGACATATTCCCTGAGCACTGAATGTTCCGTCCACGGAAcattcag                                                                                       <  2:2512527/94‑7 (MQ=255)
                      aCGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGc                                                                               <  1:3540163/100‑1 (MQ=255)
                        gTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGa                                                                             >  2:32150/1‑100 (MQ=255)
                          tttAACGCAGAAGTTGATGAACCTCGCTCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGaan                                                                           <  2:653072/100‑2 (MQ=255)
                          tttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAAc                                                                           <  1:1020033/100‑1 (MQ=255)
                           ttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt                                                                          >  2:1319156/1‑100 (MQ=255)
                            tAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTa                                                                         >  2:1474883/1‑100 (MQ=255)
                            tAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTa                                                                         >  1:707832/1‑100 (MQ=255)
                                gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                     <  1:1995148/100‑1 (MQ=255)
                                   gAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGa                                                                  <  1:3185255/100‑1 (MQ=255)
                                       ttGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGa                                                              <  2:2122843/100‑1 (MQ=255)
                                         gatgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATg                                                            >  1:348442/1‑100 (MQ=255)
                                          atgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGa                                                           >  2:2404205/1‑100 (MQ=255)
                                              agcTCGCCCAGGTTTACGGCGCGGAGATATTCCAGGAGCACTGAATTTTACGTGGAGAGAACTGGGGCGCGAATGCGAACGAAACGCGAcctgtga                                                           >  2:4197065/3‑91 (MQ=255)
                                              aCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTg                                                       >  1:3517103/1‑100 (MQ=255)
                                              aCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTg                                                       >  2:3714166/1‑100 (MQ=255)
                                               ccTCGACCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAAc                                                         >  2:3331717/1‑97 (MQ=255)
                                                 tCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                    <  2:2408812/100‑1 (MQ=255)
                                                   gCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTTCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGc                                                  >  2:1065228/1‑100 (MQ=255)
                                                   gCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTAGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGc                                                  <  2:353353/100‑1 (MQ=255)
                                                    cccAGGTTTACGTCGCGGACATATTCCCGGAGCCCTGCATGTTTCGTGGACGGAACTGGTGCGGGAAGGCGAAATAAAAAGGACCGATGAACTGGATGc                                                  >  2:821529/1‑99 (MQ=255)
                                                    cccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAAATGGATGCg                                                 >  2:729713/1‑100 (MQ=255)
                                                    cccAGGTTTACGTCGCGGACATATTCCCCCAGCGCTGCATGTTCCCTGGACGGAACTGGTGCCCGAAGCCGAACTAAAAACGACCCAGGAACTGGATGCg                                                 >  2:229400/1‑100 (MQ=255)
                                                     ccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAATACGACCGATGAACTGGATGCg                                                 <  2:3601672/99‑1 (MQ=255)
                                                      ttccgatctgTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGat                                               <  2:2197910/91‑1 (MQ=255)
                                                      cAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGat                                               <  1:1146570/100‑1 (MQ=255)
                                                      cAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCACGACGGCGAACTAAAACCGACCGATGAACTGGATGCGat                                               <  2:2163772/100‑1 (MQ=255)
                                                      cAGGTCTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACt                                                        <  2:2977620/91‑1 (MQ=255)
                                                       aGGTTTACGTCGCGGACATATTCCCGGAGAACTGACTGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGat                                               <  2:1447801/99‑1 (MQ=255)
                                                          tttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttt                                           <  1:1679907/100‑1 (MQ=255)
                                                            tACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                         <  2:926052/100‑1 (MQ=255)
                                                             aCGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                        >  1:871876/1‑100 (MQ=255)
                                                             aCGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                        <  1:1432262/100‑1 (MQ=255)
                                                              cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                       <  1:3086643/100‑1 (MQ=255)
                                                              cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                       <  1:4071868/100‑1 (MQ=255)
                                                              cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                       <  1:3277650/100‑1 (MQ=255)
                                                               gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGatagatcggaa                                      >  1:2197914/1‑92 (MQ=255)
                                                               gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTgg                                      >  2:3749070/1‑100 (MQ=255)
                                                               gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTgg                                      <  1:4104860/100‑1 (MQ=255)
                                                               gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTgg                                      >  2:2537472/1‑100 (MQ=255)
                                                                tCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                       >  2:3842447/1‑98 (MQ=255)
                                                                tCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                       <  2:453994/98‑1 (MQ=255)
                                                                 cgcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGt                                     >  2:1650011/1‑99 (MQ=255)
                                                                 cgcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTc                                    >  2:313855/1‑100 (MQ=255)
                                                                  gcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTcg                                   >  2:2593452/1‑100 (MQ=255)
                                                                  gcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTcg                                   >  2:3506064/1‑100 (MQ=255)
                                                                   cgGATATATTCCCGGACCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACAGATGAACTGGATGCGATATTTTTTGGTcgc                                  >  2:2448697/1‑100 (MQ=255)
                                                                   cgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTcgc                                  >  1:3864950/1‑100 (MQ=255)
                                                                   cgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTcgc                                  >  2:1804805/1‑100 (MQ=255)
                                                                      aCATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggc                               <  1:1362559/100‑1 (MQ=255)
                                                                         tatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACTGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTc                            >  1:1893567/1‑100 (MQ=255)
                                                                         tatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGTTCGcgtcgt                             >  2:457047/1‑99 (MQ=255)
                                                                         tatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTc                            >  1:4050227/1‑100 (MQ=255)
                                                                          atTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                           >  1:3292173/1‑100 (MQ=255)
                                                                          atTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                           <  1:450048/100‑1 (MQ=255)
                                                                          atTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                           >  1:659041/1‑100 (MQ=255)
                                                                           ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAg                          <  1:2927629/100‑1 (MQ=255)
                                                                                 gAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGa                    <  2:2931842/100‑1 (MQ=255)
                                                                                  aGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAc                   <  1:264529/100‑1 (MQ=255)
                                                                                  aGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAc                   <  1:1326648/100‑1 (MQ=255)
                                                                                   gCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACa                  >  2:4167034/1‑100 (MQ=255)
                                                                                    cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaa                 >  1:3696834/1‑100 (MQ=255)
                                                                                    cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACa                  >  2:2915528/1‑99 (MQ=255)
                                                                                     cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                >  1:3006592/8‑100 (MQ=255)
                                                                                     cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                >  1:1274177/8‑100 (MQ=255)
                                                                                     cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                >  2:1164676/8‑100 (MQ=255)
                                                                                     cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                >  1:2053191/8‑100 (MQ=255)
                                                                                     aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                >  2:2786381/1‑100 (MQ=255)
                                                                                       tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAcc              >  2:3018555/1‑100 (MQ=255)
                                                                                       caccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAcc              >  1:1016560/6‑100 (MQ=255)
                                                                                        gAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCa             >  1:3920042/1‑100 (MQ=255)
                                                                                        accagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCa             >  1:545230/5‑100 (MQ=255)
                                                                                        accagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCa             >  2:769013/5‑100 (MQ=255)
                                                                                           tGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAAtt          >  2:2655765/1‑100 (MQ=255)
                                                                                            gTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTTTTTTTTTCTTCAGCTACGACAAACCAATTa         >  2:1030922/1‑100 (MQ=255)
                                                                                            gTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTa         <  1:3730472/100‑1 (MQ=255)
                                                                                            gTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTa         >  1:2888888/1‑100 (MQ=255)
                                                                                              tCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc       >  1:1699340/1‑100 (MQ=255)
                                                                                               ccGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGCCAGCTACGACAAACCAATTATCg      <  2:1697749/100‑1 (MQ=255)
                                                                                                cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGt     <  1:3419156/100‑1 (MQ=255)
                                                                                                cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGt     >  1:3141142/1‑100 (MQ=255)
                                                                                                 gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc    <  1:378170/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAG  >  NC_000913/2653005‑2653203

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: