Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,470 C→A 10.5% L386I (CTA→ATA)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4700CA10.5% 323.8 / 22.7 95L386I (CTA→ATA) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base A (6/4);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAA  >  NC_000913/1666322‑1666617
                                                                                                                                                    |                                                                                                                                                   
tGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                                                                     <  4:273798/149‑1 (MQ=255)
  gTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                                   >  4:928232/1‑149 (MQ=255)
  gTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                                   >  4:135304/1‑149 (MQ=255)
   tCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                                                                  <  4:1147178/149‑1 (MQ=255)
    cGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACt                                                                                                                                                 <  3:927336/149‑1 (MQ=255)
        gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                                                                             <  4:769654/149‑1 (MQ=255)
        gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                                                                             <  4:1010877/149‑1 (MQ=255)
           tGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                                                                          >  1:89600/1‑149 (MQ=255)
                    ggTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                                                                                 <  1:969960/149‑1 (MQ=255)
                     gTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                                                                                >  3:762712/1‑149 (MQ=255)
                        cccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTTGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                                                                             <  3:894014/149‑1 (MQ=255)
                        cccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                                                                             <  3:863713/149‑1 (MQ=255)
                         ccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                                                                            >  1:697124/1‑149 (MQ=255)
                               ggCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                                                                      >  1:144297/1‑149 (MQ=255)
                                gCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                                                                     >  1:605880/1‑149 (MQ=255)
                                  gaagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGc                                                                                                                   >  1:10713/1‑149 (MQ=255)
                                   aagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCg                                                                                                                  >  1:448992/1‑149 (MQ=255)
                                    agaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGt                                                                                                                 <  4:293990/149‑1 (MQ=255)
                                     gaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                                                                                <  3:928232/149‑1 (MQ=255)
                                        aTTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                                                                             >  4:948131/1‑149 (MQ=255)
                                         ttACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                                                                            >  1:247520/1‑149 (MQ=255)
                                             acTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAttt                                                                                                        >  1:73818/1‑149 (MQ=255)
                                               ttttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                                                                               >  3:1097715/1‑108 (MQ=255)
                                               ttttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                                                                               <  4:1097715/108‑1 (MQ=255)
                                                tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                                                                     <  2:925575/149‑1 (MQ=255)
                                                tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                                                                     <  2:10713/149‑1 (MQ=255)
                                                tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGAGACGGTAATTTcgc                                                                                                     <  3:84459/149‑1 (MQ=255)
                                                 ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                                                                    <  2:73818/149‑1 (MQ=255)
                                                     tCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                                                                                >  2:169896/1‑149 (MQ=255)
                                                     tCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                                                                                >  1:463440/1‑149 (MQ=255)
                                                         aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                                                                            >  2:468149/1‑149 (MQ=255)
                                                            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAcc                                                                                         >  1:474734/1‑149 (MQ=255)
                                                             aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                                                                        >  3:274727/1‑149 (MQ=255)
                                                             aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                                                                        >  4:239459/1‑149 (MQ=255)
                                                             aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                                                                        >  3:534520/1‑149 (MQ=255)
                                                               cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                                                                      <  1:468149/149‑1 (MQ=255)
                                                                gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                                                                     <  4:534520/149‑1 (MQ=255)
                                                                gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                                                                     <  2:697124/149‑1 (MQ=255)
                                                                 ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                                                                    >  2:577715/1‑149 (MQ=255)
                                                                      gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                               <  2:247520/149‑1 (MQ=255)
                                                                        gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                                                                             >  1:953793/1‑149 (MQ=255)
                                                                           acTGTTGGCGGCGACCACGCACGCGCCGATTATGTCTACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                                                                     >  1:877428/1‑106 (MQ=255)
                                                                           acTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                                                                     <  2:877428/106‑1 (MQ=255)
                                                                              gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                       <  4:80984/149‑1 (MQ=255)
                                                                                       gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg                                                              >  2:518818/1‑149 (MQ=255)
                                                                                           acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgca                                                          <  3:239459/149‑1 (MQ=255)
                                                                                                 gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                                                                              <  2:672516/123‑1 (MQ=255)
                                                                                                 gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                                                                              >  1:672516/1‑123 (MQ=255)
                                                                                                   gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCt                                                  >  3:487420/1‑149 (MQ=255)
                                                                                                   gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCt                                                  <  2:953793/149‑1 (MQ=255)
                                                                                                      gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTa                                                                                             >  2:761831/1‑103 (MQ=255)
                                                                                                      gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTa                                                                                             <  1:761831/103‑1 (MQ=255)
                                                                                                      gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca                                                       <  2:185474/141‑1 (MQ=255)
                                                                                                      gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca                                                       >  1:185474/1‑141 (MQ=255)
                                                                                                           tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                          <  4:487420/149‑1 (MQ=255)
                                                                                                           tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                          <  4:762712/149‑1 (MQ=255)
                                                                                                             tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGt                                        <  2:463440/149‑1 (MQ=255)
                                                                                                             tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGt                                        <  2:448992/149‑1 (MQ=255)
                                                                                                               gACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg                                                    <  4:582611/135‑1 (MQ=255)
                                                                                                               gACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg                                                    >  3:582611/1‑135 (MQ=255)
                                                                                                               gACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTAt                                      >  3:1087541/1‑149 (MQ=255)
                                                                                                                aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTAtt                                     <  1:169896/149‑1 (MQ=255)
                                                                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGa                                           >  1:636968/1‑141 (MQ=255)
                                                                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGa                                           <  2:852802/141‑1 (MQ=255)
                                                                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGa                                           <  2:636968/141‑1 (MQ=255)
                                                                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGa                                           >  1:852802/1‑141 (MQ=255)
                                                                                                                     gatgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcccc                                >  3:612467/1‑149 (MQ=255)
                                                                                                                      atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCa                               <  1:577715/149‑1 (MQ=255)
                                                                                                                            tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTCATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCa                          <  1:1138074/148‑1 (MQ=255)
                                                                                                                            tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCa                          >  2:1138074/1‑148 (MQ=255)
                                                                                                                              tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCAc                         >  2:99704/1‑147 (MQ=255)
                                                                                                                              tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCAc                         <  1:99704/147‑1 (MQ=255)
                                                                                                                              tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGc                       <  2:474734/149‑1 (MQ=255)
                                                                                                                               gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCt                      >  2:27218/1‑149 (MQ=255)
                                                                                                                                    aaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                          >  1:460910/2‑124 (MQ=255)
                                                                                                                                    aaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                          <  2:460910/123‑1 (MQ=255)
                                                                                                                                    aaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCg                 >  4:755636/2‑149 (MQ=255)
                                                                                                                                     aCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATa                                                     >  3:310785/1‑112 (MQ=255)
                                                                                                                                     aCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATa                                                     <  4:310785/112‑1 (MQ=255)
                                                                                                                                     aCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCgg                >  4:981260/1‑149 (MQ=255)
                                                                                                                                      ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgc               >  1:923498/1‑149 (MQ=255)
                                                                                                                                       cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcg              >  2:201381/1‑149 (MQ=255)
                                                                                                                                       cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGc                                              <  3:1051801/117‑1 (MQ=255)
                                                                                                                                       cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGc                                              >  4:1051801/1‑117 (MQ=255)
                                                                                                                                       cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcg              >  2:825864/1‑149 (MQ=255)
                                                                                                                                        ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgc             >  1:248764/1‑149 (MQ=255)
                                                                                                                                         gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgcg            >  1:560761/1‑149 (MQ=255)
                                                                                                                                         gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgcg            >  1:387736/1‑149 (MQ=255)
                                                                                                                                             gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGc        <  2:144297/149‑1 (MQ=255)
                                                                                                                                              tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCa                          >  3:956965/1‑130 (MQ=255)
                                                                                                                                              tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCa                          <  4:956965/130‑1 (MQ=255)
                                                                                                                                                tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGa     >  2:849324/1‑149 (MQ=255)
                                                                                                                                                 cAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGAt    <  2:89600/149‑1 (MQ=255)
                                                                                                                                                  aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATa   >  2:172496/1‑149 (MQ=255)
                                                                                                                                                   gCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATaa  <  3:948131/149‑1 (MQ=255)
                                                                                                                                                   gCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATaa  <  1:518818/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAA  >  NC_000913/1666322‑1666617

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: ‑