Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,470 C→A 6.0% 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,4700CA 6.0% 647.8 / 23.2 183L386I (CTA→ATA) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (91/81);  new base A (7/4);  total (98/85)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-02

Only 100 of 188 total aligned reads displayed.

TAATGATCATTGCCGGGATCTTCCTTTGTAAACTGTGTGCCGTGCTGGCGAGTAGCGGTTCCGGCGCACCCGGTGGGGTCTTTACACCGACGCTATTTATCGGTCTTGCCATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCACCAATGAGCGGTGCCGGAAGTTTTTTACCCAGCACATCAATGATTTCCGCATAATGCGCCAGTCCTGGGTTGATTCGGTTAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTT  >  NC_000913/1666198‑1666739
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                             
tAATGATCATTGCCGGGATCTTCCTTTGTAAACTGTGTGCCGTGCTGGCGAGTAGCGGTTCCGGCGCACCCGGTGGGGTCTTTACACCGACGCTATTTATCGGTCTTGCCATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                                                                                                                                                             <  2:74295/275‑1 (MQ=255)
       cATTGCCGGGATCTTCCTTTGTAAACTGTGTGCCGTGCTGGCGAGTAGCGGTTCCGGCGCACCCGGTGGGGTCTTTACACCGACGCTATTTATCGGTCTTGCCATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTAATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTTTGTCGACGTTGAGGAAATGTGAAATGACCGGGACGTATCGGCGAACCCCCg                                                                                                                                                                                                                                                                      >  4:155782/1‑275 (MQ=255)
                                                                    cccGGTGGGGTCTTTACACCGACGCTATTTATCGGTCTTGCCATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGGGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTTATTGCCTGCGTAATTGCGTCGGTAATTTCGCGCACGTTACCCCGTCACccta                                                                                                                                                                                                          >  4:20071/1‑274 (MQ=255)
                                                                      cGGTGGGGTCTTTACACCGATGCTATTTATCGGTCTTGCCATTGGCAGGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCGACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                                                                                                                                                         <  2:86553/273‑1 (MQ=255)
                                                                              tCTTTACACCGACGCTAGTTATCGGTCTTGCCATTGGCAGGTGTTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                                                                                                                                                                                                <  4:137440/274‑1 (MQ=255)
                                                                                        gACGCTATTTATCGGTCTTGCCATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGCTTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgagc                                                                                                                                                                                     >  1:29200/1‑275 (MQ=255)
                                                                                                            ccATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTa                                                                                                                                                                 <  2:231515/275‑1 (MQ=255)
                                                                                                                 ggCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTACTTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCCCCCTTCCCAGCATAGCTAAGCCGGTGTATTGcc                                                                                                                                                            >  4:260991/1‑275 (MQ=255)
                                                                                                                            tGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                                                                                                                                                                                               >  7:1328861/1‑149 (MQ=255)
                                                                                                                              gTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                                                                                                                                                             >  8:134625/1‑149 (MQ=255)
                                                                                                                              gTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                                                                                                                                                             >  6:1144903/1‑149 (MQ=255)
                                                                                                                                  tAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTTCCCCAGTTCACGCTGTTCg                                                                                                                                           >  2:56428/1‑275 (MQ=255)
                                                                                                                                    gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGCTTATGTCGACGTTGATGAGATGTGACATGACCGGGCAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTTCGTCGGTAATTTCGCGGACATTACACCGTGACTCTATCTCCCGCCGCCCCACTCCGCAGCCTAGCTACGCGGCGGTGTTTCCACAGTTACCGCTGTTCGgc                                                                                                                                         >  4:332954/1‑275 (MQ=255)
                                                                                                                                    gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                                                                                                                                                                                                       <  7:646416/149‑1 (MQ=255)
                                                                                                                                    gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                                                                                                                                                                                                       <  7:134625/149‑1 (MQ=255)
                                                                                                                                      cTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                                                                                                                                                                                                     >  8:575666/1‑149 (MQ=255)
                                                                                                                                       tGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTCTCAGc                                                                                                                                                                                                                                                                               >  7:208633/1‑138 (MQ=255)
                                                                                                                                         ggATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCTCTAGCACACTGCGCAGCATAGCTCAGCGGCTGCATTGCTCCAGTTCACGCTGTTCGGcgtgcg                                                                                                                                    >  1:16225/1‑275 (MQ=255)
                                                                                                                                         ggATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGGGACTCTATCTACCGCCAGCACCCTGCGCAGCATCGCTAAGCGGATGTATTGCCCCAGTTACAGCTGTTCGgcggcgg                                                                                                                                    >  4:55943/1‑271 (MQ=255)
                                                                                                                                           aTTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTat                                                                                                                                                                                                                                                                <  5:832493/149‑1 (MQ=255)
                                                                                                                                            ttATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                                                                                                                                                                                               <  7:144020/149‑1 (MQ=255)
                                                                                                                                              aTGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattga                                                                                                                                                                                                                                                             <  5:656855/149‑1 (MQ=255)
                                                                                                                                                ggTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                                                                                                                                                                                                           >  6:612172/1‑149 (MQ=255)
                                                                                                                                                  ttCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                                                                                                                                                                                                         >  8:646826/1‑149 (MQ=255)
                                                                                                                                                  ttCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                                                                                                                                                                                                         <  8:1414675/149‑1 (MQ=255)
                                                                                                                                                   tCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGcc                                                                                                                                                                                                                                                        >  6:1162062/1‑149 (MQ=255)
                                                                                                                                                    cccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                                                                                                                                                                                                       <  8:920277/149‑1 (MQ=255)
                                                                                                                                                    cccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                                                                                                                                                         >  2:116099/1‑195 (MQ=255)
                                                                                                                                                      cTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                                                                                                                                                                                                      >  6:1052503/1‑132 (MQ=255)
                                                                                                                                                      cTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                                                                                                                                                                                                      <  5:1052503/132‑1 (MQ=255)
                                                                                                                                                       tGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                                                                                                                                                                                                    <  5:1237093/149‑1 (MQ=255)
                                                                                                                                                       tGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                                                                                                                                                                                                    >  7:945786/1‑149 (MQ=255)
                                                                                                                                                                agaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGt                                                                                                                                                                                                                                           <  6:65345/149‑1 (MQ=255)
                                                                                                                                                                    aTTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                                                                                                                                                                                                       >  8:526981/1‑149 (MQ=255)
                                                                                                                                                                      tACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCg                                                                                                                                                                                                                                                                      >  8:22729/1‑116 (MQ=255)
                                                                                                                                                                           ttttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                                                                                                                                                                                                >  5:1454600/1‑149 (MQ=255)
                                                                                                                                                                            tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                                                                                                                                                                                               <  7:111516/149‑1 (MQ=255)
                                                                                                                                                                               aCTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                                                                                                                                                                                                            <  5:556979/149‑1 (MQ=255)
                                                                                                                                                                                 tCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                                                                                                                                                                                                          >  6:1297169/1‑149 (MQ=255)
                                                                                                                                                                                    gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTa                                                                                                                                                                                                                       >  7:1243762/1‑149 (MQ=255)
                                                                                                                                                                                    gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTa                                                                                                                                                                                                                       >  5:1178192/1‑149 (MQ=255)
                                                                                                                                                                                         aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAGGTATTGACCCAGTTCCCGCTGTTCGGCGCGCGGCAGGTAAGGCAGTTCACCATTGAGCGGTTCCGTAAGTTTTTTGCCCa                                                                                    >  4:217095/1‑275 (MQ=255)
                                                                                                                                                                                          ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                                                                                                                                                                                                 <  5:1293001/149‑1 (MQ=255)
                                                                                                                                                                                          ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                                                                                                                                                                                                 <  7:442944/149‑1 (MQ=255)
                                                                                                                                                                                          ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCACCAATGAGCGGTGCCGGAAGTTTTTTACCCAg                                                                                   <  1:73609/275‑1 (MQ=255)
                                                                                                                                                                                           ctcattGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                                                                                                                                                                                                >  6:444527/6‑149 (MQ=255)
                                                                                                                                                                                            gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                                                                                                                                                                                               >  6:428545/1‑149 (MQ=255)
                                                                                                                                                                                                 ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                                                                                                                                                                                                          <  8:1293343/149‑1 (MQ=255)
                                                                                                                                                                                                  gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                                                                                                                                                         >  5:945731/1‑149 (MQ=255)
                                                                                                                                                                                                   cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                                                                                                                                                                                                        <  6:1330047/149‑1 (MQ=255)
                                                                                                                                                                                                    gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                                                                                                                                                                                                       >  6:1255368/1‑149 (MQ=255)
                                                                                                                                                                                                        cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCGCCAATGAGCGGTGCCGGAAGTTTTTTACCCAGCACATCAATGAtt                                                                      <  2:155782/274‑1 (MQ=255)
                                                                                                                                                                                                         tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                                                                                                                                                                                                  >  5:153072/1‑149 (MQ=255)
                                                                                                                                                                                                          gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                                                                                                                                                 <  5:78396/149‑1 (MQ=255)
                                                                                                                                                                                                          gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                                                                                                                                                 <  5:1396133/149‑1 (MQ=255)
                                                                                                                                                                                                           ttGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                                                                                                                                                                                                <  5:374806/149‑1 (MQ=255)
                                                                                                                                                                                                           ttGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                                                                                                                                                                                                >  6:926660/1‑149 (MQ=255)
                                                                                                                                                                                                               cggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTaa                                                                                                                                                                                                                                     >  8:730185/1‑108 (MQ=255)
                                                                                                                                                                                                               cggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTaa                                                                                                                                                                                                                                     <  7:730185/108‑1 (MQ=255)
                                                                                                                                                                                                                     ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg                                                                                                                                                                                      >  8:299027/1‑149 (MQ=255)
                                                                                                                                                                                                                       acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgca                                                                                                                                                                                    >  5:301390/1‑149 (MQ=255)
                                                                                                                                                                                                                         gcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagc                                                                                                                                                                                  >  6:591312/1‑149 (MQ=255)
                                                                                                                                                                                                                            cgcgCCGATTATGTGGACGTTTATGATATTTGAAAGGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCACCAATGAGCGGTGCCGGAAGTTTTTTACCCAGCACATCAATGATTTCCGCATAATGCGCCAGTcc                                                  <  4:127749/274‑1 (MQ=255)
                                                                                                                                                                                                                            cgcgCCGATTATGTCGACGTTGATGATATGTTAAATGACCGGGAAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCTTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATa                                                                                                                                                                               <  6:1233043/149‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg                                                                                                                                                                              <  8:391172/149‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg                                                                                                                                                                              <  5:352098/149‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg                                                                                                                                                                              <  7:1176292/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTaa                                                                                                                                                                          <  5:549312/149‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGc                                                                                                                                                                        <  6:318923/149‑1 (MQ=255)
                                                                                                                                                                                                                                       tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                                                                                                                                                    <  8:458412/149‑1 (MQ=255)
                                                                                                                                                                                                                                       tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                                                                                                                                                    <  7:526981/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGt                                                                                                                                                                  >  8:1323710/1‑149 (MQ=255)
                                                                                                                                                                                                                                          cGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTa                                                                                                                                                                 <  5:859030/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ttGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCATGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcc                                                                                                                                                            >  8:591930/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAg                                                                                                                                                        <  6:637133/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGc                                                                                                                                                 <  6:229416/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCt                                                                                                                                                >  6:889890/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTg                                                                                                                                               >  6:1189140/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcg                                                                                                                                        >  2:358788/1‑156 (MQ=255)
                                                                                                                                                                                                                                                              aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGtt                                                                                                                                             <  8:658568/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGtt                                                                                                                                             <  7:1290841/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGtt                                                                                                                                             >  7:31178/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                                                                                                                                                    <  5:208129/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                gACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                                                                                                                                                                                                <  2:132021/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                gACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCg                                                                                                                                           <  5:250332/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                gACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGcgcgcg                                                                                                                                    <  4:358788/156‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCgg                                                                                                                                          >  6:95001/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ccGGGGATTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgc                                                                                                                                         >  6:792257/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGGCCCAGTTCACGCTGTTCGgcgc                                                                                                                                       >  5:633721/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                                                                                                                                                         <  7:735398/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg                                                                                                                                                                                        <  5:1449893/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg                                                                                                                                                                                        >  6:1449893/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                      ggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGc                                                                                                                                                             <  6:681025/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGc                                                                                                                                  >  7:1112502/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCa                                                                                                                                 >  5:1316789/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGa                                                                                                                               >  8:441482/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGa                                                                                                                               >  7:975624/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCACCAATGAGCGGTGCCGGAAGTTTTTTACCCAGCACATCAATGATTTCCGCATAATGCGCCCGTCCTGGGTTGATTCGGTTAGCCACCCAGCCAATGAGCGGACACCCGTCGtt  >  4:108888/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                               gCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATaa                                                                                                                            >  8:868265/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                cTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCACCAATGAGCGGTGCCGGAAGTTTTTTACCCAGCACATCaa                                                                           >  4:214566/1‑197 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                             
TAATGATCATTGCCGGGATCTTCCTTTGTAAACTGTGTGCCGTGCTGGCGAGTAGCGGTTCCGGCGCACCCGGTGGGGTCTTTACACCGACGCTATTTATCGGTCTTGCCATTGGCATGTTGTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGGCAGTTCACCAATGAGCGGTGCCGGAAGTTTTTTACCCAGCACATCAATGATTTCCGCATAATGCGCCAGTCCTGGGTTGATTCGGTTAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTT  >  NC_000913/1666198‑1666739

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