Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 563,387 G→A 13.0% A20T (GCA→ACA)  sfmF → FimA homolog, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913563,3870GA13.0% 2294.0 / 266.1 832A20T (GCA→ACA) sfmFFimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base G (353/369);  new base A (55/53);  total (409/423)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-03

Only 100 of 848 total aligned reads displayed.

TGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGC  >  NC_000913/563239‑563535
                                                                                                                                                    |                                                                                                                                                    
tGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATggg                                                                                                                                                      >  2:4025347/1‑149 (MQ=255)
 ggCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGc                                                                                                                                                     >  4:2267216/1‑149 (MQ=255)
 ggCCTATCAGTATTACCGGGCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTATTCTTTAGCTGTTTCTGAGGCGTACTGTGGAGTTCGGGTGGATggga                                                                                                                                                     >  2:1566074/1‑148 (MQ=255)
    cTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGt                                                                                                                                                  >  2:464583/1‑149 (MQ=255)
     tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGtt                                                                                                                                                 <  1:1527244/149‑1 (MQ=255)
     tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGtt                                                                                                                                                 >  2:2207947/1‑149 (MQ=255)
     tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGtt                                                                                                                                                 >  1:6803245/1‑149 (MQ=255)
     tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGACTACTGTGGAGTTCCAGTGGATGGGCAGtt                                                                                                                                                 <  3:4030489/149‑1 (MQ=255)
      aTCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTg                                                                                                                                                >  2:5414457/1‑149 (MQ=255)
       tCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGa                                                                                                                                               >  3:2120666/1‑149 (MQ=255)
           tATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCt                                                                                                                                           >  1:5947725/1‑149 (MQ=255)
              tACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa                                                                                                                                        <  3:829240/149‑1 (MQ=255)
                ccGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAg                                                                                                                                       >  4:1542421/1‑148 (MQ=255)
                    tCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACg                                                                                                                                  >  3:6010508/1‑149 (MQ=255)
                     cAAAAACCGACCGAAGGGCCGTTTAGCTCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGACAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGa                                                                                                                                 <  2:6585799/149‑1 (MQ=255)
                     cAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGa                                                                                                                                 <  2:4487703/149‑1 (MQ=255)
                      aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAt                                                                                                                                >  4:1051743/1‑149 (MQ=255)
                      aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAt                                                                                                                                <  3:3297414/149‑1 (MQ=255)
                      aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAt                                                                                                                                <  1:460433/149‑1 (MQ=255)
                       aaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAtt                                                                                                                               >  2:7013035/1‑149 (MQ=255)
                          accgaccgAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAat                                                                                                                            >  3:5531104/1‑149 (MQ=255)
                          accgaccgAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAat                                                                                                                            <  1:3111695/149‑1 (MQ=255)
                           ccgaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTa                                                                                                                                        >  3:2913905/1‑136 (MQ=255)
                            cgaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAatat                                                                                                                          <  2:5847901/149‑1 (MQ=255)
                            cgaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAatat                                                                                                                          <  4:5587545/149‑1 (MQ=255)
                            cgaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAatat                                                                                                                          <  3:7056349/149‑1 (MQ=255)
                            cgaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAatat                                                                                                                          <  3:3501024/149‑1 (MQ=255)
                                    gggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGc                                                                                                                  <  1:2348131/149‑1 (MQ=255)
                                     ggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCa                                                                                                                 >  1:5606710/1‑149 (MQ=255)
                                      gCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCAc                                                                                                                <  1:2001901/149‑1 (MQ=255)
                                        cGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAAc                                                                                                                                   >  2:1876614/1‑128 (MQ=255)
                                           ttAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTaa                                                                                                           <  2:2734647/149‑1 (MQ=255)
                                             aGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACg                                                                                                         >  3:2434424/1‑149 (MQ=255)
                                                  gcTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACgttgtt                                                                                                    <  4:4148253/149‑1 (MQ=255)
                                                  gcTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACgttgtt                                                                                                    <  4:466307/149‑1 (MQ=255)
                                                  gcTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACattgtt                                                                                                    <  1:605229/149‑1 (MQ=255)
                                                  gcTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGAGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACgttgtt                                                                                                    <  2:3828037/149‑1 (MQ=255)
                                                    tGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTCGa                                                                                                  >  3:4739202/1‑149 (MQ=255)
                                                          aTTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTtctc                                                                                            <  1:1046272/149‑1 (MQ=255)
                                                          aTTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTtctc                                                                                            <  2:1680795/149‑1 (MQ=255)
                                                           tttACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTcc                                                                                           <  2:4224738/149‑1 (MQ=255)
                                                              aCGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGt                                                                                        >  1:4467017/1‑149 (MQ=255)
                                                                  gTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCg                                                                                    >  2:13933/1‑149 (MQ=255)
                                                                  gTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCg                                                                                    >  3:13154/1‑149 (MQ=255)
                                                                  gTCGATTATCAATGAGGTAAGGAGAATGAGAAGAGTACTCTTGAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTGGAACCTTTAGGAACGATTAATATCAATTTGCACGGTAAAGTTGTTGATTTCTCCTGTACCg                                                                                    >  2:5273606/1‑149 (MQ=255)
                                                                       ttATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAac                                                                               >  1:5051449/1‑149 (MQ=255)
                                                                             aTGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACgttgtt                                                                                                    <  4:5370787/122‑1 (MQ=255)
                                                                             aTGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACgttgtt                                                                                                    >  3:5370787/1‑122 (MQ=255)
                                                                              tGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCgg                                                                        >  3:1781494/1‑149 (MQ=255)
                                                                                aGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAt                                                                                                                                >  1:5362092/1‑91 (MQ=255)
                                                                                  gTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGatat                                                                    >  4:3436560/1‑149 (MQ=255)
                                                                                     cGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAAAGTTGTTGACTTCTCCTGTACCGTAAACACAGCGGATATTGa                                                                 >  4:5809011/1‑149 (MQ=255)
                                                                                      ggAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGAt                                                                <  1:3351189/149‑1 (MQ=255)
                                                                                      ggAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGAt                                                                >  4:2724031/1‑149 (MQ=255)
                                                                                       gagaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATa                                                               >  2:970014/1‑149 (MQ=255)
                                                                                        agaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATaa                                                              <  1:6604249/149‑1 (MQ=255)
                                                                                         gaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAg                                                             >  3:1538782/1‑149 (MQ=255)
                                                                                         gaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAg                                                             >  2:3354721/1‑149 (MQ=255)
                                                                                               gaagaGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTa                                                       <  4:1742304/149‑1 (MQ=255)
                                                                                                    gTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTtctc                                                                                            <  4:3221272/107‑1 (MQ=255)
                                                                                                    gTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGAttt                                                  <  4:13154/149‑1 (MQ=255)
                                                                                                     tACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTa                                                 >  4:3021841/1‑149 (MQ=255)
                                                                                                     tACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTa                                                 <  1:1827061/149‑1 (MQ=255)
                                                                                                     tACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTa                                                 <  2:6329471/149‑1 (MQ=255)
                                                                                                      aCTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATATGACGGTAGATTTAg                                                <  3:512591/149‑1 (MQ=255)
                                                                                                      aCTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAg                                                <  4:765787/149‑1 (MQ=255)
                                                                                                      aCTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAg                                                <  1:6301923/149‑1 (MQ=255)
                                                                                                       ctctTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAgg                                               >  1:3547253/1‑149 (MQ=255)
                                                                                                       ctctTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAgg                                               <  1:6376954/149‑1 (MQ=255)
                                                                                                        tctTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGc                                              >  4:3746104/1‑149 (MQ=255)
                                                                                                        tctTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTAc                                                                                      <  2:3127633/109‑1 (MQ=255)
                                                                                                         ctTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCa                                             <  3:374275/149‑1 (MQ=255)
                                                                                                          tttAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCttt                                                                                                                                         <  4:4460501/56‑1 (MQ=255)
                                                                                                          tttAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGa                                                                       >  1:4157203/1‑122 (MQ=255)
                                                                                                          tttAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAg                                            >  3:1177267/1‑149 (MQ=255)
                                                                                                           ttAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGTATATTGATAAGACGGTAGATTTAGGCAGa                                           <  3:5790202/149‑1 (MQ=255)
                                                                                                           ttAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGa                                           >  3:4552695/1‑149 (MQ=255)
                                                                                                            tAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGAt                                          <  2:128708/149‑1 (MQ=255)
                                                                                                             aGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATg                                         <  4:6330489/149‑1 (MQ=255)
                                                                                                             aGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATg                                         >  3:6849687/1‑149 (MQ=255)
                                                                                                              gCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGCCAGTTGACCCTTTAGTAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCGCCTGTGCCGTAAACGCGGCGGATATTGATAAGACTGGAGATTTAGGCAGATgg                                        <  4:4869081/149‑1 (MQ=255)
                                                                                                               cTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAg                                            >  3:4271692/1‑144 (MQ=255)
                                                                                                                 gTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCt                                     <  4:5150077/149‑1 (MQ=255)
                                                                                                                 gTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCt                                     >  1:4845261/1‑149 (MQ=255)
                                                                                                                  tttCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGa                                                            <  1:271498/125‑1 (MQ=255)
                                                                                                                  tttCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTa                                    >  2:4873586/1‑149 (MQ=255)
                                                                                                                     cTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGa                                 <  2:2869682/149‑1 (MQ=255)
                                                                                                                                tgtgGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTg                      >  1:1351386/1‑149 (MQ=255)
                                                                                                                                tgtgGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTg                      >  3:2706472/1‑149 (MQ=255)
                                                                                                                                   ggAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAAc                   <  4:2063452/149‑1 (MQ=255)
                                                                                                                                    gAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACg                  >  4:670404/1‑149 (MQ=255)
                                                                                                                                     aGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAAcac                 >  4:401263/1‑147 (MQ=255)
                                                                                                                                       ttCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTCACgttgtt                                                                                                    <  3:5014432/64‑1 (MQ=255)
                                                                                                                                          cAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAg                                                <  3:6333396/113‑1 (MQ=255)
                                                                                                                                          cAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCg            >  1:379155/1‑149 (MQ=255)
                                                                                                                                            gTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGTCTACGACACAACTACTGAACGCTGGCGAt          >  2:6025092/1‑149 (MQ=255)
                                                                                                                                            gTGGATGGACAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGAt          <  1:2909521/149‑1 (MQ=255)
                                                                                                                                               gATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATAcc       >  2:3517684/1‑149 (MQ=255)
                                                                                                                                                   ggCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTtctc                                                                                            <  1:764272/60‑1 (MQ=255)
                                                                                                                                                    gCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGc  >  4:611489/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGC  >  NC_000913/563239‑563535

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