Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 175,145 Δ1 bp 16.7% coding (39/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.16.7% 618.9 / 129.4 210coding (38/1422 nt)clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (96/77);  new base . (17/18);  total (115/95)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 217 total aligned reads displayed.

TTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGC  >  NC_000913/175045‑175239
                                                                                                   |                                                                                                
ttCATGGAACGTGCGAGTAAAATGCCGGCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGc                                                                                                  >  2:2187540/1‑100 (MQ=255)
  cATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcg                                                                                                >  1:1058689/1‑100 (MQ=255)
    tGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcgcg                                                                                              <  1:1264753/100‑1 (MQ=255)
     ggAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcgcgc                                                                                             >  2:706801/1‑100 (MQ=255)
      gAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑cgcgc‑gg                                                                                          <  1:1530446/100‑2 (MQ=255)
        cgatctgAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑cgcgc‑ggcc                                                                                        <  2:1559050/94‑4 (MQ=255)
        aCGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑Gcc                                                                                         <  2:1211201/100‑1 (MQ=255)
        aCGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑cgcgc‑ggcc                                                                                        <  2:816384/100‑4 (MQ=255)
          gTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTg                                                                                       <  2:329234/100‑1 (MQ=255)
             cGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATCCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGa                                                                                    <  2:3230288/100‑1 (MQ=255)
             cGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGa                                                                                    <  2:3053273/100‑1 (MQ=255)
             cGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCgg                                                                                    <  1:2707166/100‑2 (MQ=255)
             cGAGTAAAATGCCGTCATCTTATATGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAACCAACGCACGCCGCTC‑Gcc                                                                                         >  2:3194584/1‑95 (MQ=255)
              gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGAc                                                                                   >  2:1084723/1‑100 (MQ=255)
              gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggt                                                                                   <  2:3950197/100‑3 (MQ=255)
              gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggt                                                                                   <  1:2028522/100‑3 (MQ=255)
              gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggt                                                                                   <  1:1598177/100‑3 (MQ=255)
              gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggt                                                                                   <  1:2867593/100‑3 (MQ=255)
               atctaaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACg                                                                                  <  1:1160937/96‑1 (MQ=255)
               aGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACg                                                                                  <  2:573671/100‑1 (MQ=255)
               aGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggtg                                                                                  <  2:3600112/100‑4 (MQ=255)
               aGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGG‑CGCGCGGCCTGCggtg                                                                                  <  1:3116217/100‑4 (MQ=255)
                gTAAAAAGCCGTCATCTTATTTGTATGACCCATAAGTGATCACTGGACGACAACTGATCCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGcgccgc                                                                                 <  1:2150702/100‑1 (MQ=255)
                   aaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACgagat                                                                              >  2:1160928/1‑96 (MQ=255)
                    aaTGCCGTCAGCTTAGTTGTATGACCAATAAGTGAACAGTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGACGCGC‑GCCTGCGACGCAGAc                                                                             >  2:428539/1‑100 (MQ=255)
                     attccttcatcaTCTTTGTATGACCAACAAGTGATCATTGGATGCAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACa                                                                            <  2:518791/94‑1 (MQ=255)
                        ccGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACt                                                                         <  2:1299508/100‑1 (MQ=255)
                           tCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGAt                                                                      >  1:343374/1‑100 (MQ=255)
                              tCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCg                                                                   >  2:107831/1‑100 (MQ=255)
                                ttATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGcc                                                                 >  1:460472/1‑100 (MQ=255)
                                ttATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGcc                                                                 <  1:2522517/100‑1 (MQ=255)
                                 tagttGTATGACCAATAGGTGTTAATTGGATGAAAACTGATACTCCCTCTTTCGAGACACCGCAGGCCGCGG‑GACTGCGACGCAGACAACTGATTCGCCa                                                                >  2:2719807/4‑100 (MQ=255)
                                 tATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCa                                                                <  2:3121701/100‑1 (MQ=255)
                                 tATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCa                                                                <  1:480935/100‑1 (MQ=255)
                                    ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAAct                                                             <  2:1726673/100‑1 (MQ=255)
                                     tGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActt                                                            >  1:813419/1‑100 (MQ=255)
                                     tGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActt                                                            >  1:1406441/1‑100 (MQ=255)
                                     tGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActt                                                            >  1:2691440/1‑100 (MQ=255)
                                      gTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActtc                                                           >  2:3822810/1‑100 (MQ=255)
                                       tATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAActtct                                                          >  1:3188609/1‑100 (MQ=255)
                                         tGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTg                                                        <  2:1219846/100‑1 (MQ=255)
                                         tGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTg                                                        >  1:3234979/1‑100 (MQ=255)
                                           aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTTGCCAActtctcg                                                        >  2:1757364/1‑96 (MQ=255)
                                           aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAg                                                      <  1:3350778/100‑1 (MQ=255)
                                             cAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcg                                                    >  1:2752585/1‑100 (MQ=255)
                                              aaTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcgc                                                   <  1:59983/100‑1 (MQ=255)
                                               aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGcgcg                                                  >  1:1662146/1‑100 (MQ=255)
                                                 aaGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGAt                                                <  2:2549383/100‑1 (MQ=255)
                                                  aGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATa                                               >  2:1549807/1‑100 (MQ=255)
                                                   gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATaa                                              >  1:3855685/1‑100 (MQ=255)
                                                   gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATaa                                              >  1:2301207/1‑100 (MQ=255)
                                                    tGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATaaa                                             >  1:2825544/1‑100 (MQ=255)
                                                      aTCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAc                                           <  1:975551/100‑1 (MQ=255)
                                                        cATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAccc                                         >  1:2740959/1‑100 (MQ=255)
                                                          ttGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCg                                       <  1:2358962/100‑1 (MQ=255)
                                                           tGGATGAAACCTGATAGTCCCTCTTTAGGAACACCGCAGGCCGCGC‑GCCTGCGACGCCGAGCACTAATTCGCCGACTTCTTGAGCGCGCTGACAcccact                                      >  2:1941862/1‑97 (MQ=255)
                                                            ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                     >  1:2971118/1‑100 (MQ=255)
                                                            ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                     >  1:412889/1‑100 (MQ=255)
                                                            ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                     >  2:1848936/1‑100 (MQ=255)
                                                             gATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                    >  1:1312776/1‑100 (MQ=255)
                                                             gATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                    <  2:1339242/100‑1 (MQ=255)
                                                               tGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGc                                  <  2:2931015/100‑1 (MQ=255)
                                                                gAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCCCGATAAAACCCCGTTAGcc                                 >  2:3592186/1‑100 (MQ=255)
                                                                      tGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTg                           <  1:3549476/100‑1 (MQ=255)
                                                                         tACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCATTAGCCATTTTGttt                        >  1:3887204/1‑100 (MQ=255)
                                                                          aCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTTGCCATTTTGTTTa                       <  2:2441831/100‑1 (MQ=255)
                                                                          aCTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTa                       <  2:118092/100‑1 (MQ=255)
                                                                              ccTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggc                   >  1:3295519/1‑100 (MQ=255)
                                                                               ttccgatctaaCACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                  <  1:1458466/91‑1 (MQ=255)
                                                                               ctctTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                  >  2:3708089/1‑100 (MQ=255)
                                                                               ctctTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                  >  1:3713785/1‑100 (MQ=255)
                                                                               ctctTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                  >  2:3068419/1‑100 (MQ=255)
                                                                                tctTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcgg                 <  2:630412/100‑1 (MQ=255)
                                                                                 ctTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                >  2:787300/1‑100 (MQ=255)
                                                                                 ctTGAGAAACACCCCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                <  2:3429404/100‑1 (MQ=255)
                                                                                    tgtctgcgtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgt             >  2:1259512/10‑100 (MQ=38)
                                                                                    tAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgt             >  1:559307/1‑100 (MQ=255)
                                                                                    cttctgcgtcGCAGG‑CGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgt             <  1:3675497/91‑1 (MQ=37)
                                                                                     aGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtc            <  2:2825526/100‑1 (MQ=255)
                                                                                     aGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtc            >  2:2744564/1‑100 (MQ=255)
                                                                                      gAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg           >  2:2956276/1‑100 (MQ=255)
                                                                                      gAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg           >  1:109311/1‑100 (MQ=255)
                                                                                      gAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg           <  2:3759940/100‑1 (MQ=255)
                                                                                      gAAACACCGCAGGACGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg           >  1:2655410/1‑100 (MQ=255)
                                                                                        aaCACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcgagatcggaa         >  2:1458457/1‑91 (MQ=255)
                                                                                        aaCACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgtc         >  1:1738746/1‑100 (MQ=255)
                                                                                         gcgtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGtcgtcg        >  1:3401840/5‑100 (MQ=255)
                                                                                         gcgtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGtcgtcg        >  2:2058735/5‑100 (MQ=255)
                                                                                          cacCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCgg       <  1:2290904/100‑1 (MQ=255)
                                                                                          cacCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCgg       <  1:3571821/100‑1 (MQ=255)
                                                                                           gtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc      >  2:3114284/3‑100 (MQ=255)
                                                                                           gtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc      >  1:1908787/3‑100 (MQ=255)
                                                                                           gtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc      >  2:1982260/3‑100 (MQ=255)
                                                                                           gtcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc      >  1:633268/3‑100 (MQ=255)
                                                                                           acCGCGGGCCGCGC‑GCCTGCGACGCAGACAACTGAGTCGTGAACTTCTTGAGCGCGATAAGACCACGTTAGCCATTTTGTTTATGGCGGCAGTAGTCGGc      <  2:2650394/100‑1 (MQ=255)
                                                                                           acCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc      <  1:3616048/100‑1 (MQ=255)
                                                                                           acCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc      >  2:235311/1‑100 (MQ=255)
                                                                                            tcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCa     >  1:940615/2‑100 (MQ=255)
                                                                                            tcGCAGG‑CGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCa     >  1:2039972/2‑100 (MQ=255)
                                                                                               cAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGc  >  1:707226/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGC‑GCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGC  >  NC_000913/175045‑175239

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