Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,098,166 T→A 11.9% I46F (ATT→TTT)  wzzB ← regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,098,1660TA11.9% 207.3 / 3.2 74I46F (ATT→TTT) wzzBregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base T (42/21);  new base A (3/7);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 3.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATA  >  NC_000913/2097932‑2098385
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                             
aaGGTTCGATTGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑c                                                                                                                                                                                                                           <  2:255150/238‑1 (MQ=255)
aaGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑c                                                                                                                                                                                                                           >  1:255150/1‑238 (MQ=255)
aaGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGGGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACG‑GCCAAAAAACA‑AAAAG‑cc                                                                                                                                                                                                                          >  2:457413/1‑239 (MQ=255)
  ggTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAg                                                                                                                                                                                                                        >  1:649695/1‑239 (MQ=255)
     tcaaTGGTAAGTCTTTCTCGTTCTTCTTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTA‑AAGCGACCAATTAACTTCTCCTGCAAATCCGAAACTTTCGGTGCAGGCTATCCTTGGATAACATTCATAGCATTGTTTTAGCCAGCCATTTGCCCCCCGGCGGCCTGCGTGAAAACTTCTAGTTG‑CGTACACTTTTCCTTCGCTACT‑GCCAAATCTCC‑AAAAG‑CCAcgga                                                                                                                                                                                                                     >  2:851078/4‑238 (MQ=255)
     tCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATACAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTAACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAA‑CCAGGGc                                                                                                                                                                                                                     >  1:393124/1‑239 (MQ=255)
             aaGTTTTTCTCGTTCTTCCTGATTCTCCAGCGCTTTTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGCCTGAGTGATAATTGCT‑TTTGACTTCATTTTTCCCTACGCTACT‑CGCAAAAATCA‑AAAAG‑c                                                                                                                                                                                                                           >  2:142571/1‑225 (MQ=255)
                   ttttcGTTCTTCCTTATTACCCAGCTTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGa                                                                                                                                                                                                      <  2:911404/236‑1 (MQ=255)
                   tttcggTTTTTCCTGGTTATCCCGGTTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCGGCAAAGCCGAAACTTTCGGTGCAGCCTGACCAGAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTGCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGa                                                                                                                                                                                                      <  2:807811/234‑1 (MQ=255)
                   ttCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATTCTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGa                                                                                                                                                                                                      <  2:191154/239‑1 (MQ=255)
                   ttCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTAATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CAAGGGC‑AATAGCCAAACTGa                                                                                                                                                                                                      >  1:247657/1‑239 (MQ=255)
                                    tATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTc                                                                                                                                                                                     >  1:676089/1‑239 (MQ=255)
                                            gTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTTTCATCTTgc                                                                                                                                                                              >  1:756683/1‑238 (MQ=255)
                                            gTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCATTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCAACTTgcc                                                                                                                                                                             >  1:472974/1‑239 (MQ=255)
                                                    gaaTGCAGAGAAGGCAGAACAG‑AAGGCACCAGTTAACGTCGCATGAAAGTCCAATGCTTGCGTGGAAGCCTAACCAAAGATAACATGCATGTCATTGTTAGACACAGAAATTGGCCCCACGGCGGGCTAAGTGAGATTTGCT‑GTGGACGTCCATTTCTCCTGCGTCACT‑ACCAAATCACC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATAATGATTGTCTTTTGcccacgccac                                                                                                                                                                      <  1:41433/237‑1 (MQ=255)
                                                     aatgcaagaaaGGCA‑AACTGAAAGCCAACAATAAACGTCCCCTGAAAATCGTATACTTTGGGTGCAGCCTGCCCATAGAAAACATTCATGGCATTGTTATGGCCAGCATTTTGCCCCACGTCGGGCTGAGTAATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCA‑AATAG‑CCAGGGC‑AATAGCCACAATGACGAAAATGATGATTGTCATCTTGCCACGCCACaa                                                                                                                                                                    <  1:471856/230‑1 (MQ=255)
                                                       tGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCcaca                                                                                                                                                                     >  1:1086989/1‑236 (MQ=255)
                                                       tGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                                                                                                                    >  1:959614/1‑237 (MQ=255)
                                                             gaAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATTACGGAAATGATAATTGTCATCTTGCCACCCCAAAAATGCAcaa                                                                                                                                                            >  2:8656/1‑237 (MQ=255)
                                                             gaAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTa                                                                                                                                                            >  1:1077938/1‑239 (MQ=255)
                                                             gaAGGCAGAACTG‑AAGCCACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGCGTGATAATTTCT‑GTTGACGTCCATTTCTCCTTCCCTACG‑GCCAAATATCC‑AATAG‑CCCGGTC‑AATAGCCACAATGACGGAAATGATGATTTCCATTTTccaacgccaaacccacccaa                                                                                                                                                            >  2:884425/1‑228 (MQ=255)
                                                                            atcGACAAATTAGCGTCTCCGGCAAATCCGATCCTTTCGGGTCAGCCTGACAATAGATAACATTCATGGATTGGTTAAACCCAGCATTTGCCCCCACTGCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCCCCTTCGCTACT‑GCCAAATATCC‑AAAAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaa                                                                                                                                              <  1:216140/237‑1 (MQ=255)
                                                                            aGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGAGGATTGCCTTCTTTCCCCGACCCAACACAACTAGtaaa‑caaacaaaa                                                                                                                                             >  2:747732/1‑238 (MQ=255)
                                                                                aCCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATACCATTCATGGCATTGTTCTAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGTGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACTCCACAACTGCACTAGTaaatcaatcaaatcaa                                                                                                                                          >  1:424351/1‑239 (MQ=255)
                                                                                aCCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGCCATCTTGCCACGCCCCAACTGCACTAGTaaatcaatcaaataaa                                                                                                                                          >  1:162296/1‑239 (MQ=255)
                                                                                 ccAATTAACGTTTCCTGAAAATCCGATTTTTGCGGTGCAGCCTACCCAAAGAAAACATTCATGGTTTTGTTAAAGCCAGAAATTTGCCCCAGGTCGGGCCGGGTTATAATTGGT‑GTGGACGTCCATTTCTCCTCCGCTACT‑GCAAAATATCC‑AATAG‑CAAGGGC‑AATAGCCACAATGACGGAAATTATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaatcaat                                                                                                                                         <  1:140815/239‑1 (MQ=255)
                                                                                                   aaaCCCGAGACTTGCGGGTCAGCCTGACCATAGAAAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCGGAGTTATAATGGCT‑GTTGACGTCCTTTTCTCCTTCGCTACT‑GCAAAATATCC‑AATAC‑CCAGGGC‑AATAGCCAAAATGACGAAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGtt                                                                                                                       <  1:624360/239‑1 (MQ=255)
                                                                                                            ttttggggTGCGGCCTGACCATAGATAACTTTCTGTGCTTTGTTAGAGCCAGCAATGTGCCCCCTGGCGGGTGTAGGAATAATTGCT‑GTTGCCGCCCATTTCTCCTTCGCTTCT‑GCAAAATATCC‑AATAG‑CCGGGGA‑AATAGCCACAATGACGGAAATGATGATTTCCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGa                                                                                                              <  2:89309/233‑1 (MQ=255)
                                                                                                                  ggTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCa                                                                                                                                                               >  1:799976/1‑184 (MQ=255)
                                                                                                                      cacccggcccaaaaaTAACTTTCAGGGCATTGTTAAAGCCAGCAATTTCCCCCACGGTGGGTTGAGGAATAATCTCT‑GTTGACGCCAATTTCCCCTTCGCTTCG‑GCCAAATTCCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACattatt                                                                                                     <  1:187287/224‑1 (MQ=255)
                                                                                                                          tttataagacaccgcttttgtGGTATTTTTAGACCCAGCATTTTCCCCCGGTTGGGGCAGAGTAATAAGTGCT‑GTTGACGTCCATTCCCCCTCCGCAACG‑CCAAAATACCC‑AAAAG‑CCGGGCC‑AAAAGCCACAATGACTGAAATGATGATTGTCATCTGCCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTAtttt                                                                                                <  2:344013/219‑1 (MQ=255)
                                                                                                                          ttgacatcagacacgtttCAGGGCATTTTTACGCGCAGAAATTTGCCCCACGCCGGGCTGAGTAATAATTGGT‑GTCAGAGTCCATTTCGTCTGCGCTCCT‑CTCAAATATCC‑AATAG‑CCAGGAC‑AATAGCCAAAGTGACGGAAATGATTATTTTCATCTGGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTAtttt                                                                                                <  1:1006833/223‑1 (MQ=255)
                                                                                                                               cATAGTTAACTTTCGTGGATTTGTTAAAGCCGACAATTGTCCCCACGTCGGACTGGTTGATATTTGCT‑GTTGACGCCCATTTCACCTTCGCTCCT‑GCCAAATATCC‑AATGG‑CCAGGCA‑AATAGCCAAAATGAAGAAAATGATGATTGTCATCTTGCCCCGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTAct                                                                                           <  1:949185/239‑1 (MQ=255)
                                                                                                                                 aaaaTAACATTCATGGCTTTTTAATAGCCAGAAATTCCCCCCGGGTCGGGCTGAGTGATAATTGTT‑GTTGACTTCCATTTCCCCGTGGCTATT‑GCCAAATATCC‑AATAG‑CCAGGGA‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTActct                                                                                         <  1:8656/236‑1 (MQ=255)
                                                                                                                                       aCATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTcttcat                                                                                   >  1:623120/1‑239 (MQ=255)
                                                                                                                                             aTGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCt                                                                             >  1:539697/1‑239 (MQ=255)
                                                                                                                                               ggCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCAt                                                                                                                                                                                   >  2:334544/1‑135 (MQ=255)
                                                                                                                                               ggCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCAt                                                                                                                                                                                   <  1:334544/135‑1 (MQ=255)
                                                                                                                                                   ttGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCtcat                                                                                      >  1:344013/1‑224 (MQ=255)
                                                                                                                                                    tGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTc                                                                                                                                   >  1:434441/1‑178 (MQ=255)
                                                                                                                                                           gCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAACTCAATCTGTTCCGGGCCATGGTTTTGCCCAGACCAATTATTATTTTCACCTCTCTTCATTACCCTAAAGCACGGAAAtt                                                               >  2:591408/1‑239 (MQ=255)
                                                                                                                                                           gCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTAACGGAAAtt                                                               >  1:532795/1‑239 (MQ=255)
                                                                                                                                                           gCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTAACGGAAAtt                                                               >  2:712922/1‑239 (MQ=255)
                                                                                                                                                              aGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTcc                                                            >  1:978601/1‑239 (MQ=255)
                                                                                                                                                                cAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCCCAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCCGGGTTTGTCCCAGAAACATTATAATTTTTTACTCTCATCATTACCCTAACTGCAGGAAATTTCaaa                                                          >  2:288008/1‑239 (MQ=255)
                                                                                                                                                                cAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCaa                                                          >  1:555257/1‑239 (MQ=255)
                                                                                                                                                                    tggcccccgggcGGGCTGGGGTATAATTGCT‑TTTGACGTCCATTTCTCCTTCCTTACG‑GCCAAAT‑TCC‑AAAAGCCCAGGGA‑AATAGCCAAAATGACGGAAATGATGATTTTCATTTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGt                                                        <  1:591408/226‑1 (MQ=255)
                                                                                                                                                                      gcccccggcggggctggtgtaTAATTGCT‑GTTGACGTCCATTTCTCCTGCGCTACG‑GCCAAATACCC‑AAAAC‑CCAGGGA‑AATAGCCAAAAGGACGAAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGG‑TTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCa                                                   <  1:716687/219‑1 (MQ=255)
                                                                                                                                                                       ccccccGTCGGGCTGATTTATATTTGCT‑GGTGACGTCCTTTCCTCCTTCGCTACG‑GCCAATTACCA‑AATAG‑CCAGGGA‑AATCGCCACAATGACGGAAATGAGTATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCa                                                   <  2:1077938/239‑1 (MQ=255)
                                                                                                                                                                           aCGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAAGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGCCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCCTTACCCTAACTGACGGAAATTTCCaaacccccaaagt                                               >  2:397838/1‑228 (MQ=255)
                                                                                                                                                                                 ggCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTTACGGAAATTTCAAAGTCCCCAAAGTTTATCg                                         >  2:28330/1‑239 (MQ=255)
                                                                                                                                                                                   ctggtgaatatTTGTT‑GTGCACGTCCATTTC‑CCTTTCCTACTGGCAAAATTCCC‑AATAG‑CCAGGGC‑AATAGCCAAAAGGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGa                                        <  1:457413/231‑1 (MQ=255)
                                                                                                                                                                                            atttGCT‑GTTGACGCCAATTTCCCTTTCGCTACT‑GCAAAATTGCA‑AATAG‑CCAGGGC‑AATACCCACAATGACGGAAATGTTTATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGt                              <  1:462514/237‑1 (MQ=255)
                                                                                                                                                                                            aaTTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGt                              >  1:318311/1‑239 (MQ=255)
                                                                                                                                                                                             aattcgttttGACTTCCATTCCCCCTTCGCGTCT‑GCAAAATACCC‑AATAG‑CCAGGGA‑AATAGCCAAAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGt                              <  1:57065/231‑1 (MQ=255)
                                                                                                                                                                                               tGCT‑GTTGACGTCCTTTTCCACTTCGCTACT‑GCAAATTTTCA‑AAAGG‑CCAGGCA‑AATGACCACAATGACGAAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGaa                           <  1:747732/239‑1 (MQ=255)
                                                                                                                                                                                                tttttttgcgcctctttttcccttccttACT‑GACAAATATCC‑AATAG‑CAGGGGA‑AATAACCAAAATGAAGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGaa                           <  2:1086989/219‑1 (MQ=255)
                                                                                                                                                                                                gCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaat                                                                                                                                                  >  1:959868/1‑119 (MQ=255)
                                                                                                                                                                                                gCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTCCTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTACCGAAAAACCTccggaaa                          >  2:547825/1‑233 (MQ=255)
                                                                                                                                                                                                gCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAAt                          >  1:72788/1‑239 (MQ=255)
                                                                                                                                                                                                      tgccGTCCTGTTCCCCTCCGCTAC‑‑GCCCAATATCCAAATAG‑CCGGGCA‑AATGGCCAAATTGACGGAAATGATATTTGTCATCTGGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCta                      <  1:884425/235‑1 (MQ=255)
                                                                                                                                                                                                       gACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACattatt                                                                                                     >  1:837479/1‑158 (MQ=255)
                                                                                                                                                                                                       gACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGCCATgg                                                                                                                         >  1:38673/1‑138 (MQ=255)
                                                                                                                                                                                                          gTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACatta                                                                                                       >  2:827628/1‑153 (MQ=255)
                                                                                                                                                                                                          gTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AAGAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACatta                                                                                                       <  1:827628/153‑1 (MQ=255)
                                                                                                                                                                                                             cATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATaa              >  1:179898/1‑239 (MQ=255)
                                                                                                                                                                                                               aatttcccctgccttaaggcAAAAATCCAAAAAG‑CCGGGGA‑AAAAGCCACAATGAGGAAAATGATGATTGCCATTGTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTGCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATaa              <  1:989088/220‑1 (MQ=255)
                                                                                                                                                                                                                  ctcCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGAt         >  1:575898/1‑239 (MQ=255)
                                                                                                                                                                                                                        cGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGa                 >  1:963293/1‑225 (MQ=255)
                                                                                                                                                                                                                         gCTACT‑ACAAAAT‑TTC‑AAAAA‑CCAGGGCAAATAGCCAAAATAAAGAAAATTGTGTTTGCCATCTGGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAaata  <  2:318311/239‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaTATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTc                                                             >  2:596606/1‑171 (MQ=255)
                                                                                                                                                                                                                                   aaaTATCC‑AATAG‑CCAGGGA‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTc                                                             <  1:596606/171‑1 (MQ=255)
                                                                                                                                                                                                                                      tatCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGtt                                              >  1:412883/1‑183 (MQ=255)
                                                                                                                                                                                                                                            aaTAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAggc                 >  1:813272/1‑206 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                             
AAGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTG‑AAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCT‑GTTGACGTCCATTTCTCCTTCGCTACT‑GCCAAATATCC‑AATAG‑CCAGGGC‑AATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTATTATTTTCTACTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATA  >  NC_000913/2097932‑2098385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑