Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 175,144 Δ1 bp 8.9% coding (38/1422 nt) clcA → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913175,1440C. 8.9% 312.6 / 24.8 90coding (38/1422 nt)clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (43/39);  new base . (5/3);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AGAGTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAAC  >  NC_000913/174996‑175292
                                                                                                                                                    |                                                                                                                                                     
agagTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGc                                                                                                                                                       <  2:1555886/149‑1 (MQ=255)
 gagTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGcc                                                                                                                                                      >  2:393306/1‑149 (MQ=255)
        aaaGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGcc                                                                                                                                              <  1:393306/149‑1 (MQ=255)
                gtgCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGc                                                                                                                                      <  3:779471/149‑1 (MQ=255)
                     tttGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGa                                                                                                                                         >  4:780510/1‑141 (MQ=255)
                     tttGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGa                                                                                                                                         <  3:780510/141‑1 (MQ=255)
                                        cAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtctt                                                                                                              >  4:84100/1‑149 (MQ=255)
                                        cAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtctt                                                                                                              <  3:84100/149‑1 (MQ=255)
                                         aaaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAtt                                                                                                                          <  1:1689169/136‑1 (MQ=255)
                                         aaaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAtt                                                                                                                          >  2:1689169/1‑136 (MQ=255)
                                          aaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAt                                                                                                                           >  3:1063937/1‑134 (MQ=255)
                                          aaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAt                                                                                                                           <  4:1063937/134‑1 (MQ=255)
                                                       gAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGAc                                                                                                                                        <  3:1782082/108‑1 (MQ=255)
                                                       gAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGAc                                                                                                                                        >  4:1782082/1‑108 (MQ=255)
                                                        aaCGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAccc                                                                                              >  2:241174/1‑149 (MQ=255)
                                                             gCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                                                                         <  1:1865941/149‑1 (MQ=255)
                                                                    aaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAttt                                                                                  >  2:529921/1‑149 (MQ=255)
                                                                    aaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACAACGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAttt                                                                                  >  4:1226269/1‑149 (MQ=255)
                                                                       tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAtt                                                                                                                          <  4:252528/106‑1 (MQ=255)
                                                                       tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAtt                                                                                                                          >  3:252528/1‑106 (MQ=255)
                                                                       tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGt                                                                               <  1:241174/149‑1 (MQ=255)
                                                                          cGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCgcgcgc                                                                                                                                               <  4:993687/82‑1 (MQ=255)
                                                                          cGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCgcgcgc                                                                                                                                               >  3:993687/1‑82 (MQ=255)
                                                                           gTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTAt                                                                           >  2:899901/1‑149 (MQ=255)
                                                                            tCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                                                                          <  3:956071/133‑1 (MQ=255)
                                                                            tCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                                                                          >  4:956071/1‑133 (MQ=255)
                                                                            tCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATg                                                                          >  3:484838/1‑149 (MQ=255)
                                                                              aTCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                                                               <  1:412305/110‑1 (MQ=255)
                                                                              aTCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                                                               >  2:412305/1‑110 (MQ=255)
                                                                                 ttATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAcc                                                                                               <  1:1261556/123‑1 (MQ=255)
                                                                                 ttATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAcc                                                                                               >  2:1261556/1‑123 (MQ=255)
                                                                                     ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGa                                                                                                                                   >  4:1518725/1‑83 (MQ=255)
                                                                                     ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGa                                                                                                                                   <  3:1518725/83‑1 (MQ=255)
                                                                                     ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGtt                                                                              <  4:1704827/136‑1 (MQ=255)
                                                                                     ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGtt                                                                              >  3:1704827/1‑136 (MQ=255)
                                                                                       gTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATaaaa                                                                                                 >  2:205939/1‑115 (MQ=255)
                                                                                       gTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATaaaa                                                                                                 <  1:205939/115‑1 (MQ=255)
                                                                                           gACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc                                                           >  3:1077635/1‑149 (MQ=255)
                                                                                                    gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATgg                                                                         <  4:1394458/126‑1 (MQ=255)
                                                                                                    gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATgg                                                                         >  3:1394458/1‑126 (MQ=255)
                                                                                                           ttGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTgg                                           >  2:229197/1‑149 (MQ=255)
                                                                                                           ttGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTgg                                           >  2:984237/1‑149 (MQ=255)
                                                                                                             ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcgg                                                                      <  1:1058848/120‑1 (MQ=255)
                                                                                                             ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcgg                                                                      >  2:1058848/1‑120 (MQ=255)
                                                                                                                tGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCg                                      <  4:484838/149‑1 (MQ=255)
                                                                                                                  aaaaCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGt                                    <  1:899901/149‑1 (MQ=255)
                                                                                                                    aaCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTggg                                               <  2:219039/136‑1 (MQ=255)
                                                                                                                    aaCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTggg                                               >  1:219039/1‑136 (MQ=255)
                                                                                                                     aCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGc                                 >  3:1818054/1‑149 (MQ=255)
                                                                                                                       tGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTGTATGGCGGCAGTCGTCTGCACGCTTGTTGGGCTGGCAGCGGTTGCtt                               >  3:1884964/1‑149 (MQ=255)
                                                                                                                         aTACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             >  1:561846/1‑149 (MQ=255)
                                                                                                                             tCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcgg                                                                      <  4:1182684/104‑1 (MQ=255)
                                                                                                                             tCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACa                         >  1:540287/1‑149 (MQ=255)
                                                                                                                             tCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACa                         >  3:714504/1‑149 (MQ=255)
                                                                                                                             tCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTCGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcgg                                                                      >  3:1182684/1‑104 (MQ=255)
                                                                                                                               ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGt                                                                               <  1:1716359/93‑1 (MQ=255)
                                                                                                                               ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGt                                                                               >  2:1716359/1‑93 (MQ=255)
                                                                                                                               ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATgg                                                                         >  4:1752676/1‑99 (MQ=255)
                                                                                                                               ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATgg                                                                         <  3:1752676/99‑1 (MQ=255)
                                                                                                                               ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACaaa                       <  1:529921/149‑1 (MQ=255)
                                                                                                                               ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACaaa                       <  3:1226269/149‑1 (MQ=255)
                                                                                                                                ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAg                      <  4:1884964/149‑1 (MQ=255)
                                                                                                                                ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAg                      <  2:561846/149‑1 (MQ=255)
                                                                                                                                 tctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAAc                                                                                                                               >  4:1653911/1‑43 (MQ=255)
                                                                                                                                 tctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAAc                                                                                                                               <  3:1653911/43‑1 (MQ=255)
                                                                                                                                   tttAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGgtg                   <  4:714504/149‑1 (MQ=255)
                                                                                                                                     tAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTc                 >  2:420209/1‑149 (MQ=255)
                                                                                                                                        aaaCACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAgg                     >  1:1594106/1‑142 (MQ=255)
                                                                                                                                        aaaCACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAgg                     <  2:1594106/142‑1 (MQ=255)
                                                                                                                                          acacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                              <  4:1473555/131‑1 (MQ=255)
                                                                                                                                          acacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                              >  3:1473555/1‑131 (MQ=255)
                                                                                                                                           cacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGa                           >  4:1015760/1‑133 (MQ=255)
                                                                                                                                           cacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGa                           <  3:1015760/133‑1 (MQ=255)
                                                                                                                                            gtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             >  2:736719/3‑130 (MQ=255)
                                                                                                                                            gtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             >  4:1448934/3‑130 (MQ=255)
                                                                                                                                            gtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             >  2:50194/3‑130 (MQ=255)
                                                                                                                                            gtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             <  1:50194/128‑1 (MQ=255)
                                                                                                                                            gtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             <  1:736719/128‑1 (MQ=255)
                                                                                                                                            gtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                             <  3:1448934/128‑1 (MQ=255)
                                                                                                                                             tcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGtt         >  1:306775/2‑149 (MQ=255)
                                                                                                                                             ccGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGtt         >  1:1873523/1‑149 (MQ=255)
                                                                                                                                              cGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTgg                                                >  4:1874668/1‑109 (MQ=255)
                                                                                                                                              cGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTgg                                                <  3:1874668/109‑1 (MQ=255)
                                                                                                                                              cGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCt                                             <  3:445619/112‑1 (MQ=255)
                                                                                                                                              cGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCt                                             >  4:445619/1‑112 (MQ=255)
                                                                                                                                              cGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTg        >  3:1191078/1‑149 (MQ=255)
                                                                                                                                               gCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGcc              >  3:1057648/1‑142 (MQ=255)
                                                                                                                                               gCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGcc              <  4:1057648/142‑1 (MQ=255)
                                                                                                                                                cAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCa      <  1:229197/149‑1 (MQ=255)
                                                                                                                                                cAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCa      <  1:984237/149‑1 (MQ=255)
                                                                                                                                                   gCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGaa   >  3:186838/1‑149 (MQ=255)
                                                                                                                                                   gCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGaa   >  3:725783/1‑149 (MQ=255)
                                                                                                                                                    ccGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAAc  <  4:1818054/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGAGTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAAC  >  NC_000913/174996‑175292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: