Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,452 G→T 16.5% intergenic (‑96/+20) cutA ← / ← dcuA divalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,4520GT16.5% 670.9 / 98.9 250intergenic (‑96/+20)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base G (97/111);  new base T (22/20);  total (119/131)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

GCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAT  >  NC_000913/4365353‑4365551
                                                                                                   |                                                                                                   
   tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                  <  2:151454/100‑1 (MQ=255)
   tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                  <  1:1849445/100‑1 (MQ=255)
     aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGAACGc                                                                                                >  2:4309856/1‑100 (MQ=255)
       cacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGccc                                                                                              <  1:1526997/100‑1 (MQ=255)
         ccTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcg                                                                                            >  2:2745129/1‑100 (MQ=255)
          cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGGCGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                           <  2:326497/100‑1 (MQ=255)
          cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                           >  1:3101668/1‑100 (MQ=255)
          cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                           >  1:3521767/1‑100 (MQ=255)
            cGCAAGAAGAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTGGTCGGCGGGAAAGAACGCGGGGCGGTCACGccc                                                                                              <  1:2519880/95‑1 (MQ=255)
            cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGa                                                                                         >  2:2382949/1‑100 (MQ=255)
             gCAAGAAGAGATGGAACCGCAAGAGAGAGGTCACTTTACAACGGACGGGGGAAGGATGTTAGTCGGTGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAt                                                                                        <  1:2910055/100‑1 (MQ=255)
               aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                      <  1:1729534/100‑6 (MQ=255)
               aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                      <  1:4377081/100‑6 (MQ=255)
                aGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                     <  1:805773/100‑1 (MQ=255)
                aGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                     <  2:559963/100‑1 (MQ=255)
                    cAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                 >  1:28029/1‑100 (MQ=255)
                    cAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                 >  2:104010/1‑100 (MQ=255)
                        tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                             >  2:1936385/1‑100 (MQ=255)
                         ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGTGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                            >  1:2351706/1‑100 (MQ=255)
                         ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                            <  1:4059027/100‑1 (MQ=255)
                         ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                            <  1:775588/100‑1 (MQ=255)
                         ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                            <  1:3714053/100‑1 (MQ=255)
                          gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGc                                                                           >  1:1274926/1‑100 (MQ=255)
                           aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                          <  2:3582853/100‑1 (MQ=255)
                           aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                          >  1:636272/1‑100 (MQ=255)
                           aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                          >  1:3045184/1‑100 (MQ=255)
                                cAAGAGAGAGGTCACTATACAACGGACGGGGTAAGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                     <  2:4108564/100‑1 (MQ=255)
                                cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                     <  1:3250349/100‑1 (MQ=255)
                                cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                     >  2:1789895/1‑100 (MQ=255)
                                cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                     <  2:1356395/100‑1 (MQ=255)
                                cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                     >  1:1142402/1‑100 (MQ=255)
                                    agggagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAc                                                                 <  2:1234846/97‑1 (MQ=255)
                                     gagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAcc                                                                <  2:289806/100‑1 (MQ=255)
                                      agagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAccc                                                               <  2:1697818/100‑1 (MQ=255)
                                      agagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAccc                                                               <  1:926721/100‑1 (MQ=255)
                                      agagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAccc                                                               <  2:3748649/100‑1 (MQ=255)
                                         ggTCACTATACAACGGACGGGGGCAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGATCATTACAGCATGAAGCTACCCa                                                              >  1:1212635/1‑98 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                         >  2:817802/1‑100 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                         <  2:726490/100‑1 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGAAcg                                                         >  1:1003329/1‑100 (MQ=255)
                                              cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaa                                                       >  2:3537838/1‑100 (MQ=255)
                                               tataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaag                                                      <  1:2731594/100‑1 (MQ=255)
                                                ataCAACGGACGGCAGAAGGATGTTAGTCTGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGATGCTACCCAGCAcgccg                                                      <  1:2974047/99‑4 (MQ=255)
                                                 taCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagcc                                                    <  1:381318/100‑1 (MQ=255)
                                                   cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccga                                                  >  2:896148/1‑100 (MQ=255)
                                                   cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccga                                                  >  2:1894508/1‑100 (MQ=255)
                                                    aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                 >  2:3237078/1‑100 (MQ=255)
                                                    aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCACGCGATTGATCATTAGAGCATGAAGCTACCCAGCAcgaagccgca                                                 >  1:216976/1‑98 (MQ=255)
                                                     acggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                >  2:786647/1‑100 (MQ=255)
                                                     acggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                <  1:1389308/100‑1 (MQ=255)
                                                      cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                               <  2:3290117/100‑1 (MQ=255)
                                                      cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGAGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                               <  1:1970709/100‑1 (MQ=255)
                                                       ggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                              >  2:3906479/1‑100 (MQ=255)
                                                           gggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAAGAAGCTACCCAGCACGAAGCCGCAGCAAACg                                          >  2:1332781/1‑100 (MQ=255)
                                                            ggggAAGGGTGGTAGTCGGCGGGAAAGAAAGCGGGGCGTGGGCGCGCCGCGGGTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAGCgg                                         >  1:3664914/1‑100 (MQ=255)
                                                            ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCacccgg                                         >  1:4303745/1‑95 (MQ=255)
                                                            ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                         >  2:4012862/1‑100 (MQ=255)
                                                            ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                         >  1:1661091/1‑100 (MQ=255)
                                                            ggggAACGATGTTAGCGCGCGGGCTCGAAAGCGGGGCTTGAACGCCCCCCGATTGATCATTACAGCATGTAGCTACCGATCAcgaagccgattca                                              <  1:2559275/95‑5 (MQ=255)
                                                            gaggAAGGATGTTGGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                         >  2:2989379/3‑100 (MQ=255)
                                                               gAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCa                                      >  1:2857478/1‑100 (MQ=255)
                                                                 aGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                    <  1:3908990/100‑1 (MQ=255)
                                                                   gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                  <  2:1183924/100‑1 (MQ=255)
                                                                    aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                 <  1:3201979/100‑1 (MQ=255)
                                                                     tGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCaa                                <  2:3805758/100‑1 (MQ=255)
                                                                     tGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                 <  1:2134328/99‑1 (MQ=255)
                                                                       ttAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAaca                              >  1:357506/1‑100 (MQ=255)
                                                                       ttAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCACGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAaca                              <  2:1094393/100‑1 (MQ=255)
                                                                          gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                           <  1:499741/100‑1 (MQ=255)
                                                                            cggcggGAAAGAAAGCGGGGCGTGATCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                         <  2:3404684/100‑1 (MQ=255)
                                                                            cggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                         >  2:3196530/1‑100 (MQ=255)
                                                                              gcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagan                       <  1:979690/100‑2 (MQ=255)
                                                                              gcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                       <  1:3371929/100‑1 (MQ=255)
                                                                               cggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                      >  2:1989100/1‑100 (MQ=255)
                                                                                 ggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAc                    >  2:2090073/1‑100 (MQ=255)
                                                                                  gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAcc                   >  2:3855843/1‑100 (MQ=255)
                                                                                  gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAcc                   >  1:410703/1‑100 (MQ=255)
                                                                                    tcaatcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                 >  1:1539958/7‑100 (MQ=39)
                                                                                    aagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                 <  1:580038/100‑1 (MQ=255)
                                                                                    aagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                 >  1:480986/1‑100 (MQ=255)
                                                                                      gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg               >  1:4393844/1‑100 (MQ=255)
                                                                                      gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg               <  2:3135318/100‑1 (MQ=255)
                                                                                      gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg               >  2:1672217/1‑100 (MQ=255)
                                                                                      aatcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg               >  2:249825/5‑100 (MQ=255)
                                                                                       atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa              >  1:2465952/4‑100 (MQ=255)
                                                                                       atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa              >  2:2340315/4‑100 (MQ=255)
                                                                                       aaagCGGGGCGTGAATGCCGCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCACCACACAGAGTACCAggca              >  1:854755/1‑98 (MQ=255)
                                                                                        tcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGGCAGGGCAACACCCAGAGTACACGGGAt             >  2:253698/3‑100 (MQ=255)
                                                                                        aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt             >  1:1433277/1‑100 (MQ=255)
                                                                                         cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg            >  1:313808/2‑100 (MQ=255)
                                                                                         cgCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg            >  2:3506823/2‑100 (MQ=255)
                                                                                         agCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg            <  1:340531/100‑1 (MQ=255)
                                                                                          gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga           >  2:3499854/1‑100 (MQ=255)
                                                                                          gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga           >  2:1599571/1‑100 (MQ=255)
                                                                                          gCGGGGCGTTCACGCCCCGCGATTGAACATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga           >  2:2682832/1‑100 (MQ=255)
                                                                                          gCGGGGCATGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg               <  1:2612599/96‑1 (MQ=255)
                                                                                            ggggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaag         >  1:1632205/1‑100 (MQ=255)
                                                                                            ggggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaag         <  2:449209/100‑1 (MQ=255)
                                                                                               gCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc      <  2:3936732/100‑1 (MQ=255)
                                                                                               gCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc      <  2:3432214/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAT  >  NC_000913/4365353‑4365551

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: