Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,353 A→T 6.1% L174Q (CTG→CAG)  gatD ← galactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,3530AT 6.1% 360.6 / 10.2 98L174Q (CTG→CAG) gatDgalactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (48/44);  new base T (3/3);  total (51/47)
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.99e-01

TGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAT  >  NC_000913/2172212‑2172494
                                                                                                                                             |                                                                                                                                             
tGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                        <  4:378740/149‑1 (MQ=255)
tGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                        >  3:1210615/1‑149 (MQ=255)
              gcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTGTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                                                                          <  4:1210615/149‑1 (MQ=255)
               cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaat                                                                                                                         >  1:882815/1‑149 (MQ=255)
                  aTTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataat                                                                                                                      >  2:1692170/1‑149 (MQ=255)
                       gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAAtaataataa                                                                                                                 >  4:607563/1‑149 (MQ=255)
                            tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                            >  3:1488656/1‑149 (MQ=255)
                               tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                         <  2:882815/149‑1 (MQ=255)
                                 aaTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAtt                                                                                                       >  1:1024840/1‑149 (MQ=255)
                                   tGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                                                                      >  4:994577/1‑148 (MQ=255)
                                   tGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                                                                      <  3:994577/148‑1 (MQ=255)
                                           tCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCt                                                                                             >  2:1807302/1‑149 (MQ=255)
                                            cGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCtt                                                                                            <  4:1488656/149‑1 (MQ=255)
                                                 cGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCt                                                                                       >  2:1015348/1‑149 (MQ=255)
                                                  gAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAt                                                                                                        >  1:1343722/1‑131 (MQ=255)
                                                  gAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAt                                                                                                        <  2:1343722/131‑1 (MQ=255)
                                                   aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                <  1:1063131/122‑1 (MQ=255)
                                                   aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                >  2:1063131/1‑122 (MQ=255)
                                                       aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                >  1:994263/1‑118 (MQ=255)
                                                       aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                <  2:994263/118‑1 (MQ=255)
                                                       aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcac                                                                                                  >  1:1062239/1‑132 (MQ=255)
                                                       aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcac                                                                                                  <  2:1062239/132‑1 (MQ=255)
                                                       aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGa                                                                                 >  4:356015/1‑149 (MQ=255)
                                                        ttttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaat                                                                                                                         <  4:1426474/108‑1 (MQ=255)
                                                        ttttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaat                                                                                                                         >  3:1426474/1‑108 (MQ=255)
                                                            gccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGc                                                                            <  1:1015348/149‑1 (MQ=255)
                                                             ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATCAAAAGCa                                                                           <  2:1024840/149‑1 (MQ=255)
                                                                   ccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                        >  3:335552/1‑82 (MQ=255)
                                                                   ccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                        <  4:335552/82‑1 (MQ=255)
                                                                   ccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGa                                                                     >  1:1529593/1‑149 (MQ=255)
                                                                    cagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGc                                                                            >  3:718044/1‑141 (MQ=255)
                                                                    cagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGc                                                                            <  4:718044/141‑1 (MQ=255)
                                                                     agtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtctc                                                                                                    <  1:1410106/116‑1 (MQ=255)
                                                                     agtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtctc                                                                                                    >  2:1410106/1‑116 (MQ=255)
                                                                           tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGt                                                             <  1:1807302/149‑1 (MQ=255)
                                                                           tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGt                                                             <  2:1529593/149‑1 (MQ=255)
                                                                            cTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTa                                                            <  1:1692170/149‑1 (MQ=255)
                                                                             tGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCatt                                                                                                                                        <  4:1783167/72‑3 (MQ=255)
                                                                             tGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCatt                                                                                                                                        >  3:1783167/1‑70 (MQ=255)
                                                                               aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                         >  1:1095992/1‑149 (MQ=255)
                                                                               aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                         >  4:628879/1‑149 (MQ=255)
                                                                               aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                         >  3:1194493/1‑149 (MQ=255)
                                                                               aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                         <  3:628879/149‑1 (MQ=255)
                                                                                 cTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAAt                                                                                   >  3:1794050/1‑121 (MQ=255)
                                                                                 cTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAAt                                                                                   <  4:1794050/121‑1 (MQ=255)
                                                                                  tAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                                                                          >  2:1403560/1‑129 (MQ=255)
                                                                                  tAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                                                                          <  1:1403560/129‑1 (MQ=255)
                                                                                   aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGc                                                      <  1:242160/148‑1 (MQ=255)
                                                                                   aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGc                                                      >  2:242160/1‑148 (MQ=255)
                                                                                    aTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTc                                                    <  3:607563/149‑1 (MQ=255)
                                                                                             gCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                >  4:608920/1‑128 (MQ=255)
                                                                                             gCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                <  3:608920/128‑1 (MQ=255)
                                                                                                 tCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCa                                       >  1:1119402/1‑149 (MQ=255)
                                                                                                     acTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCt                                   >  2:368106/1‑149 (MQ=255)
                                                                                                       tctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCGCAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCg                                                        >  1:1374926/1‑126 (MQ=255)
                                                                                                       tctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCg                                                        <  2:1374926/126‑1 (MQ=255)
                                                                                                       tctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCa                                 >  2:1628567/1‑149 (MQ=255)
                                                                                                         tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaa                                                                                <  3:656517/100‑1 (MQ=255)
                                                                                                         tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaa                                                                                >  4:656517/1‑100 (MQ=255)
                                                                                                             cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaa                                                                                >  3:868158/1‑96 (MQ=255)
                                                                                                             cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaa                                                                                <  4:868158/96‑1 (MQ=255)
                                                                                                             cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGc                           <  2:1095992/149‑1 (MQ=255)
                                                                                                                ccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtctc                                                                                                    <  1:1671724/73‑1 (MQ=255)
                                                                                                                ccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtctc                                                                                                    >  2:1671724/1‑73 (MQ=255)
                                                                                                                ccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCata                        >  2:592211/1‑149 (MQ=255)
                                                                                                                 cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAtt                                                                                                       <  4:1350098/69‑1 (MQ=255)
                                                                                                                 cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAtt                                                                                                       >  3:1350098/1‑69 (MQ=255)
                                                                                                                 cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                                                              <  1:836639/110‑1 (MQ=255)
                                                                                                                 cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                                                              >  2:836639/1‑110 (MQ=255)
                                                                                                                   gtgcgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                         >  2:1740655/3‑129 (MQ=255)
                                                                                                                   gcgcGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                         <  1:1740655/129‑1 (MQ=255)
                                                                                                                      cgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTg                  <  1:1628567/149‑1 (MQ=255)
                                                                                                                       gACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTgg                 <  4:1194493/149‑1 (MQ=255)
                                                                                                                          aGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTa              >  3:674774/1‑149 (MQ=255)
                                                                                                                           gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                     <  3:1170754/61‑1 (MQ=255)
                                                                                                                           gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                     >  4:1170754/1‑61 (MQ=255)
                                                                                                                             gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTagag           >  4:1117605/1‑149 (MQ=255)
                                                                                                                             gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTACTCTCACAACCtt                                                                                            <  2:1285370/68‑1 (MQ=255)
                                                                                                                             gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTACTCTCACAACCtt                                                                                            >  1:1285370/1‑68 (MQ=255)
                                                                                                                               aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                >  3:469720/1‑94 (MQ=255)
                                                                                                                               aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                <  4:469720/94‑1 (MQ=255)
                                                                                                                               aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCa         <  1:592211/149‑1 (MQ=255)
                                                                                                                                 tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa       <  4:674774/149‑1 (MQ=255)
                                                                                                                                 tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa       <  1:368106/149‑1 (MQ=255)
                                                                                                                                  gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAg      >  4:797381/1‑149 (MQ=255)
                                                                                                                                  cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa       >  1:173714/4‑148 (MQ=255)
                                                                                                                                  cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa       <  2:173714/145‑1 (MQ=255)
                                                                                                                                   attGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCat                         >  2:1496047/3‑129 (MQ=255)
                                                                                                                                   attGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCat                         <  1:1496047/127‑1 (MQ=255)
                                                                                                                                   aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                                                                   <  1:1649746/55‑1 (MQ=255)
                                                                                                                                   aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                                                                   >  2:1649746/1‑55 (MQ=255)
                                                                                                                                   aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa       <  2:811274/147‑1 (MQ=255)
                                                                                                                                   aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa       >  1:811274/1‑147 (MQ=255)
                                                                                                                                   aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa     <  3:1117605/149‑1 (MQ=255)
                                                                                                                                    aTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCt                                                     >  4:1440200/1‑100 (MQ=255)
                                                                                                                                    aTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCt                                                     <  3:1440200/100‑1 (MQ=255)
                                                                                                                                      ggCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAt  >  2:1865213/1‑149 (MQ=255)
                                                                                                                                      ggCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAt  >  2:769820/1‑149 (MQ=255)
                                                                                                                                           gcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCatat                       <  1:743507/123‑1 (MQ=255)
                                                                                                                                           gcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCatat                       >  2:743507/1‑123 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAT  >  NC_000913/2172212‑2172494

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