Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,446 A→T 17.9% intergenic (‑35/+24) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4460AT17.9% 422.3 / 83.9 179intergenic (‑35/+24)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (71/76);  new base T (20/12);  total (91/88)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

Only 100 of 184 total aligned reads displayed.

AAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATC  >  NC_000913/4356347‑4356544
                                                                                                   |                                                                                                  
aaGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCaa                                                                                                    >  1:3716844/1‑100 (MQ=255)
  gTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGc                                                                                                  <  1:734229/100‑1 (MQ=255)
    aCTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCa                                                                                                >  1:894679/1‑100 (MQ=255)
       cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                             >  1:2219633/1‑100 (MQ=255)
       cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                             >  2:2176410/1‑100 (MQ=255)
        ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            >  2:1373203/1‑100 (MQ=255)
        ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  2:3180632/100‑1 (MQ=255)
        ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  2:3088876/100‑1 (MQ=255)
        ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  1:26861/100‑1 (MQ=255)
        ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  1:2177573/100‑1 (MQ=255)
         cAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTcc                                                                                           <  1:3824963/100‑1 (MQ=255)
          aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                          <  1:2609535/100‑1 (MQ=255)
          aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                          <  1:3023198/100‑1 (MQ=255)
           gATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCt                                                                                         >  1:284958/1‑100 (MQ=255)
            aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCATTTCCCtt                                                                                        <  1:2073502/100‑1 (MQ=255)
            aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                        <  1:3923687/100‑1 (MQ=255)
             tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTg                                                                                       >  2:1105141/1‑100 (MQ=255)
             tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGGGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCttt                                                                                       <  1:3769322/100‑2 (MQ=255)
                ggATCGCCACGATATGGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttttt                                                                                    <  2:2547751/100‑5 (MQ=38)
                 gATCGCCACGATATATTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCCCCTTATGAGCAACAAAGGCACGTGGCAAGCCACTTCCCTTGTACg                                                                                   <  2:792624/100‑1 (MQ=255)
                  aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGTTTGCCACTTCCCttttttgc                                                                                  <  2:2410403/100‑7 (MQ=38)
                  aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgc                                                                                  <  1:100657/100‑7 (MQ=38)
                   tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAg                                                                                 <  1:2987825/100‑1 (MQ=255)
                     gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                               >  2:3770959/1‑100 (MQ=255)
                     gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                               >  1:1158524/1‑100 (MQ=255)
                     gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                               >  1:3202719/1‑100 (MQ=255)
                     gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGGTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                                >  1:252049/1‑99 (MQ=255)
                          gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                          >  1:3845962/1‑100 (MQ=255)
                           atatAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                          >  2:3511581/1‑99 (MQ=255)
                               aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                     <  1:1274643/100‑1 (MQ=255)
                                gTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                    >  1:859149/1‑100 (MQ=255)
                                gTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                    <  2:2288551/100‑1 (MQ=255)
                                 tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTg                                                                   <  1:3642941/100‑1 (MQ=255)
                                    atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                                >  2:3666994/1‑100 (MQ=255)
                                     tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCttt                                                               <  1:1767947/100‑1 (MQ=255)
                                     tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCttt                                                               >  2:2970759/1‑100 (MQ=255)
                                      cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttc                                                              >  1:1818094/1‑100 (MQ=255)
                                       gcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                             <  2:2918664/100‑1 (MQ=255)
                                           gcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                         <  2:1343477/100‑1 (MQ=255)
                                            cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                        <  1:1291249/100‑1 (MQ=255)
                                             gtcTGTGAGGGTGTTTTCATGTGTTCTCATAATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  2:354174/98‑1 (MQ=255)
                                             ggatgtgAGGGTGTTTTCATGTGTGCTCCTTATGAGCAAAAAAGGGATGTGGCAAGCCACCGCCCTTGTACGCGCTAATTATTTTTTGCTAtcttct                                                          <  2:2220772/94‑1 (MQ=255)
                                             ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  1:2350451/100‑1 (MQ=255)
                                              gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                      >  1:20064/1‑100 (MQ=255)
                                                tgtgATGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAATGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAAt                                                    <  2:3056325/100‑1 (MQ=255)
                                                tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAAt                                                    <  2:1048906/100‑1 (MQ=255)
                                                  tgAGGGTGTTTTCATGGGTTCGCCGTATGAGCAAAAAAGGGAAGTGGCAAGCTACTTCCCTTGTAAGAGCGCATTATTTTTTTCTTTATTCTTTTAATAc                                                  >  2:2880628/1‑100 (MQ=255)
                                                   gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTCTTTTTTTCTTTCAATAcc                                                 >  2:598254/1‑100 (MQ=255)
                                                     gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCagatcgga                                               >  1:2927767/1‑93 (MQ=255)
                                                     gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                               >  1:1459518/1‑100 (MQ=255)
                                                     gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                               >  1:3166505/1‑100 (MQ=255)
                                                       gtgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                             >  1:749078/1‑100 (MQ=255)
                                                       gtgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTa                                              <  2:1122349/99‑1 (MQ=255)
                                                        tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                            <  2:568498/100‑1 (MQ=255)
                                                        tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                            <  2:1169907/100‑1 (MQ=255)
                                                          ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                          >  1:1236291/1‑100 (MQ=255)
                                                                tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                    >  1:3565443/1‑100 (MQ=255)
                                                                tgtgTTCTCATTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                    >  1:1087722/1‑100 (MQ=255)
                                                                  tgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                  <  1:353580/100‑1 (MQ=255)
                                                                   gtTGTCCTTATGCGCAAAAAAGGGAAGTGTCAAGCCACTTCCCTTGGACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                  <  2:3421006/99‑1 (MQ=255)
                                                                   gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                 >  1:3346486/1‑100 (MQ=255)
                                                                    ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                <  1:472050/100‑1 (MQ=255)
                                                                    ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                <  2:1262101/100‑1 (MQ=255)
                                                                     tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCa                               >  1:64105/1‑100 (MQ=255)
                                                                       tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATc                             <  2:717704/100‑1 (MQ=255)
                                                                        ccTTATGAGCAAAAAAGGGAAGTGTCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            <  1:2097033/100‑1 (MQ=255)
                                                                        ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            >  2:2797374/1‑100 (MQ=255)
                                                                        ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            <  1:2315453/100‑1 (MQ=255)
                                                                        ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            >  1:3700975/1‑100 (MQ=255)
                                                                          ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                          <  1:1669945/100‑1 (MQ=255)
                                                                          ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                          <  1:3388402/100‑1 (MQ=255)
                                                                          ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                          <  1:2920615/100‑1 (MQ=255)
                                                                           tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGGGGCCATCa                            >  2:3230363/1‑97 (MQ=255)
                                                                             tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                        >  1:1023302/1‑99 (MQ=255)
                                                                             tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                       >  2:2727314/1‑100 (MQ=255)
                                                                                 cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTTTTTTAATACCTTAACTGTATAGTTGCCATCatttt                        >  2:1920671/1‑91 (MQ=255)
                                                                                 cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    >  2:418062/1‑99 (MQ=255)
                                                                                 cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATGATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                   <  2:1631025/100‑1 (MQ=255)
                                                                                 aaaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    <  2:3845933/98‑1 (MQ=255)
                                                                                  gtacaaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                  >  1:3479897/5‑100 (MQ=38)
                                                                                  aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                  >  1:816256/1‑100 (MQ=255)
                                                                                  aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                  >  2:1654332/1‑100 (MQ=255)
                                                                                   aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                     >  2:782797/1‑96 (MQ=255)
                                                                                   aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                 >  2:2462964/1‑100 (MQ=255)
                                                                                     caaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg               >  2:3700057/2‑100 (MQ=255)
                                                                                     aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg               <  1:1888491/100‑1 (MQ=255)
                                                                                      aaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc              >  1:1547393/1‑100 (MQ=255)
                                                                                       aGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCa             >  1:589541/1‑100 (MQ=255)
                                                                                       aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCa             <  2:1459509/100‑1 (MQ=255)
                                                                                        gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc            >  1:10377/1‑100 (MQ=255)
                                                                                        gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc            >  1:9814/1‑100 (MQ=255)
                                                                                        gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc            >  1:345533/1‑100 (MQ=255)
                                                                                         ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTGTTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:2783523/1‑100 (MQ=255)
                                                                                         ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTTATTCAATACCTTAACGGTATAACTGCCATCATCCTTACGATATGCAcc           >  2:2052734/1‑100 (MQ=255)
                                                                                         ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCATTACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  2:678282/1‑100 (MQ=255)
                                                                                         ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCACTACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  2:2653885/1‑100 (MQ=255)
                                                                                         ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  2:1287155/1‑100 (MQ=255)
                                                                                          gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg          <  1:1813749/100‑1 (MQ=255)
                                                                                             gTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTTTTTCAATACCTTAACGGTATAGCGGTCATCAGCTTGACGGTATGCACCGTg        >  2:308862/1‑99 (MQ=255)
                                                                                                  aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATc  >  2:64944/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATC  >  NC_000913/4356347‑4356544

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: