Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,358 T→A 5.5% L17Q (CTG→CAG)  ybgE → putative inner membrane protein in cydABX‑ybgE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,3580TA 5.5% 569.4 / 17.3 163L17Q (CTG→CAG) ybgEputative inner membrane protein in cydABX‑ybgE operon
Reads supporting (aligned to +/- strand):  ref base T (69/85);  new base A (4/5);  total (73/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 165 total aligned reads displayed.

ATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGT  >  NC_000913/774210‑774505
                                                                                                                                                    |                                                                                                                                                    
aTGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCt                                                                                                                                                      >  3:21649/1‑149 (MQ=255)
 tGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTg                                                                                                                                                     >  1:842333/1‑149 (MQ=255)
    aTTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCgg                                                                                                                                                  <  1:319427/149‑1 (MQ=255)
      tCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGgc                                                                                                                                                >  4:484604/1‑149 (MQ=255)
          ggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCGCAGGGGGCGCtt                                                                                                                                            <  4:1430483/149‑1 (MQ=255)
          ggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCGCAGGGGGCGCtt                                                                                                                                            <  2:439521/149‑1 (MQ=255)
            aaCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCGCAGGGGGCGCttgt                                                                                                                                          <  3:870281/149‑3 (MQ=255)
                cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTc                                                                                                                                     <  1:530839/149‑1 (MQ=255)
                cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTc                                                                                                                                     <  3:1067589/149‑1 (MQ=255)
                 ttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGTAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCTGGCGCTT‑TCCTGCg                                                                                                                                    >  4:967773/1‑149 (MQ=255)
                    ttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGa                                                                                                                                 >  4:1261127/1‑149 (MQ=255)
                    ttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGa                                                                                                                                 >  1:759167/1‑149 (MQ=255)
                        cTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGgc                                                                                                                             >  3:1265842/1‑149 (MQ=255)
                         tGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGgcg                                                                                                                            >  3:568777/1‑149 (MQ=255)
                           ttCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGTCTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCt                                                                                                                          >  1:1032994/1‑149 (MQ=255)
                           ttCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCt                                                                                                                          >  2:1562975/1‑149 (MQ=255)
                              gTTTGGAGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCt                                                                                                                       <  2:940568/149‑1 (MQ=255)
                               tttGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTg                                                                                                                      >  1:937896/1‑149 (MQ=255)
                                    ggTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                 <  2:389674/149‑1 (MQ=255)
                                    ggTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                 <  2:1032994/149‑1 (MQ=255)
                                          cACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGatg                                                                                                           >  3:263241/1‑149 (MQ=255)
                                             cgcgCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGatgtat                                                                                                        <  1:1568433/149‑1 (MQ=255)
                                             cgcgCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGatgtat                                                                                                        <  3:776701/149‑1 (MQ=255)
                                                 cTGGCGCTTGACCACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATTAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGCCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCGGGATGTATGttt                                                                                                    <  4:1567616/149‑1 (MQ=255)
                                                 cTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATGATCGCGACTTTGTATGCGGTAATTGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGttt                                                                                                    <  2:404087/149‑1 (MQ=255)
                                                      gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGa                                                                                               >  3:1483789/1‑149 (MQ=255)
                                                        ttGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAcc                                                                                             >  3:552972/1‑149 (MQ=255)
                                                                gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGt                                                                                     >  3:824746/1‑149 (MQ=255)
                                                                  cGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGttt                                                                                   <  1:1334640/149‑1 (MQ=255)
                                                                   gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGtttt                                                                                  <  3:187261/149‑1 (MQ=255)
                                                                     aTCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtctc                                                                                       <  4:1141258/142‑1 (MQ=255)
                                                                     aTCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGc                                                                                <  4:709601/149‑1 (MQ=255)
                                                                        aGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGt                                                                                                                     <  3:506214/109‑1 (MQ=255)
                                                                           cAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                          <  4:47771/149‑1 (MQ=255)
                                                                              aGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGAc                                                                       <  2:657055/149‑1 (MQ=255)
                                                                                ccGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCa                                                                     <  1:737035/149‑1 (MQ=255)
                                                                                  ggTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGt                                                                   <  2:490020/149‑1 (MQ=255)
                                                                                     cAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGaa                                                                <  4:816951/149‑1 (MQ=255)
                                                                                       agaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                              <  1:543571/149‑1 (MQ=255)
                                                                                        gaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTg                                                             >  2:1595851/1‑149 (MQ=255)
                                                                                        gaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTg                                                             >  1:250389/1‑149 (MQ=255)
                                                                                        gaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTg                                                             >  4:652536/1‑149 (MQ=255)
                                                                                         aagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                            >  3:1291134/1‑149 (MQ=255)
                                                                                          agaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGa                                                           >  1:1052950/1‑149 (MQ=255)
                                                                                           gaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGaa                                                          >  1:389579/1‑149 (MQ=255)
                                                                                            aCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGaaa                                                         <  4:1554951/149‑1 (MQ=255)
                                                                                              aTCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGAc                                                                       >  4:1322148/1‑133 (MQ=255)
                                                                                              aTCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGAc                                                                       <  3:1322148/133‑1 (MQ=255)
                                                                                              aTCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGa                                                           >  1:1526925/1‑145 (MQ=255)
                                                                                                   atgaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGt                                                                                                                                   <  2:710151/68‑1 (MQ=255)
                                                                                                      aGTAAGATTATGGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTgg                                                                                                 >  3:1375331/1‑99 (MQ=255)
                                                                                                      aGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTgg                                                                                                 <  4:1375331/99‑1 (MQ=255)
                                                                                                        tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                             >  1:1528493/1‑149 (MQ=255)
                                                                                                        tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                             <  4:620533/149‑1 (MQ=255)
                                                                                                        tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                             >  3:1191779/1‑149 (MQ=255)
                                                                                                          aGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                            >  2:1246080/1‑100 (MQ=255)
                                                                                                          aGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                            <  1:1246080/100‑1 (MQ=255)
                                                                                                          aGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTa                                           <  2:700105/149‑1 (MQ=255)
                                                                                                          aGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTa                                           <  4:1483789/149‑1 (MQ=255)
                                                                                                           gATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAt                                          >  3:565843/1‑149 (MQ=255)
                                                                                                            aTTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAtt                                         <  3:1036817/149‑1 (MQ=255)
                                                                                                            aTTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAtt                                         >  4:1036817/1‑149 (MQ=255)
                                                                                                             ttATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTgg                                                    <  1:1209115/137‑1 (MQ=255)
                                                                                                             ttATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTgg                                                    >  2:1209115/1‑137 (MQ=255)
                                                                                                              tATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGc                                       <  3:899642/149‑1 (MQ=255)
                                                                                                               aTCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCt                                      <  2:76227/149‑1 (MQ=255)
                                                                                                               aTCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCt                                      <  1:1478042/149‑1 (MQ=255)
                                                                                                                tCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTg                                     >  2:143295/1‑149 (MQ=255)
                                                                                                                 cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGa                                    <  4:552972/149‑1 (MQ=255)
                                                                                                                 cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGa                                    >  3:1125973/1‑149 (MQ=255)
                                                                                                                 cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGa                                    <  2:749869/149‑1 (MQ=255)
                                                                                                                   cgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtg                                  >  2:980472/1‑149 (MQ=255)
                                                                                                                    gACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtgt                                 >  2:182913/1‑149 (MQ=255)
                                                                                                                      cTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                               >  3:506948/1‑149 (MQ=255)
                                                                                                                       tttGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGttt                                            >  1:342369/1‑135 (MQ=255)
                                                                                                                       tttGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGttt                                            <  2:342369/135‑1 (MQ=255)
                                                                                                                       tttGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTggg                              <  2:1294764/149‑1 (MQ=255)
                                                                                                                        ttGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGc                             <  3:1437208/149‑1 (MQ=255)
                                                                                                                           tATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGATCTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGt                          >  4:124801/1‑149 (MQ=255)
                                                                                                                           tATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGt                          >  3:912536/1‑149 (MQ=255)
                                                                                                                             tGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCt                        <  2:759167/149‑1 (MQ=255)
                                                                                                                             tGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCt                        <  2:1492128/149‑1 (MQ=255)
                                                                                                                                     tGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgt                <  3:1454489/149‑1 (MQ=255)
                                                                                                                                     tGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgt                <  1:1595851/149‑1 (MQ=255)
                                                                                                                                      gggcAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtg               <  4:844078/146‑1 (MQ=255)
                                                                                                                                       ggaaaGCGCCCGCAGGGGGCGCTTGTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtg               >  3:434230/4‑149 (MQ=255)
                                                                                                                                       ggaaaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt              <  1:760792/146‑1 (MQ=255)
                                                                                                                                       gACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt              <  4:246590/149‑1 (MQ=255)
                                                                                                                                       gACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt              <  2:842333/149‑1 (MQ=255)
                                                                                                                                        gaaaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTg                  >  4:113822/3‑144 (MQ=255)
                                                                                                                                        gaaaGCGCCCGCAGGGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTg                  <  3:113822/142‑1 (MQ=255)
                                                                                                                                         cAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGa            <  3:652536/149‑1 (MQ=255)
                                                                                                                                          aaGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAt           <  1:661552/149‑1 (MQ=255)
                                                                                                                                           aGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAtt          <  1:182913/149‑1 (MQ=255)
                                                                                                                                            gcgcCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCt                        <  1:961883/134‑1 (MQ=255)
                                                                                                                                                ccccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                              <  1:1468097/92‑1 (MQ=255)
                                                                                                                                                ccccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                              >  2:1468097/1‑92 (MQ=255)
                                                                                                                                                 cccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAAc                                                               <  2:739558/90‑1 (MQ=255)
                                                                                                                                                 cccTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAAc                                                               >  1:739558/1‑90 (MQ=255)
                                                                                                                                                   cTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGt  <  1:1419564/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTT‑TCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGT  >  NC_000913/774210‑774505

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