Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 175,144 Δ1 bp 5.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. 5.9% 287.9 / 12.1 85coding (38/1422 nt)clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (33/47);  new base . (2/3);  total (35/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAA  >  NC_000913/175002‑175291
                                                                                                                                              |                                                                                                                                                    
taAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCcgcgc                                                                                                                                                 >  2:821891/1‑148 (MQ=255)
   aaGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGcc                                                                                                                                             >  2:311431/1‑148 (MQ=255)
   aaGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGcc                                                                                                                                             <  1:247355/148‑1 (MQ=255)
     gTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTg                                                                                                                                           <  1:577731/148‑1 (MQ=255)
        ttGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCgg                                                                                                                                        <  1:874312/148‑2 (MQ=255)
           tgCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggtgtt                                                                                                                                    <  4:542651/149‑6 (MQ=255)
             cATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAg                                                                                                                                   <  2:375563/148‑1 (MQ=255)
                       aTTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAtt                                                                                                                         <  1:900378/148‑1 (MQ=255)
                       aTTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAtt                                                                                                                         <  1:887383/148‑1 (MQ=255)
                        ttACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTGGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCg                                                                                                                       >  4:837112/1‑149 (MQ=255)
                          aCGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGcc                                                                                                                     <  3:487918/149‑1 (MQ=255)
                          aCGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGcc                                                                                                                     <  3:341521/149‑1 (MQ=255)
                            gCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCaa                                                                                                                   >  4:942741/1‑149 (MQ=255)
                            gCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCa                                                                                                                    <  2:482959/148‑1 (MQ=255)
                             cTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACa                                                                                                                                <  2:386387/135‑1 (MQ=255)
                             cTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACa                                                                                                                                >  1:386387/1‑135 (MQ=255)
                                tGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtc                                                                                                               <  3:574487/149‑1 (MQ=255)
                                 gCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCTGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtc                                                                                                               <  1:151813/148‑1 (MQ=255)
                                 gCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                                                              >  3:60031/1‑149 (MQ=255)
                                  cAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                                                              <  2:709158/148‑1 (MQ=255)
                                  cAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                                                              <  1:748073/148‑1 (MQ=255)
                                   aaaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtctt                                                                                                             >  2:503477/1‑148 (MQ=255)
                                    aaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTc                                                                                                                        >  3:666263/1‑136 (MQ=255)
                                    aaGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTc                                                                                                                        <  4:666263/136‑1 (MQ=255)
                                      ggATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgc                                                                                                         >  4:674339/1‑149 (MQ=255)
                                         ttttCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGcgcg                                                                                                      >  3:781103/1‑149 (MQ=255)
                                               tGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAGCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATaaaa                                                                                                <  4:870752/149‑1 (MQ=255)
                                                ggAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAc                                                                                               <  3:674339/149‑1 (MQ=255)
                                                       gCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                                                                         <  1:53651/148‑1 (MQ=255)
                                                       gCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                                                                        <  4:941012/149‑1 (MQ=255)
                                                           gTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGcc                                                                                     <  2:290910/148‑1 (MQ=255)
                                                              aaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTTATACTCCCTCTTTAGAAGCACCGCCGGCTGCGCGCCTCCTACGCAGACAACTGATTCGCCAACTTCTTGAGCGGGATAAAACCCCGGTAGCCAtg                                                                                  <  2:496876/148‑2 (MQ=255)
                                                                aTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTg                                                                               >  3:791528/1‑149 (MQ=255)
                                                                 tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAtttt                                                                                <  3:613535/147‑1 (MQ=255)
                                                                 tGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAtttt                                                                                >  4:613535/1‑147 (MQ=255)
                                                                  gCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGt                                                                              <  2:571894/148‑1 (MQ=255)
                                                                     gTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTa                                                                           <  1:800344/148‑1 (MQ=255)
                                                                           ttATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                                                                    >  3:30219/1‑149 (MQ=255)
                                                                            tATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                                                                    >  2:406976/1‑148 (MQ=255)
                                                                            tATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGAGTCGCCAACTTCATGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcggc                                                                    >  1:94553/1‑148 (MQ=255)
                                                                             aTTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAg                                                                  <  4:781103/149‑1 (MQ=255)
                                                                               ttGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgt                                                                 >  2:720641/1‑148 (MQ=255)
                                                                                  tATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgtc                                                             <  3:647141/149‑1 (MQ=255)
                                                                                  tATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgt                                                              <  2:626778/148‑1 (MQ=255)
                                                                                     gACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGGCAACTGATTCGGCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCgg                                                           >  2:690040/1‑148 (MQ=255)
                                                                                      aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                                                                        <  4:891145/118‑1 (MQ=255)
                                                                                      aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                                                                        >  3:891145/1‑118 (MQ=255)
                                                                                      aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc                                                          >  2:903138/1‑148 (MQ=255)
                                                                                      aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc                                                          >  2:286292/1‑148 (MQ=255)
                                                                                          aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGa                                                                                                                           >  4:421921/1‑79 (MQ=255)
                                                                                          aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGa                                                                                                                           <  3:421921/79‑1 (MQ=255)
                                                                                          aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCt                                                     >  3:215739/1‑149 (MQ=255)
                                                                                           tAAGTTATCATTGGATGAAGACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCtt                                                    <  4:215739/149‑1 (MQ=255)
                                                                                              gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgtt                                                 >  3:945847/1‑149 (MQ=255)
                                                                                              gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgtt                                                 <  4:395713/149‑1 (MQ=255)
                                                                                              gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgt                                                  <  1:286292/148‑1 (MQ=255)
                                                                                              gTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgt                                                  <  1:515509/148‑1 (MQ=255)
                                                                                               tGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgtt                                                 >  1:584921/1‑148 (MQ=255)
                                                                                                gATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgttg                                                <  1:821891/148‑1 (MQ=255)
                                                                                                  tCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGc                                             <  3:91490/149‑1 (MQ=255)
                                                                                                    aTTGGATGAAAACTGATACTCCATCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTg                                           <  3:541407/149‑1 (MQ=255)
                                                                                                      tGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTgg                                          >  1:942991/1‑148 (MQ=255)
                                                                                                      tGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTgg                                          >  1:720870/1‑148 (MQ=255)
                                                                                                          tGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGc                                      >  1:322775/1‑148 (MQ=255)
                                                                                                          tGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCg                                     >  4:950889/1‑149 (MQ=255)
                                                                                                              aaCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTg                                 >  4:530280/1‑149 (MQ=255)
                                                                                                               aCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTg                                 <  1:608120/148‑1 (MQ=255)
                                                                                                                   aTACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                             <  2:322775/148‑1 (MQ=255)
                                                                                                                   aTACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                             <  2:942991/148‑1 (MQ=255)
                                                                                                                     aCTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTg                           <  1:2458/148‑1 (MQ=255)
                                                                                                                     aCTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTg                           <  1:503243/148‑1 (MQ=255)
                                                                                                                          ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAg                     >  4:317719/1‑149 (MQ=255)
                                                                                                                          ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAg                     <  4:30219/149‑1 (MQ=255)
                                                                                                                          ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAg                     <  4:132342/149‑1 (MQ=255)
                                                                                                                           tctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAATTTCTTGAGCGCGATAAAACCCCGTTa                                                                                        <  4:578957/81‑1 (MQ=255)
                                                                                                                           tctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAATTTCTTGAGCGCGATAAAACCCCGTTa                                                                                        >  3:578957/1‑81 (MQ=255)
                                                                                                                                aGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCg               >  4:319881/1‑149 (MQ=255)
                                                                                                                                  ctgcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                             >  3:354514/7‑133 (MQ=255)
                                                                                                                                  ctgcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCttt                             <  4:354514/127‑1 (MQ=255)
                                                                                                                                  aaaCACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGc              <  2:401128/148‑1 (MQ=255)
                                                                                                                                    gcgtcGCAGGCGCGCGGCCTGCGGCGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCt            >  1:968382/5‑148 (MQ=255)
                                                                                                                                    acacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTg           >  3:190570/1‑149 (MQ=255)
                                                                                                                                     cacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTAACAAAGGTGTCGCCTgg          <  4:694870/149‑1 (MQ=255)
                                                                                                                                         gCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAGCTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGCCGCCTGGTTGc      <  3:837112/149‑1 (MQ=255)
                                                                                                                                             gCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGaa  <  4:60031/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TAAAAGTTTTGTGCATTTGAAACATTACGCTTTGCAAAGGATTTTCATGGAACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAA  >  NC_000913/175002‑175291

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