Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,133,453 T→G 5.1% intergenic (‑56/+37) wcaA ← / ← wzc putative glycosyl transferase/colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,133,4530TG 5.1% 512.7 / 14.2 146intergenic (‑56/+37)wcaA/wzcputative glycosyl transferase/colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (68/70);  new base G (5/3);  total (73/73)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

TCCTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGC  >  NC_000913/2133306‑2133601
                                                                                                                                                   |                                                                                                                                                    
tCCTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                     >  2:577509/1‑149 (MQ=255)
  cTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                   <  2:951972/149‑1 (MQ=255)
   tGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                  >  1:345382/1‑149 (MQ=255)
   tGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCGCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGGGTGAACGCCTTATc                                                                                                                                                  >  2:1033747/1‑149 (MQ=255)
    ggCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 >  2:1128268/1‑149 (MQ=255)
       gcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                              >  2:1420586/1‑149 (MQ=255)
       gcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                              <  2:298185/149‑1 (MQ=255)
       gcAAAACCGATGTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                              <  1:654211/149‑1 (MQ=255)
         aaaaCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaa                                                                                                                                            <  1:90232/149‑3 (MQ=255)
            aCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                         >  3:1427363/1‑149 (MQ=255)
               gATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgggg                                                                                                                                      >  2:166250/1‑149 (MQ=255)
               gATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgggg                                                                                                                                      >  3:921922/1‑149 (MQ=255)
                         cccGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAAt                                                                                                                            <  4:257483/149‑1 (MQ=255)
                          ccGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATg                                                                                                                           >  1:170044/1‑149 (MQ=255)
                            gAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCa                                                                                                                         <  3:542508/149‑1 (MQ=255)
                             aaTCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAg                                                                                                                        <  3:1184311/149‑1 (MQ=255)
                             aaTCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAg                                                                                                                        <  3:17807/149‑1 (MQ=255)
                               tCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGt                                                                                                                                  <  3:387617/137‑1 (MQ=255)
                               tCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGt                                                                                                                                  >  4:387617/1‑137 (MQ=255)
                               tCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGc                                                                                                                      <  4:882917/149‑1 (MQ=255)
                                       tgttgGCGGTTGCAGGTCGGCATATAGATTGAGATCAGCGGAGTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTGGCTTTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTTCGAATGCAGGCGTCGTTAt                                                                                                              <  4:1463509/149‑1 (MQ=39)
                                         ttgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTAttt                                                                                                            <  1:1511479/149‑1 (MQ=255)
                                         ttgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTAttt                                                                                                            >  3:1271258/1‑149 (MQ=255)
                                         atgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTAGTTTGCTTGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTAttt                                                                                                            >  4:1462565/2‑149 (MQ=255)
                                           ggCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCg                                                                                                          >  1:1092704/1‑149 (MQ=255)
                                                ttCCAGGTCGACATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATcc                                                                                                     >  2:1259456/1‑149 (MQ=255)
                                                  ccAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGa                                                                                                   <  4:1276231/149‑1 (MQ=255)
                                                  ccAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGa                                                                                                   <  1:737049/149‑1 (MQ=255)
                                                   cAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGAc                                                                                                  <  3:579007/149‑1 (MQ=255)
                                                   cAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGAATCCGAc                                                                                                  <  2:919500/149‑1 (MQ=255)
                                                    agttCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACt                                                                                                 >  2:799970/4‑149 (MQ=255)
                                                    aGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACt                                                                                                 >  3:732687/1‑149 (MQ=255)
                                                    aGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACt                                                                                                 >  2:412955/1‑149 (MQ=255)
                                                    aGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGCATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACt                                                                                                 >  4:1215214/1‑149 (MQ=255)
                                                      gTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTta                                                                                               <  3:804202/149‑1 (MQ=255)
                                                           cATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattc                                                                                          <  1:438408/149‑1 (MQ=255)
                                                               taGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcgtat                                                                                      >  4:892722/1‑149 (MQ=255)
                                                                  aTTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcgtattcg                                                                                   >  4:450733/1‑149 (MQ=255)
                                                                   ttGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATattcgtattcgt                                                                                  >  2:19649/1‑149 (MQ=255)
                                                                        aTCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTAttt                                                                                                            <  3:25928/118‑1 (MQ=255)
                                                                        aTCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTAttt                                                                                                            >  4:25928/1‑118 (MQ=255)
                                                                        aTCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAg                                                                             <  3:1107734/149‑1 (MQ=255)
                                                                          cAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGcc                                                                           <  2:345382/149‑1 (MQ=255)
                                                                          cAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGcc                                                                           <  1:392016/149‑1 (MQ=255)
                                                                           aGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCa                                                                          >  1:953165/1‑149 (MQ=255)
                                                                              ggATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATaa                                                                       >  4:302467/1‑149 (MQ=255)
                                                                               gATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcgtatt                                                                                     >  2:970443/1‑134 (MQ=255)
                                                                                aTTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATc                                                                     <  4:921922/149‑1 (MQ=255)
                                                                                aTTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATc                                                                     >  2:1297078/1‑149 (MQ=255)
                                                                                 ttGTTTTTCATAGGTAGTCCTCCCGGCTTATGACCTGGTTTGCTCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATcc                                                                    <  4:236993/149‑1 (MQ=255)
                                                                                    tttttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAg                                                                                                                        <  3:377263/94‑1 (MQ=255)
                                                                                    tttttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAg                                                                                                                        >  4:377263/1‑94 (MQ=255)
                                                                                    tttttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATa                                                                              >  1:306796/1‑136 (MQ=255)
                                                                                    tttttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTAATTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGa                                                                 >  1:568310/1‑149 (MQ=255)
                                                                                     ttttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGAt                                                                <  2:1436917/149‑1 (MQ=255)
                                                                                      tttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATa                                                               <  2:1366434/149‑1 (MQ=255)
                                                                                        tCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATAcg                                                             >  3:181679/1‑149 (MQ=255)
                                                                                          attggTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATAcgcg                                                           >  3:948354/4‑149 (MQ=255)
                                                                                          aTAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATAcgcg                                                           >  2:336023/1‑149 (MQ=255)
                                                                                           tAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACgcgc                                                          <  1:166250/149‑1 (MQ=255)
                                                                                            aGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCt                                                         <  2:1092704/149‑1 (MQ=255)
                                                                                               ttgttCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGgc                                                      >  2:1706339/1‑149 (MQ=255)
                                                                                                     tcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcgtattcg                                                                                   <  1:1358516/114‑1 (MQ=255)
                                                                                                     tcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtatt                                                                                           >  4:468094/1‑106 (MQ=255)
                                                                                                     tcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtatt                                                                                           <  3:468094/106‑1 (MQ=255)
                                                                                                     tcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcggcg                                                >  2:589035/1‑149 (MQ=255)
                                                                                                     tcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcggcg                                                >  1:540280/1‑149 (MQ=255)
                                                                                                     tcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcggcg                                                >  4:213607/1‑149 (MQ=255)
                                                                                                      cccGGCTTACGAGCTGGTTGGCGTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCCAGcgcgtatc                                               <  1:917928/149‑5 (MQ=39)
                                                                                                       ccGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcgg                                                  <  1:83053/145‑1 (MQ=255)
                                                                                                       ccGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGa                                              >  4:84861/1‑149 (MQ=255)
                                                                                                       ccGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGa                                              >  4:1201746/1‑149 (MQ=255)
                                                                                                        cGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGaa                                             <  1:1128268/149‑1 (MQ=255)
                                                                                                        cGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGaa                                             <  4:1271258/149‑1 (MQ=255)
                                                                                                           cTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGAt                                          <  1:767287/149‑1 (MQ=255)
                                                                                                            ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATa                                                                              >  4:262614/1‑112 (MQ=255)
                                                                                                            ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATa                                                                              <  3:262614/112‑1 (MQ=255)
                                                                                                            ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATg                                         <  4:1638508/149‑1 (MQ=255)
                                                                                                               tGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAg                                      <  1:1436855/149‑1 (MQ=255)
                                                                                                                 aGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGcgcg                                                   >  1:1305393/1‑134 (MQ=255)
                                                                                                                 aGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGcgcg                                                   <  2:1305393/134‑1 (MQ=255)
                                                                                                                   cTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCa                                  <  3:451111/149‑1 (MQ=255)
                                                                                                                     ggTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGa                                >  3:150972/1‑149 (MQ=255)
                                                                                                                      gTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGaa                               >  3:1478005/1‑149 (MQ=255)
                                                                                                                      gTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGaa                               <  3:1201746/149‑1 (MQ=255)
                                                                                                                      gTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGaa                               <  1:1264349/149‑1 (MQ=255)
                                                                                                                         tgcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCa                            <  1:1262634/149‑1 (MQ=255)
                                                                                                                          gcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCa                            <  2:248006/148‑1 (MQ=255)
                                                                                                                            ttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCAccc                         <  2:540280/149‑1 (MQ=255)
                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCtt                      >  4:1681075/1‑149 (MQ=255)
                                                                                                                                    taaggcgTTCACGCCGCATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGt                                     >  2:815638/4‑129 (MQ=255)
                                                                                                                                    taaggcgTTCACGCCGCATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGt                                     <  1:815638/126‑1 (MQ=255)
                                                                                                                                        gcgTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGaa             <  4:1514619/149‑1 (MQ=255)
                                                                                                                                           tGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATAc          >  1:1096140/1‑149 (MQ=255)
                                                                                                                                               cGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGtt      <  4:726089/149‑1 (MQ=255)
                                                                                                                                               cGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGtt      <  1:39635/149‑1 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGttt     >  2:773001/1‑149 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTaa                                                                                <  4:287326/72‑1 (MQ=255)
                                                                                                                                                   ttATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGc  <  1:966937/149‑1 (MQ=255)
                                                                                                                                                   ttATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGc  >  1:1118946/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGC  >  NC_000913/2133306‑2133601

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: