Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,350 A→T 5.7% L175Q (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,3500AT 5.7% 303.1 / 7.4 87L175Q (CTG→CAG) gatDgalactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (44/38);  new base T (3/2);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

TAAAACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTT  >  NC_000913/2172202‑2172496
                                                                                                                                                    |                                                                                                                                                  
tAAAACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCca                                                                                                                                                    >  4:852187/1‑149 (MQ=255)
tAAAACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCca                                                                                                                                                    >  3:57613/1‑149 (MQ=255)
  aaaCGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagc                                                                                                                                                  >  3:497445/1‑149 (MQ=255)
       ctctGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGc                                                                                                                                              >  2:800140/1‑148 (MQ=255)
        tctGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCa                                                                                                                                            >  4:72272/1‑149 (MQ=255)
              tttGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcag                                                                                                                                       <  1:800140/148‑6 (MQ=255)
               ttGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGtt                                                                                                                                      <  2:6509/148‑1 (MQ=255)
                        gcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGATCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                                                                            >  4:625733/1‑149 (MQ=255)
                         cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                                                                            >  1:839561/1‑148 (MQ=255)
                             tttCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataata                                                                                                                       <  3:852187/149‑1 (MQ=255)
                               tCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataataa                                                                                                                      >  1:487981/1‑148 (MQ=255)
                                cGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAataataata                                                                                                                    <  4:299924/149‑1 (MQ=255)
                                 gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATCGTTCCGGCACCAAtaataataa                                                                                                                   >  3:554245/1‑149 (MQ=255)
                                   tACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACg                                                                                                                 >  3:503480/1‑149 (MQ=255)
                                     cTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                          >  1:322082/1‑122 (MQ=255)
                                     cTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                          <  2:322082/122‑1 (MQ=255)
                                     cTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtt                                                                                                               <  4:385351/149‑1 (MQ=255)
                                       gTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                              >  2:171897/1‑148 (MQ=255)
                                         tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCa                                                                                                                                            <  2:629590/116‑1 (MQ=255)
                                         tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCa                                                                                                                                            >  1:629590/1‑116 (MQ=255)
                                         tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttttt                                                                                                            >  2:795616/1‑148 (MQ=255)
                                         tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttttt                                                                                                            >  1:652235/1‑148 (MQ=255)
                                         tAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttttt                                                                                                            >  2:273484/1‑148 (MQ=255)
                                           agtGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAt                                                                                                          >  2:687825/3‑148 (MQ=255)
                                                   cATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcacac                                                                                                  <  1:795616/148‑2 (MQ=255)
                                                   cATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                                  >  2:457385/1‑148 (MQ=255)
                                                      cGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCAGACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCt                                                                                               <  1:687825/148‑1 (MQ=255)
                                                      cGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCt                                                                                               >  2:65822/1‑148 (MQ=255)
                                                        cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTg                                                                                             <  2:211412/148‑1 (MQ=255)
                                                             aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTcc                                                                                                                                    <  4:546371/104‑1 (MQ=255)
                                                             aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTcc                                                                                                                                    >  3:546371/1‑104 (MQ=255)
                                                              aaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                                                                      <  4:630462/149‑1 (MQ=255)
                                                                   tttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaa                                                                                  >  2:287704/1‑148 (MQ=255)
                                                                    ttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaaaa                                                                                <  4:497445/149‑1 (MQ=255)
                                                                      gccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGTCAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGc                                                                              <  3:625733/149‑1 (MQ=255)
                                                                         agtgcgagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCCTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGCTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                           <  4:585665/149‑1 (MQ=255)
                                                                          gtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                           <  2:384160/148‑1 (MQ=255)
                                                                          gtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCGCACAACCTTGCGCTAAATGAAAAGCATg                                                                           >  2:884481/1‑148 (MQ=255)
                                                                              cagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGa                                                                       <  1:884481/148‑1 (MQ=255)
                                                                              cagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGa                                                                       <  1:273484/148‑1 (MQ=255)
                                                                                  ttttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCaa                                                                   <  2:839561/148‑1 (MQ=255)
                                                                                   tttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                  <  1:536761/148‑1 (MQ=255)
                                                                                   tttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                  <  1:486255/148‑1 (MQ=255)
                                                                                       tGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGACACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTaa                                                             >  4:829526/1‑149 (MQ=255)
                                                                                        gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCTCCAATAATATTAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTaa                                                             <  1:171897/148‑1 (MQ=255)
                                                                                        gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTGTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTaa                                                             >  2:593454/1‑148 (MQ=255)
                                                                                         aaCTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                           >  3:896182/1‑149 (MQ=255)
                                                                                             aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                                                                            <  2:89611/128‑1 (MQ=255)
                                                                                             aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                                                                            >  1:89611/1‑128 (MQ=255)
                                                                                             aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGc                                                        <  2:487981/148‑1 (MQ=255)
                                                                                              aTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTc                                                      <  3:286863/149‑1 (MQ=255)
                                                                                                gTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa                                                    >  4:809873/1‑149 (MQ=255)
                                                                                                 tCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAAt                                                   >  4:290936/1‑149 (MQ=255)
                                                                                                     tCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaaa                                               >  3:816938/1‑149 (MQ=255)
                                                                                                     tCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                >  2:803961/1‑148 (MQ=255)
                                                                                                     tCGACGTCACACTCTTTGCTCCCACCGCGACAGCGCATTTAACTGCCAGCAGGCCAATGGTTCCGACACCAATAATCATAACGTTTTTATTCTCACATCCTTTCGCTAAATGAAAGTCATGCAGACCAACGGTAATCGGGACAATaaa                                                <  2:421732/148‑1 (MQ=255)
                                                                                                        ccGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                             <  1:503990/148‑1 (MQ=255)
                                                                                                        ccGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                             >  2:503990/1‑148 (MQ=255)
                                                                                                         cGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                           <  3:402911/149‑1 (MQ=255)
                                                                                                         cGTCACACTCGTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCCATGGTTCCGGCACCGATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAAAGAAAAGCATGCAGACCCACGGTAATCGGCTCAATAAAAGccc                                           >  4:786328/1‑149 (MQ=255)
                                                                                                             acacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATc                                       >  4:336217/1‑149 (MQ=255)
                                                                                                               acTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                           >  4:369200/1‑143 (MQ=255)
                                                                                                               acTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                           <  3:369200/143‑1 (MQ=255)
                                                                                                                ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCt                                     >  1:270264/1‑148 (MQ=255)
                                                                                                                 tctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCa                                   >  3:496474/1‑149 (MQ=255)
                                                                                                                 tctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCa                                   >  3:846100/1‑149 (MQ=255)
                                                                                                                  ctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCa                                   >  2:392006/1‑148 (MQ=255)
                                                                                                                   tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCGCCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTATTCGGCTCAATAAAAGCCCCATCCTCAAt                                 <  4:554245/149‑1 (MQ=255)
                                                                                                                           cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCata                          >  2:909249/1‑148 (MQ=255)
                                                                                                                                 gACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTg                    <  1:803961/148‑1 (MQ=255)
                                                                                                                                   cAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa                                                    >  1:444140/1‑114 (MQ=255)
                                                                                                                                   cAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa                                                    <  2:444140/114‑1 (MQ=255)
                                                                                                                                     gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTag               >  3:89359/1‑149 (MQ=255)
                                                                                                                                        cACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCAt                                        <  1:152121/121‑1 (MQ=255)
                                                                                                                                        cACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCAt                                        >  2:152121/1‑121 (MQ=255)
                                                                                                                                         aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTc                                                      <  3:218568/106‑1 (MQ=255)
                                                                                                                                         aCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTc                                                      >  4:218568/1‑106 (MQ=255)
                                                                                                                                          cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaaa                                               >  1:89462/1‑112 (MQ=255)
                                                                                                                                          cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaaa                                               <  2:89462/112‑1 (MQ=255)
                                                                                                                                          cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                >  2:427747/1‑111 (MQ=255)
                                                                                                                                          cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                <  1:427747/111‑1 (MQ=255)
                                                                                                                                           tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa         <  3:336217/149‑1 (MQ=255)
                                                                                                                                            gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                >  2:154175/1‑109 (MQ=255)
                                                                                                                                            gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                <  1:154175/109‑1 (MQ=255)
                                                                                                                                            gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa         >  1:604818/1‑148 (MQ=255)
                                                                                                                                            gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAg        >  4:112363/1‑149 (MQ=255)
                                                                                                                                             aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa       <  4:485501/149‑1 (MQ=255)
                                                                                                                                              ttGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa       >  1:218703/2‑148 (MQ=255)
                                                                                                                                                 gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAtt   <  3:72272/149‑1 (MQ=255)
                                                                                                                                                  ccAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAttt  >  3:890429/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TAAAACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTT  >  NC_000913/2172202‑2172496

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