Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,452 G→T 15.9% 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,4520GT15.9% 1020.9 / 142.7 381intergenic (‑96/+20)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base G (177/141);  new base T (35/26);  total (212/169)
Fisher's exact test for biased strand distribution p-value = 8.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 387 total aligned reads displayed.

GCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAT  >  NC_000913/4365353‑4365551
                                                                                                   |                                                                                                   
gCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTg                                                                                                     <  1:1380678/100‑1 (MQ=255)
   tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                   >  1:1283155/1‑99 (MQ=255)
   tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                  <  2:2230433/100‑1 (MQ=255)
     aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACg                                                                                                 <  1:2024618/99‑1 (MQ=255)
      acacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGcc                                                                                               >  1:2557674/1‑100 (MQ=255)
       cacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGcc                                                                                               >  1:1672318/1‑99 (MQ=255)
       cacCTCGCAAGAACAGATGGAACCGAAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGccc                                                                                              <  1:2084168/100‑1 (MQ=255)
         ccTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcg                                                                                            >  2:767873/1‑100 (MQ=255)
         ccTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcg                                                                                            >  2:2470408/1‑100 (MQ=255)
           tCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgct                                                                                          <  1:1829647/100‑2 (MQ=255)
           tCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                          >  2:128495/1‑100 (MQ=255)
             gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                        <  1:786667/100‑4 (MQ=255)
             gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                        <  1:922801/100‑4 (MQ=255)
             gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAt                                                                                        <  2:2542538/100‑1 (MQ=255)
              cAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                       <  1:2877642/100‑1 (MQ=255)
              cAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                       <  1:459183/100‑1 (MQ=255)
              cAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGCCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGa                                                                                         >  1:1918402/1‑98 (MQ=255)
               aaGAGCAGAGGGAGCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                       <  1:1986689/99‑1 (MQ=255)
               aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                      <  1:1230917/100‑6 (MQ=255)
               aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                      <  2:1714725/100‑6 (MQ=255)
               aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                      <  2:226535/100‑1 (MQ=255)
               aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                      <  2:395582/100‑1 (MQ=255)
                aGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                     >  2:473470/1‑100 (MQ=255)
                 gAACAGATGGAAGCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGGGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                     <  1:101775/99‑1 (MQ=255)
                  naCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                   >  1:2341599/2‑100 (MQ=255)
                  aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                   <  1:2135680/100‑1 (MQ=255)
                  aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                   <  2:139352/100‑1 (MQ=255)
                    cAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                 <  2:1660838/100‑1 (MQ=255)
                      gATGGAGACGCAGGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                 <  1:2432527/98‑1 (MQ=255)
                      gATGGAACCTCAAGAGAGAGGTCACTATACAACGGACGGGGAAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                               >  1:1587533/1‑100 (MQ=255)
                      gATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAACGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                               >  2:1173065/1‑100 (MQ=255)
                        tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGACAGAAAACGGGGCGTGAACGCCCCGCGATTGATGATTACa                                                                             >  1:1766892/1‑100 (MQ=255)
                        tGGAACCGCAACAGAGAGGTCACTATACAACGGACGGGGGCAGGATGTTAGTCGGCGGAAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGATCATTACa                                                                             >  1:688592/1‑100 (MQ=255)
                            aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAt                                                                         >  2:2176791/1‑100 (MQ=255)
                            aCCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                          <  1:117196/99‑1 (MQ=255)
                             ccGCAAGAGAGAGGTCACTATACAACGGACGGGGGGAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAt                                                                         >  1:2202143/1‑99 (MQ=255)
                             ccGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAt                                                                         <  2:2330824/99‑1 (MQ=255)
                               gCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                      >  2:1221180/1‑100 (MQ=255)
                               gCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                      >  1:875138/1‑100 (MQ=255)
                                cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                     <  2:1895171/100‑1 (MQ=255)
                                 aaGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                    <  2:273594/100‑1 (MQ=255)
                                 aaGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                    >  1:1619512/1‑100 (MQ=255)
                                 aaGAGAGAGGTCACTATACAACGGACGGGCGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                    <  2:87973/100‑1 (MQ=255)
                                  agagagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                    <  1:2381294/99‑1 (MQ=255)
                                       gagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCa                                                              >  2:1211386/1‑100 (MQ=255)
                                         ggTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                            >  1:1106505/1‑100 (MQ=255)
                                         ggTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                            >  1:1570778/1‑100 (MQ=255)
                                         ggTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                            <  2:1557505/100‑1 (MQ=255)
                                         ggTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                            <  2:357603/100‑1 (MQ=255)
                                         ggTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCACCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                            >  1:2586610/1‑100 (MQ=255)
                                          gTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTCATCATTACAGCATGAAGCTACccacca                                                           >  1:247/1‑100 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                         >  1:1205836/1‑100 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                         <  2:1001838/100‑1 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                         >  2:1049707/1‑100 (MQ=255)
                                            cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                         >  1:495834/1‑100 (MQ=255)
                                               tataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaag                                                      >  2:2022502/1‑100 (MQ=255)
                                               tataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCAGCGATTGAACATTACAGCATGAAGCTACCCAGCAcgaag                                                      >  1:1658038/1‑100 (MQ=255)
                                                 taCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagcc                                                    <  1:1478227/100‑1 (MQ=255)
                                                   cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACACAGCAcgaagccga                                                  >  2:851191/1‑100 (MQ=255)
                                                     acggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                <  2:1130666/100‑1 (MQ=255)
                                                      cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                               <  1:1032019/100‑1 (MQ=255)
                                                      cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                               >  1:106176/1‑100 (MQ=255)
                                                        gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCaa                                             >  1:2278865/1‑100 (MQ=255)
                                                          cggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAAc                                           >  2:539319/1‑100 (MQ=255)
                                                            gggtAAGGATGTCAGACGGCGAGCAAGAAACCGAGGCGTGAACGCCCCGCGCTGGATCATTACAGCATGAAGCTACCCAGCACGACGCCGAAGCAAAcag                                         >  1:2661149/1‑98 (MQ=255)
                                                            ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagctccgat                                         >  2:1023578/1‑94 (MQ=255)
                                                            ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                          >  1:6418/1‑99 (MQ=255)
                                                               gAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCa                                      >  2:2361411/1‑100 (MQ=255)
                                                                 aGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                    >  2:1643690/1‑100 (MQ=255)
                                                                  ggATGTGAGTCGGCGGGAAAGAAAGCGGGGCGTCCACGCCCCGAGATTGATCATTACAGCATGAAGCTACCCAGCACGAGGCCGAAGCAAACGGCCAg                                     <  1:533764/98‑1 (MQ=255)
                                                                   gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACGGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgggn                                  <  2:281633/100‑2 (MQ=255)
                                                                    aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                 >  1:413602/1‑100 (MQ=255)
                                                                        tAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAacac                             >  1:2846648/1‑100 (MQ=255)
                                                                          gTCGGCGGGAAAAAAAGCGGGGCGTGAACTCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCCCGAAGCCGAAGCAAACGGACAGGGCAACAccc                           >  1:2782684/1‑100 (MQ=255)
                                                                           tCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                          <  2:1598457/100‑1 (MQ=255)
                                                                           tCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                          <  1:2150097/100‑1 (MQ=255)
                                                                            cggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                         <  2:2903145/100‑1 (MQ=255)
                                                                              gcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGTAAACGGCCAGGGCAACACCCagag                       <  1:369658/100‑1 (MQ=255)
                                                                               cggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACTCAGAGt                      >  2:2187165/1‑100 (MQ=255)
                                                                               cggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                      >  1:325496/1‑100 (MQ=255)
                                                                                 ggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAc                    >  2:2851762/1‑100 (MQ=255)
                                                                                  gatcaatcacGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAACCGAAGCAAACGACCAGGGCAACACCCAGAATa                     >  1:837523/10‑98 (MQ=37)
                                                                                   atcaatcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAccc                  >  2:585054/8‑100 (MQ=38)
                                                                                    aagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                 <  1:174391/100‑1 (MQ=255)
                                                                                     agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCgg                <  2:1788778/100‑1 (MQ=255)
                                                                                     agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCgg                <  2:2442611/100‑1 (MQ=255)
                                                                                     agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                 >  2:1996246/1‑99 (MQ=255)
                                                                                      gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg               >  2:2662262/1‑100 (MQ=255)
                                                                                      aatcgCGGAGCGTTCACGCCAAGCGATTGATCATTAAAGCATGAAGCTAGCCCGCACGAAGCTGAATCAACCTGCCAGGGCAACACCCAGAGTACCCg                 >  1:2743605/5‑98 (MQ=255)
                                                                                       atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa              >  2:916577/4‑100 (MQ=255)
                                                                                       atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa              >  1:1088420/4‑100 (MQ=255)
                                                                                        aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt             >  1:754767/1‑100 (MQ=255)
                                                                                        aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt             <  2:1090142/100‑1 (MQ=255)
                                                                                         cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg            >  2:146417/2‑100 (MQ=255)
                                                                                         cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg            <  2:2621349/99‑1 (MQ=255)
                                                                                          gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga           >  2:1797716/1‑100 (MQ=255)
                                                                                             gggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaaga        <  1:2424534/100‑1 (MQ=255)
                                                                                             gggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaaga        >  2:2081329/1‑100 (MQ=255)
                                                                                              ggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa       >  2:498558/1‑100 (MQ=255)
                                                                                                cGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg     <  1:1078186/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: