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% 332.3 / 8.3 96L174Q (CTG→CAG) gatDgalactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (44/46);  new base T (2/4);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTC  >  NC_000913/2172205‑2172500
                                                                                                                                                    |                                                                                                                                                   
aaCGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagca                                                                                                                                                     >  1:822582/1‑149 (MQ=255)
       tGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                              <  2:890004/149‑1 (MQ=255)
         cATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattca                                                                                                                                            <  4:47994/149‑5 (MQ=255)
         cATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattca                                                                                                                                            <  4:642485/149‑5 (MQ=255)
         cATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATgg                                                                                                                                            >  3:187285/1‑149 (MQ=255)
           tttGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGtt                                                                                                                                          >  4:256368/1‑149 (MQ=255)
            ttGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGCTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcagtg                                                                                                                                         <  2:770228/149‑8 (MQ=255)
             tGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTcc                                                                                                                                        <  4:532438/149‑1 (MQ=255)
                       tCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACAGTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaata                                                                                                                              >  4:768983/1‑149 (MQ=255)
                       tCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTTCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaata                                                                                                                              >  1:171245/1‑149 (MQ=255)
                         aTTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataat                                                                                                                            <  2:699799/149‑1 (MQ=255)
                          tttCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataata                                                                                                                           <  4:663777/149‑1 (MQ=255)
                              gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAAtaataataa                                                                                                                       >  3:234513/1‑149 (MQ=255)
                              gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAAtaataataa                                                                                                                       <  2:232739/149‑1 (MQ=255)
                                tACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACg                                                                                                                     >  1:120951/1‑149 (MQ=255)
                                tACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACg                                                                                                                     >  3:835828/1‑149 (MQ=255)
                                  cTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtt                                                                                                                   <  1:581559/149‑1 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                                  >  3:377076/1‑149 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                                  <  4:508959/149‑1 (MQ=255)
                                     ttAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttttt                                                                                                                >  4:136562/1‑149 (MQ=255)
                                      tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                               <  3:381616/149‑1 (MQ=255)
                                         aTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAttg                                                                                                            >  1:827452/1‑148 (MQ=255)
                                         aTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                                                                            <  4:724057/149‑1 (MQ=255)
                                         aTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                                                                            <  4:293559/149‑1 (MQ=255)
                                          tGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                           <  1:623611/149‑1 (MQ=255)
                                           gTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtctc                                                                                                          >  4:138088/1‑149 (MQ=255)
                                               gCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                                      <  2:351513/149‑1 (MQ=255)
                                                     cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGCATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgc                                                                                                <  4:251930/149‑1 (MQ=255)
                                                        cGAAAGATTTTGTCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                      >  2:18588/1‑124 (MQ=255)
                                                        cGAAAGATTTTGTCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                      <  1:18588/124‑1 (MQ=255)
                                                        cGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttttt                                                                                                                <  2:571401/130‑1 (MQ=255)
                                                        cGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttttt                                                                                                                >  1:571401/1‑130 (MQ=255)
                                                          aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaa                                                                                           >  1:524091/1‑149 (MQ=255)
                                                          aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaa                                                                                           >  2:464097/1‑149 (MQ=255)
                                                           aaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                                                                          <  3:778929/149‑1 (MQ=255)
                                                           aaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                                                                          >  4:647798/1‑149 (MQ=255)
                                                           aaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                                                                          <  2:259829/149‑1 (MQ=255)
                                                           aaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                                                                          <  4:377076/149‑1 (MQ=255)
                                                               ttttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaa                                                                                      >  4:397497/1‑149 (MQ=255)
                                                                tttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaaa                                                                                     <  3:83481/149‑1 (MQ=255)
                                                                    ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCa                                                                                 <  3:852878/149‑1 (MQ=255)
                                                                     cagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                           >  3:390394/1‑122 (MQ=255)
                                                                     cagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                           <  4:390394/122‑1 (MQ=255)
                                                                      agtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGTCAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                               >  3:438531/1‑149 (MQ=255)
                                                                          ccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGa                                                                           >  4:638048/1‑149 (MQ=255)
                                                                           cagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAc                                                                          <  3:758710/149‑1 (MQ=255)
                                                                            agcttttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAcc                                                                         >  3:175828/4‑149 (MQ=255)
                                                                              tttttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCaa                                                                       >  1:763953/1‑149 (MQ=255)
                                                                              tttttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCaa                                                                       <  4:119394/149‑1 (MQ=255)
                                                                               ttttCTTAACTAATGTCGATCGCCGTCACACCCCTTGCTCCCAGCGCGACAGCGCAATGAAAGGCCAGCAGGCCAATGGGTCCGGCCCCAATAATTATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                      <  3:491436/149‑1 (MQ=255)
                                                                                  tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGt                                                                   <  3:199480/149‑1 (MQ=255)
                                                                                     gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGTCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAAt                                                                <  2:370815/149‑1 (MQ=255)
                                                                                      aaCTAATTTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                               >  3:370476/1‑149 (MQ=255)
                                                                                      aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                               >  4:622159/1‑149 (MQ=255)
                                                                                      aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                               >  2:131311/1‑149 (MQ=255)
                                                                                        cTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg                                                             >  1:926678/1‑149 (MQ=255)
                                                                                         tAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGc                                                            >  2:497061/1‑149 (MQ=255)
                                                                                           aTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGACACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTc                                                          <  3:638048/149‑1 (MQ=255)
                                                                                              tCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAAt                                                       <  2:871495/149‑1 (MQ=255)
                                                                                                 aTCGCCGTCACACTCTTTGCTCCCAGCGCGCCAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                    >  2:301941/1‑149 (MQ=255)
                                                                                                 aTCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCTGCTCAATaaa                                                    >  1:92951/1‑149 (MQ=255)
                                                                                                  tCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaaa                                                   <  4:234513/149‑1 (MQ=255)
                                                                                                  tCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCCGCTCAATaaaa                                                   <  3:525768/149‑1 (MQ=255)
                                                                                                    gCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                                 <  1:111010/149‑1 (MQ=255)
                                                                                                       gTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                                  >  3:487855/1‑81 (MQ=255)
                                                                                                       gTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                                  <  4:487855/81‑1 (MQ=255)
                                                                                                          acacTCTTTGCTCCCTGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGATAATCGGCTCAATAAAAGCCCCATc                                           >  2:627590/1‑149 (MQ=255)
                                                                                                          acacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcgc                                                                                              >  1:829887/1‑98 (MQ=255)
                                                                                                          acacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcgc                                                                                              <  2:829887/98‑1 (MQ=255)
                                                                                                            acTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAGATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCt                                         >  2:770689/1‑149 (MQ=255)
                                                                                                              tctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCa                                       >  1:90929/1‑149 (MQ=255)
                                                                                                                tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAAt                                     >  3:698341/1‑149 (MQ=255)
                                                                                                                tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAAt                                     >  3:188477/1‑149 (MQ=255)
                                                                                                                 ttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATa                                    <  1:464097/149‑1 (MQ=255)
                                                                                                                    cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGc                                 <  2:776271/149‑1 (MQ=255)
                                                                                                                       ccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCata                              >  3:627008/1‑149 (MQ=255)
                                                                                                                         aGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATc                            >  1:257532/1‑149 (MQ=255)
                                                                                                                         aGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATc                            >  4:645556/1‑149 (MQ=255)
                                                                                                                         aGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATc                            >  2:636711/1‑149 (MQ=255)
                                                                                                                         aGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATc                            >  2:877765/1‑149 (MQ=255)
                                                                                                                             cgtcAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTg                        <  4:438531/146‑1 (MQ=255)
                                                                                                                             cgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTg                        <  3:136562/149‑1 (MQ=255)
                                                                                                                              gACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTgg                       <  4:302358/149‑1 (MQ=255)
                                                                                                                                    gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTagag                 <  3:228709/149‑1 (MQ=255)
                                                                                                                                      aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCa               >  2:318219/1‑149 (MQ=255)
                                                                                                                                      aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCa               >  4:857235/1‑149 (MQ=255)
                                                                                                                                        tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             <  3:647798/149‑1 (MQ=255)
                                                                                                                                        tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             <  3:397497/149‑1 (MQ=255)
                                                                                                                                        tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             <  2:926678/149‑1 (MQ=255)
                                                                                                                                        tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             <  1:636711/149‑1 (MQ=255)
                                                                                                                                        tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAACCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             <  3:622159/149‑1 (MQ=255)
                                                                                                                                         cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATcc                           <  2:33259/131‑1 (MQ=255)
                                                                                                                                         cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATcc                           >  1:33259/4‑134 (MQ=255)
                                                                                                                                         cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaaa            >  4:178406/4‑148 (MQ=255)
                                                                                                                                          aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa           >  2:33184/1‑149 (MQ=255)
                                                                                                                                           aTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGAc          <  1:493527/149‑1 (MQ=255)
                                                                                                                                             ggCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAt        <  4:188477/149‑1 (MQ=255)
                                                                                                                                                   cagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTc  <  4:355675/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTC  >  NC_000913/2172205‑2172500

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