Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,002,332 A→G 15.6% R425R (CGT→CGC fliC ← flagellar filament structural protein (flagellin)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,002,3320AG15.6% 2016.0 / 313.6 790R425R (CGT→CGCfliCflagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (363/303);  new base G (65/59);  total (428/362)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

Only 100 of 824 total aligned reads displayed.

GGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTT  >  NC_000913/2002184‑2002480
                                                                                                                                                    |                                                                                                                                                    
     atCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGaa                                                                                                                                                 <  2:2712302/149‑1 (MQ=255)
      tCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAAt                                                                                                                                                >  3:561430/1‑149 (MQ=255)
        tGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAAttt                                                                                                                                              <  1:502153/149‑1 (MQ=255)
           gcTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTc                                                                                                                                           >  4:1530380/1‑149 (MQ=255)
           gcTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTc                                                                                                                                           <  4:3136668/149‑1 (MQ=255)
            cTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGa                                                                                                                                                  >  2:4950558/1‑141 (MQ=255)
             tttCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCt                                                                                                                                          <  1:3196321/148‑1 (MQ=255)
             tttCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTa                                                                                                                                         >  2:7282400/1‑149 (MQ=255)
             tttCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTa                                                                                                                                         >  1:5146063/1‑149 (MQ=255)
                 gACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGaa                                                                                                                                                 <  4:351642/137‑1 (MQ=255)
                  aCATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGAt                                                                                                                                    >  2:7033137/1‑149 (MQ=255)
                    atatTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGc                                                                                                                                  <  3:1336266/149‑1 (MQ=255)
                     tatTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGAt                                                                                                                                    >  4:1264164/1‑146 (MQ=255)
                      atTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGa                                                                                                                                <  4:2037476/149‑1 (MQ=255)
                      atTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGa                                                                                                                                <  1:6717993/149‑1 (MQ=255)
                      atTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGa                                                                                                                                >  3:2037476/1‑149 (MQ=255)
                       ttGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGGCAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGAt                                                                                                                               >  3:41304/1‑149 (MQ=255)
                       ttGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACGGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGAt                                                                                                                               >  4:2912813/1‑149 (MQ=255)
                       ttGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGAt                                                                                                                               >  2:32506/1‑149 (MQ=255)
                        tGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATa                                                                                                                              <  1:3472451/149‑1 (MQ=255)
                          gACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTAc                                                                                                                                        <  4:6959505/137‑1 (MQ=255)
                          gACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGc                                                                                                                            <  1:6470332/149‑1 (MQ=255)
                             acTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGaa                                                                                                                                                 >  4:4444174/1‑125 (MQ=255)
                             acTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATc                                                                                                                         >  4:2565100/1‑149 (MQ=255)
                             acTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATc                                                                                                                         >  1:1288428/1‑149 (MQ=255)
                             acTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATc                                                                                                                         <  4:5666762/149‑1 (MQ=255)
                               ttCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTAc                                                                                                                                        >  1:1036598/1‑132 (MQ=255)
                               ttCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGt                                                                                                                       >  1:242854/1‑149 (MQ=255)
                               ttCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGt                                                                                                                       >  4:5471026/1‑149 (MQ=255)
                                  ggTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCa                                                                                                                    <  2:2892280/149‑1 (MQ=255)
                                  ggTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCa                                                                                                                    <  3:4350427/149‑1 (MQ=255)
                                            tCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAg                                                                                                          >  3:6091219/1‑149 (MQ=255)
                                             cGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGc                                                                                                         <  3:1611951/149‑1 (MQ=255)
                                               gCGTCCTGAATACTGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCgg                                                                                                       >  4:3584138/1‑149 (MQ=255)
                                                cGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGTGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGa                                                                                                      <  4:5222346/149‑1 (MQ=255)
                                                    cTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCg                                                                                                  >  3:5172378/1‑149 (MQ=255)
                                                      gAATACGGGACTGCTCTTCATACAGGGTGTTAGTGGTGTTGTTCAGTTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTg                                                                                                >  2:6719837/1‑149 (MQ=255)
                                                       aaTACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTgg                                                                                               >  3:4820055/1‑149 (MQ=255)
                                                       aaTACGGGACTGCGCCTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCCTTCAGCGGATCCGTgg                                                                                               <  1:6027944/149‑1 (MQ=255)
                                                         tACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGtt                                                                                             >  1:1897992/1‑149 (MQ=255)
                                                         tACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGtt                                                                                             >  1:4351566/1‑149 (MQ=255)
                                                           cGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGtttt                                                                                           >  1:7133315/1‑149 (MQ=255)
                                                            gggACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTa                                                                                          <  2:4553982/149‑1 (MQ=255)
                                                            gggACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTa                                                                                          <  3:3388359/149‑1 (MQ=255)
                                                             ggACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTAc                                                                                         <  2:3965279/149‑1 (MQ=255)
                                                              gACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTAc                                                                                         >  3:6808171/1‑148 (MQ=255)
                                                                cTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCAt                                                                                      >  1:6275647/1‑149 (MQ=255)
                                                                 tGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCAtt                                                                                     <  2:6471731/149‑1 (MQ=255)
                                                                 tGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCAtt                                                                                     <  2:138221/149‑1 (MQ=255)
                                                                 tGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCAtt                                                                                     >  1:5804651/1‑149 (MQ=255)
                                                                    gcTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTc                                                                                                            <  4:1059110/123‑1 (MQ=255)
                                                                    gcTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGc                                                                                  >  4:3830448/1‑149 (MQ=255)
                                                                     cTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGa                                                                                                                                                  <  3:6412342/84‑1 (MQ=255)
                                                                     cTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGa                                                                                                                                                  >  4:6412342/1‑84 (MQ=255)
                                                                     cTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCa                                                                                 >  1:68266/1‑149 (MQ=255)
                                                                      ttCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCaa                                                                                <  4:2293208/149‑1 (MQ=255)
                                                                      ttCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCaa                                                                                >  2:3544035/1‑149 (MQ=255)
                                                                         aGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAAcag                                                                             <  4:1835689/149‑1 (MQ=255)
                                                                         aGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAAcag                                                                             >  2:5239592/1‑149 (MQ=255)
                                                                          gACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAAcagc                                                                            >  3:7020046/1‑149 (MQ=255)
                                                                          gACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGaa                                                                                                                                                 <  4:4936747/80‑1 (MQ=255)
                                                                          gACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAAcagc                                                                            >  2:1474566/1‑149 (MQ=255)
                                                                             aGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCACACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGtt                                                                                             >  2:5937774/1‑129 (MQ=255)
                                                                                ttGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCag                                                                      >  2:5548361/1‑149 (MQ=255)
                                                                                 tGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCaga                                                                     <  3:354151/149‑1 (MQ=255)
                                                                                   gTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCt                                                                                                                                          <  2:2270925/78‑1 (MQ=255)
                                                                                   gTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCaga                                                                     >  3:2178611/1‑147 (MQ=255)
                                                                                     aGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCt                                                                 >  3:2398818/1‑149 (MQ=255)
                                                                                       tggtgTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTca                                                               <  1:1526370/149‑1 (MQ=255)
                                                                                        ggtgTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGAAAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTcac                                                              <  2:4224312/149‑1 (MQ=255)
                                                                                         gtgtTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTcacc                                                             <  1:778655/149‑1 (MQ=255)
                                                                                           gttgttCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTcaccac                                                           <  3:2912813/149‑1 (MQ=255)
                                                                                           gttgttCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTcacca                                                            <  1:986835/148‑1 (MQ=255)
                                                                                            ttgttCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCATCGAGGGAAGAACGGAATTTGTCTAc                                                                                                                                        <  2:4766231/71‑1 (MQ=255)
                                                                                            ttgttCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCaccacc                                                          >  2:1694039/1‑149 (MQ=255)
                                                                                              gttCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACcag                                                        <  3:2607730/149‑1 (MQ=255)
                                                                                                  aGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGt                                                    >  3:1829835/1‑149 (MQ=255)
                                                                                                         tAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCt                                             <  1:2128435/149‑1 (MQ=255)
                                                                                                          aaCCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAAcagca                                                                           <  1:3827939/118‑1 (MQ=255)
                                                                                                           aCCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTc                                                                        <  3:1118550/120‑1 (MQ=255)
                                                                                                             cgcgGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGttt                                         >  2:909834/1‑149 (MQ=255)
                                                                                                              gcgGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTg                                        >  1:5372894/1‑149 (MQ=255)
                                                                                                              gcgGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTg                                        >  4:4681008/1‑149 (MQ=255)
                                                                                                                 gAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGc                                       >  2:510797/1‑147 (MQ=255)
                                                                                                                       aGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACATCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTAct                               <  2:6174864/149‑2 (MQ=255)
                                                                                                                       aGACGGTTTTGCACCGCACCGAGGGAAGAGCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTAcc                               >  3:5466759/1‑149 (MQ=255)
                                                                                                                       aGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCg                                                                                                        <  1:4748948/76‑1 (MQ=255)
                                                                                                                           ggTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCa                           <  4:41304/149‑1 (MQ=255)
                                                                                                                           ggTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCa                           >  1:2907614/1‑149 (MQ=255)
                                                                                                                                   accgcaccgAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATc                   >  4:3292360/1‑149 (MQ=255)
                                                                                                                                   accgcaccgAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATc                   >  2:2483974/1‑149 (MQ=255)
                                                                                                                                      gcaccgAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGc                >  3:6602786/1‑149 (MQ=255)
                                                                                                                                      gcaccgAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGc                <  1:2689622/149‑1 (MQ=255)
                                                                                                                                             gggAAGAACGGAATTTGTCTACAGAAGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCg         >  2:1969094/1‑149 (MQ=255)
                                                                                                                                              ggAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGt        >  4:6647803/1‑149 (MQ=255)
                                                                                                                                               gaagaaCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTa       <  3:3830448/149‑1 (MQ=255)
                                                                                                                                               gaagaaCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTa       <  1:4159639/149‑1 (MQ=255)
                                                                                                                                               gaagaaCGGAATTTGTCTACAGATGCGATAACATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTa       >  4:730786/1‑149 (MQ=255)
                                                                                                                                                aagaaCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTAt      >  4:6141412/1‑149 (MQ=255)
                                                                                                                                                    aCGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGttt  >  3:3153218/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAGCATCGTCCAGCGCTTTCAGCGGATCCGTGGTTTTACCATTTGCAACAGCAGTCAGAGCCTCACCACCAGCAGTCAAACCTGTTTGCAGATTACCGCCATTTAAATCGGCAGAATCGTATGTTT  >  NC_000913/2002184‑2002480

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