Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 175,144 Δ1 bp 6.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. 6.9% 392.2 / 23.5 131coding (38/1422 nt)clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (58/63);  new base . (5/4);  total (63/68)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-02

Only 100 of 133 total aligned reads displayed.

GTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACC  >  NC_000913/174999‑175293
                                                                                                                                                 |                                                                                                                                                      
gTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcgcc                                                                                                                                                    >  2:723120/1‑149 (MQ=255)
gTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTATCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcgcg                                                                                                                                                    <  1:255439/150‑1 (MQ=255)
 tataAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAACGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcgcgc                                                                                                                                                   <  2:657060/150‑1 (MQ=255)
   taAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCgcgc‑gc                                                                                                                                                >  2:412768/1‑150 (MQ=255)
    aaaaGTTTTGTGCATTTAAAACATTACGCTTTGCAAAGGATTTTCATGGAAAGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACGAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑Gcc                                                                                                                                               <  2:545177/150‑1 (MQ=255)
     aaaGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCt                                                                                                                                              >  2:1165179/1‑150 (MQ=255)
       aGTTTTGTGCATTTGAAACATTACGCGTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGc                                                                                                                                            <  2:1533943/150‑1 (MQ=255)
        gTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCg                                                                                                                                           >  1:1282399/1‑150 (MQ=255)
         ttttGTGCATTTGAAACAATACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCgg                                                                                                                                          <  2:1510468/150‑2 (MQ=255)
          tttGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGAc                                                                                                                                         <  1:703264/150‑1 (MQ=255)
          tttGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggt                                                                                                                                         <  2:956576/150‑3 (MQ=255)
           ttGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggtg                                                                                                                                        <  2:982275/150‑4 (MQ=255)
             gtgCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCATGCCGCGC‑GCCTGCGACGCa                                                                                                                                      >  1:349259/1‑150 (MQ=255)
             gtgCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCa                                                                                                                                      <  2:1282399/150‑1 (MQ=255)
             gtgCATTTGAAACATTACGCTTGGCAAAGGATTATCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCa                                                                                                                                      <  2:470506/150‑1 (MQ=255)
              tgaaTTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAg                                                                                                                                     <  2:741913/147‑1 (MQ=255)
              tgCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGTAg                                                                                                                                     >  1:1711512/1‑150 (MQ=255)
              tgCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAg                                                                                                                                     >  2:853817/1‑150 (MQ=255)
                cATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGAc                                                                                                                                   >  1:1847754/1‑150 (MQ=255)
                cATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGAc                                                                                                                                   <  1:855647/150‑1 (MQ=255)
                cATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGAc                                                                                                                                   >  1:1301329/1‑150 (MQ=255)
                cATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATATTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTGTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGAc                                                                                                                                   <  2:1581405/150‑1 (MQ=255)
                    tGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACt                                                                                                                               <  2:1847754/150‑1 (MQ=255)
                     gAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCTC‑GCCTGCGACGCAGACAACTg                                                                                                                              >  1:1917767/1‑150 (MQ=255)
                     gAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTg                                                                                                                              >  1:1527297/1‑150 (MQ=255)
                        aCATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGGCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGAtt                                                                                                                           <  2:1464444/150‑1 (MQ=255)
                         cATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCTACGCAGACAACTGATTc                                                                                                                          >  1:1233704/1‑150 (MQ=255)
                           ttACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAATACTGCTACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCTACGCAGACAACTGATTCGc                                                                                                                        >  1:133937/1‑150 (MQ=255)
                                  ttGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActtc                                                                                                                 >  1:1631140/1‑150 (MQ=255)
                                   tGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGACACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                                                                >  2:907059/1‑150 (MQ=255)
                                        aGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgc                                                                                                           <  2:498182/150‑1 (MQ=255)
                                         ggATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCGTTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcg                                                                                                          >  2:764880/1‑150 (MQ=255)
                                         ggATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcg                                                                                                          >  2:390023/1‑150 (MQ=255)
                                         ggATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcg                                                                                                          <  2:1271202/150‑1 (MQ=255)
                                             tttCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGAt                                                                                                      <  1:897681/150‑1 (MQ=255)
                                              ttCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATa                                                                                                     <  2:1538505/150‑1 (MQ=255)
                                                 aTGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACACCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTAGCCAACTTCTTGAGCGCGATaaaa                                                                                                  >  1:1769787/1‑150 (MQ=255)
                                                  tGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAc                                                                                                 <  1:1881256/150‑1 (MQ=255)
                                                     aaCGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAcccc                                                                                              >  1:501569/1‑150 (MQ=255)
                                                        gTGCGAGTAAAATGCCGTCATCTTAGTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                                                                           <  1:1014853/150‑1 (MQ=255)
                                                          gCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAg                                                                                         <  2:1112361/150‑1 (MQ=255)
                                                          gCGAGTAAAATGCCGTCATATTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAg                                                                                         <  1:1680848/150‑1 (MQ=255)
                                                           cGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACAGCAGGCCGCGC‑GCCTGCGACCCAGACAACTTATTCGCCAACATCATGAGCGAGATAACACCCCGATatc                                                                                        >  1:338750/1‑148 (MQ=255)
                                                                 aaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAtttt                                                                                  <  2:1213569/150‑1 (MQ=255)
                                                                  aaTGCTGTCATCTTATTTGTATGACCAATAAGTGAGCATTGGATGAAAACTGATACTCCCTCTTTATAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTg                                                                                 <  1:963608/150‑1 (MQ=255)
                                                                    tGCCGTCTTCTTATTTTTATGAGAAATAAGTTATCATTGGTTTATATCTGTTACTCCCTCTTTAGATACACCGCAGGACGCTC‑GCCTGGGACGCAGACAACTGATTCGCGAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGtt                                                                               <  1:1536680/150‑1 (MQ=255)
                                                                    tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGtt                                                                               <  1:195661/150‑1 (MQ=255)
                                                                    tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGtt                                                                               <  1:102894/150‑1 (MQ=255)
                                                                     gCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGttt                                                                              <  2:205859/150‑1 (MQ=255)
                                                                     gCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACACCGTTAGACATTTTGttt                                                                              >  2:1363424/1‑150 (MQ=255)
                                                                        gTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTTCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATg                                                                           <  1:539169/150‑1 (MQ=255)
                                                                        gTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTCAGAAACACCGCAGGCCGCGC‑GCCTGAGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATg                                                                           <  2:1193485/150‑1 (MQ=255)
                                                                        gTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAAATGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAAATGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATg                                                                           <  2:1204379/150‑1 (MQ=255)
                                                                         tCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATgg                                                                          <  1:1162959/150‑1 (MQ=255)
                                                                          cATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggc                                                                         >  2:1042562/1‑150 (MQ=255)
                                                                          cATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggc                                                                         <  2:1426788/150‑1 (MQ=255)
                                                                           gtCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTTGCCATTTTGTTTATggcg                                                                        >  1:488989/2‑150 (MQ=255)
                                                                             cTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                                                                      <  1:1100/150‑1 (MQ=255)
                                                                             cTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                                                                      >  2:1685255/1‑150 (MQ=255)
                                                                                  ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTCGCCAGTTTGTTTATGGCGGCAgtcg                                                                 >  1:29219/1‑150 (MQ=255)
                                                                                   tGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgt                                                                <  1:1935013/150‑1 (MQ=255)
                                                                                      aTGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCgg                                                             <  1:1480055/150‑1 (MQ=255)
                                                                                        gACCAATAAGTGATCATTGGATGAAAGCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCa                                                           >  2:1590317/1‑150 (MQ=255)
                                                                                          ccacttaGTGATATTTGGATGAAAAGTGATACTCCCTCATTAGAAATACCGAAGGCCGCGG‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACg                                                         <  2:1471043/144‑1 (MQ=255)
                                                                                           cAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGc                                                        >  2:1195696/1‑150 (MQ=255)
                                                                                           cAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGc                                                        >  2:904117/1‑150 (MQ=255)
                                                                                               aaGTGATCATTGGATGAAAACTGATTCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgt                                                    >  2:1636793/1‑150 (MQ=255)
                                                                                               aaGTGATCAGTGGATGAAAACTGATTCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGCTTCGCCATCTTCTTGAGCGCGATAAACCCGCGTTAGCCTTTTTGTTTATGGCGGCTGTCGTCGGCACGGttgg                                                    >  2:411320/1‑149 (MQ=255)
                                                                                                 gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgttg                                                  >  2:675753/1‑150 (MQ=255)
                                                                                                      cATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTg                                             <  1:1110304/150‑1 (MQ=255)
                                                                                                       aTTGGATGAAAAATGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTgg                                            <  2:133937/150‑1 (MQ=255)
                                                                                                           gATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATATAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGc                                        >  1:1732797/1‑150 (MQ=255)
                                                                                                              gAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAGTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGt                                     >  2:1830600/1‑150 (MQ=255)
                                                                                                                 aaCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGc                                  <  2:691038/150‑1 (MQ=255)
                                                                                                                 aaCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGc                                  <  1:898776/150‑1 (MQ=255)
                                                                                                                    tGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                               <  2:1872684/150‑1 (MQ=255)
                                                                                                                     gATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCtttt                              <  1:1673096/150‑1 (MQ=255)
                                                                                                                        aCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGAc                           <  2:843469/150‑1 (MQ=255)
                                                                                                                            ccTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAg                       <  2:1887680/150‑1 (MQ=255)
                                                                                                                             ctctTTCGAAACACCGCAGGCCGCGC‑GCCTGCGACGCACACAACTGATTCGCCAACTTCTTGAGCGCGTTATAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGAGAAAgg                      <  2:949510/150‑1 (MQ=255)
                                                                                                                             ctctTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCGATTTTGTTTATGGCGGCAGTCGTCGGCGCGCTTGTTCGGCAGGCAGCGGTTGCTTTTGACTAAgg                      >  2:798551/1‑150 (MQ=255)
                                                                                                                                tttAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGgtgt                   >  1:1610929/1‑150 (MQ=255)
                                                                                                                                  tAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCCCGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCg                 >  1:30470/1‑150 (MQ=255)
                                                                                                                                   aGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGc                <  1:466693/150‑1 (MQ=255)
                                                                                                                                    tctgcgtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCAGGCAGCGGTTGCTTTTGACAAAGGTGTCGcc               >  2:1739444/8‑150 (MQ=255)
                                                                                                                                    gAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGcc               >  1:1585979/1‑150 (MQ=255)
                                                                                                                                    gAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGcc               >  2:797681/1‑150 (MQ=255)
                                                                                                                                      tgcgtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTg             >  1:1365817/6‑150 (MQ=255)
                                                                                                                                      aaCACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTg             <  2:1492751/150‑1 (MQ=255)
                                                                                                                                       acacCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACGAAGGTGTCGCCTgg            <  2:1527297/150‑1 (MQ=255)
                                                                                                                                       acacCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTgg            <  2:501569/150‑1 (MQ=255)
                                                                                                                                        cgtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCGTTTGACAAAGGTGTCGCCTGGt           >  1:1668161/4‑150 (MQ=255)
                                                                                                                                         gtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGtt          >  2:578819/3‑150 (MQ=255)
                                                                                                                                          ccGCAGGCCGCGC‑GCCTGCGACGCAGACGACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTg         >  1:373976/1‑150 (MQ=255)
                                                                                                                                           cGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGc        <  2:540924/150‑1 (MQ=255)
                                                                                                                                            gCAGGCCGCGC‑GCCTGCGACGCAGACAACTTATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCa       <  2:1913083/150‑1 (MQ=255)
                                                                                                                                               ggCCGCGA‑GCCTGCGACGCCGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTCTGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGaa    <  2:1997679/150‑1 (MQ=255)
                                                                                                                                                gCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCACCTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCCCGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAAc   >  1:1240923/1‑150 (MQ=255)
                                                                                                                                                 ccGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACTAAGGTGTCGCCTGGTTGCAGCAcc  >  2:362085/1‑150 (MQ=255)
                                                                                                                                                 ccGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAAcc  >  2:330973/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GTATAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACC  >  NC_000913/174999‑175293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑