Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 18.0% intergenic (‑34/+25) 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,4450AT18.0% 422.4 / 85.9 179intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (70/76);  new base T (20/12);  total (91/88)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

Only 100 of 186 total aligned reads displayed.

GAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATC  >  NC_000913/4356346‑4356544
                                                                                                   |                                                                                                   
gAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                     >  1:752627/1‑100 (MQ=255)
 aaGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCaa                                                                                                    >  1:3716844/1‑100 (MQ=255)
   gTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGc                                                                                                  <  1:734229/100‑1 (MQ=255)
    tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                 <  1:3477761/100‑1 (MQ=255)
    tACTCCCAGATTTGGATCGCCACGATATAGTATAGCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                 <  2:2054856/100‑1 (MQ=255)
    tACTCCCAGATTTGGATAGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTACAAAAAAAAAAAAAAGGGAAGTGGCAAGcc                                                                                                 <  2:1840475/100‑1 (MQ=255)
         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  1:2383276/100‑1 (MQ=255)
         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  1:26861/100‑1 (MQ=255)
         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  1:2177573/100‑1 (MQ=255)
         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  2:3088876/100‑1 (MQ=255)
         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  2:3180632/100‑1 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                          <  2:194522/100‑1 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                          <  1:2609535/100‑1 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                          <  1:3023198/100‑1 (MQ=255)
            gATTTGGATCGCCACGTTATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                          <  2:1639492/99‑1 (MQ=255)
             aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCATTTCCCtt                                                                                        <  1:2073502/100‑1 (MQ=255)
             aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                        <  1:3923687/100‑1 (MQ=255)
              tttGGATCGCCACGATATAGTATATCTCGCGCGGCTGTGAGGGTGTTTTCATCTGTTCTCCTCATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTg                                                                                       <  2:1157155/100‑1 (MQ=255)
              tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGGGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCttt                                                                                       <  1:3769322/100‑2 (MQ=255)
                 ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttttg                                                                                    <  2:1251985/100‑5 (MQ=38)
                 ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGAAATAAAGGCAAGTTGCAAGCCACTTCCCTTGTAc                                                                                    <  2:3427610/100‑1 (MQ=255)
                  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)
                      gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGGTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                                >  1:252049/1‑99 (MQ=255)
                      gCCACGATATAGTATATCGCGCGCGGCTGTGAGAGAGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCATGTACGAGCt                                                                               >  2:2607845/1‑100 (MQ=255)
                        cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                                <  2:489214/97‑1 (MQ=255)
                           gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                          <  1:991890/100‑1 (MQ=255)
                           gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                          >  1:3845962/1‑100 (MQ=255)
                               taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                      <  1:86623/100‑1 (MQ=255)
                                aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                     <  1:1274643/100‑1 (MQ=255)
                                 gTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                    <  2:2288551/100‑1 (MQ=255)
                                 gTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                    <  2:1049823/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)
                                        gcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                             <  2:2918664/100‑1 (MQ=255)
                                              tccgatctgggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  2:2927749/92‑1 (MQ=255)
                                              ggatgtgAGGGTGTTTTCATGTGTGCTCCTTATGAGCAAAAAAGGGATGTGGCAAGCCACCGCCCTTGTACGCGCTAATTATTTTTTGCTAtcttct                                                          <  2:2220772/94‑1 (MQ=255)
                                              ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  2:2713065/100‑1 (MQ=255)
                                               gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                      >  1:20064/1‑100 (MQ=255)
                                                   tgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                  >  1:2969462/1‑100 (MQ=255)
                                                   tgAGGGTGTTTTCATGGGTTCGCCGTATGAGCAAAAAAGGGAAGTGGCAAGCTACTTCCCTTGTAAGAGCGCATTATTTTTTTCTTTATTCTTTTAATAc                                                  >  2:2880628/1‑100 (MQ=255)
                                                    gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTCTTTTTTTCTTTCAATAcc                                                 >  2:598254/1‑100 (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                                            <  1:856058/100‑1 (MQ=255)
                                                          gtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                           >  1:3522038/1‑100 (MQ=255)
                                                          gtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                           >  2:1267569/1‑100 (MQ=255)
                                                           ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                          >  1:1236291/1‑100 (MQ=255)
                                                              tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTCTTTTTTCaaaaactcacg                                            >  2:3052/1‑86 (MQ=255)
                                                                 tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                    <  2:3522013/100‑1 (MQ=255)
                                                                 tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                    >  1:3565443/1‑100 (MQ=255)
                                                                 tgtgTTCTCATTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                    >  1:1087722/1‑100 (MQ=255)
                                                                    gtTGTCCTTATGCGCAAAAAAGGGAAGTGTCAAGCCACTTCCCTTGGACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                  <  2:3421006/99‑1 (MQ=255)
                                                                     ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                <  1:472050/100‑1 (MQ=255)
                                                                      tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCa                               >  1:2235558/1‑100 (MQ=255)
                                                                      tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCa                               >  1:64105/1‑100 (MQ=255)
                                                                        tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATc                             <  2:717704/100‑1 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGTCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            <  1:2097033/100‑1 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            >  1:3700975/1‑100 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            >  2:3077906/1‑100 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            >  2:2797374/1‑100 (MQ=255)
                                                                           ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                          <  1:3388402/100‑1 (MQ=255)
                                                                           ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                          <  1:1669945/100‑1 (MQ=255)
                                                                            tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGGGGCCATCa                            >  2:3230363/1‑97 (MQ=255)
                                                                            tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                         <  2:1370958/100‑1 (MQ=255)
                                                                              tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                        >  1:1023302/1‑99 (MQ=255)
                                                                              tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                       >  1:231920/1‑100 (MQ=255)
                                                                              tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                       >  2:1320408/1‑100 (MQ=255)
                                                                                 gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAATGGTATAGCGGCCATCAGCCTGAc                     >  2:2500533/1‑99 (MQ=255)
                                                                                  cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    >  2:418062/1‑99 (MQ=255)
                                                                                  cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATGATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                   <  2:1631025/100‑1 (MQ=255)
                                                                                  aaaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    <  2:3845933/98‑1 (MQ=255)
                                                                                    aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                     >  2:782797/1‑96 (MQ=255)
                                                                                     aaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                <  2:246034/100‑1 (MQ=255)
                                                                                      caaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg               >  1:2927079/2‑100 (MQ=255)
                                                                                      caaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg               >  2:3700057/2‑100 (MQ=255)
                                                                                       aaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc              >  1:1547393/1‑100 (MQ=255)
                                                                                        aGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCa             >  1:1448466/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)
                                                                                         gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc            >  1:345533/1‑100 (MQ=255)
                                                                                         gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATGATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc            <  1:1842464/100‑1 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTGTTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:2783523/1‑100 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCATTACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  2:678282/1‑100 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:3784199/1‑100 (MQ=255)
                                                                                          ggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:1941323/1‑100 (MQ=255)
                                                                                           gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg          >  1:3353864/1‑100 (MQ=255)
                                                                                           gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg          <  1:1813749/100‑1 (MQ=255)
                                                                                              gTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTg        >  1:962126/1‑99 (MQ=255)
                                                                                               tGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa      <  1:3277107/100‑1 (MQ=255)
                                                                                                 gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAata    <  1:2864235/100‑1 (MQ=255)
                                                                                                   aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATc  >  2:64944/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATC  >  NC_000913/4356346‑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: