Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,171,904 C→G 7.3% E324Q (GAA→CAA)  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,171,9040CG 7.3% 327.4 / 10.9 93E324Q (GAA→CAA) gatDgalactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (48/38);  new base G (2/5);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATT  >  NC_000913/2171764‑2172049
                                                                                                                                            |                                                                                                                                                 
cGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTc                                                                                                                                           <  4:1261372/149‑1 (MQ=255)
cGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTc                                                                                                                                           <  1:1860499/149‑1 (MQ=255)
                gTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGGGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGa                                                                                                                                   <  3:1028469/141‑1 (MQ=255)
                gTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCCCGGTGa                                                                                                                                   >  4:1028469/1‑141 (MQ=255)
                 tGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTc                                                                                                                                           <  3:1777753/132‑1 (MQ=255)
                 tGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTc                                                                                                                                           >  4:1777753/1‑132 (MQ=255)
                         cTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAg                                                                                                                  <  1:677676/149‑1 (MQ=255)
                                  acaATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTc                                                                                                                                           <  3:1891482/115‑1 (MQ=255)
                                  acaATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTc                                                                                                                                           >  4:1891482/1‑115 (MQ=255)
                                              cAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAg                                                                                             >  2:67427/1‑149 (MQ=255)
                                              cAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAg                                                                                             >  1:1678536/1‑149 (MQ=255)
                                               aGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAAc                                                                                                           >  4:296844/1‑134 (MQ=255)
                                               aGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAAc                                                                                                           <  3:296844/134‑1 (MQ=255)
                                                  aTCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAAc                                                                                         >  3:105577/1‑149 (MQ=255)
                                                   tCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAAcc                                                                                        >  1:754416/1‑149 (MQ=255)
                                                     aCGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATgg                                                                                                                        <  1:180357/115‑1 (MQ=255)
                                                     aCGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATgg                                                                                                                        >  2:180357/1‑115 (MQ=255)
                                                       gCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTcgc                                                                                                                                         >  2:1303585/1‑94 (MQ=255)
                                                       gCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTcgc                                                                                                                                         <  1:1303585/96‑3 (MQ=255)
                                                       gCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCt                                                                                    >  4:1542036/1‑149 (MQ=255)
                                                             ccGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGt                                                                              >  2:927610/1‑149 (MQ=255)
                                                                   tttCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCa                                                                        >  3:46876/1‑149 (MQ=255)
                                                                       aGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCtt                                                                                                                                            <  4:829648/77‑1 (MQ=255)
                                                                       aGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCtt                                                                                                                                            >  3:829648/1‑77 (MQ=255)
                                                                       aGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTAc                                                                                                       <  1:1339624/114‑1 (MQ=255)
                                                                       aGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTAc                                                                                                       >  2:1339624/1‑114 (MQ=255)
                                                                       aGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTcc                                                                    >  4:1052080/1‑149 (MQ=255)
                                                                       aGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTcc                                                                    >  4:29928/1‑149 (MQ=255)
                                                                          gAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCTGCTCGCCGTTTCCCACTCCAGc                                                                 >  1:1612973/1‑149 (MQ=255)
                                                                             tGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTAc                                                                                                       >  4:141076/1‑108 (MQ=255)
                                                                             tGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTAc                                                                                                       >  4:1793260/1‑108 (MQ=255)
                                                                             tGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTAc                                                                                                       <  3:1793260/108‑1 (MQ=255)
                                                                             tGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTc                                                                     >  1:1738749/1‑142 (MQ=255)
                                                                             tGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTc                                                                     <  2:1738749/142‑1 (MQ=255)
                                                                             tGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCAACGGTGAGCGATTAATGGCTCCAGGCTTAACTTAc                                                                                                       <  3:141076/108‑1 (MQ=255)
                                                                              gAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCg                                                             <  2:1678536/149‑1 (MQ=255)
                                                                                gCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGc                                                           <  3:1052080/149‑1 (MQ=255)
                                                                                 cAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCa                                                                                              <  2:1754414/113‑1 (MQ=255)
                                                                                 cAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCa                                                                                              >  1:1754414/1‑113 (MQ=255)
                                                                                  aaCACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCa                                                         >  1:3860/1‑149 (MQ=255)
                                                                                       tttGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGc                                                    >  3:811247/1‑149 (MQ=255)
                                                                                       tttGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGc                                                    >  1:690328/1‑149 (MQ=255)
                                                                                       tttGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGc                                                    >  4:597253/1‑149 (MQ=255)
                                                                                         tGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTg                                                  <  4:46876/149‑1 (MQ=255)
                                                                                         tGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTg                                                  <  1:927610/149‑1 (MQ=255)
                                                                                           cccGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGa                                                <  2:754416/149‑1 (MQ=255)
                                                                                           cccGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGa                                                <  3:1542036/149‑1 (MQ=255)
                                                                                             cGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGagt                                              <  1:67427/149‑1 (MQ=255)
                                                                                              ggCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGagta                                             <  4:811247/149‑1 (MQ=255)
                                                                                              ggCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGagta                                             <  4:105577/149‑1 (MQ=255)
                                                                                                cATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTCTCCCACTCCTGCCCCGGCCAAGGGCTGGagtagt                                           >  1:479240/1‑149 (MQ=255)
                                                                                                 aTAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGtt                                          >  4:575296/1‑149 (MQ=255)
                                                                                                 aTAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGtt                                          >  1:485172/1‑149 (MQ=255)
                                                                                                 aTAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGtt                                          >  3:1824984/1‑149 (MQ=255)
                                                                                                  tAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGt                                                                              >  4:1510561/1‑112 (MQ=255)
                                                                                                  tAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGt                                                                              <  3:1510561/112‑1 (MQ=255)
                                                                                                  tAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTc                                         <  3:29928/149‑1 (MQ=255)
                                                                                                     cATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTCTGCACGCTTCCAAGGTGAGCGATTAATAGCTCCAGGCTTAACTTACGTTCTGTCAGCAGCCGTCTCGCCGTTTCCCCCTCCTGCCCCGGCCAAGGCTTTGTGGAGTTCATc                                      <  2:3860/149‑1 (MQ=255)
                                                                                                      aTTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAAt                                                                                                                          <  2:950039/64‑1 (MQ=255)
                                                                                                      aTTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAAt                                                                                                                          >  1:950039/1‑64 (MQ=255)
                                                                                                      aTTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATcc                                     >  2:1195337/1‑149 (MQ=255)
                                                                                                        tACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCaa                                                        >  3:1846830/1‑128 (MQ=255)
                                                                                                        tACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCaa                                                        <  4:1846830/128‑1 (MQ=255)
                                                                                                             gCGATGTCGCGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTg                                                  >  3:373612/1‑129 (MQ=255)
                                                                                                             gCGATGTCGCGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTg                                                  <  4:373612/129‑1 (MQ=255)
                                                                                                                 tGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTc                                         >  3:1329124/1‑134 (MQ=255)
                                                                                                                 tGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTc                                         <  4:1329124/134‑1 (MQ=255)
                                                                                                                   tCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTACCTTACGTTCt                                                                                                  <  4:1228672/75‑1 (MQ=255)
                                                                                                                   tCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCt                                                                                                  >  3:1228672/1‑75 (MQ=255)
                                                                                                                   tCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAAc                        >  3:582962/1‑149 (MQ=255)
                                                                                                                     aCGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCAt                                       >  3:1125179/1‑132 (MQ=255)
                                                                                                                     aCGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCAt                                       <  4:1125179/132‑1 (MQ=255)
                                                                                                                        cACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCa                   <  2:485172/149‑1 (MQ=255)
                                                                                                                          ccGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGc                               <  1:1870728/135‑1 (MQ=255)
                                                                                                                          ccGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGc                               >  2:1870728/1‑135 (MQ=255)
                                                                                                                           cGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGct                >  2:432949/1‑149 (MQ=255)
                                                                                                                              cTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCtt             <  3:597253/149‑1 (MQ=255)
                                                                                                                              cTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCtt             <  4:582962/149‑1 (MQ=255)
                                                                                                                                 ggCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCa                                                                                                                   <  2:560844/44‑1 (MQ=255)
                                                                                                                                 ggCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCa                                                                                                                   >  1:560844/1‑44 (MQ=255)
                                                                                                                                   cGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGcc                              <  3:909853/127‑1 (MQ=255)
                                                                                                                                   cGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGcc                              >  4:909853/1‑127 (MQ=255)
                                                                                                                                     aaGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAggg                                                     <  2:1566585/102‑1 (MQ=255)
                                                                                                                                     aaGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAggg                                                     >  1:1566585/1‑102 (MQ=255)
                                                                                                                                       gCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAgg                                                      <  4:509017/99‑1 (MQ=255)
                                                                                                                                       gCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAga                                                      >  3:509017/1‑98 (MQ=255)
                                                                                                                                       gCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAata    >  1:360471/1‑149 (MQ=255)
                                                                                                                                         tttCAAATCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATAtt  >  2:1031539/1‑149 (MQ=255)
                                                                                                                                         tttCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATc                                      >  2:19573/1‑113 (MQ=255)
                                                                                                                                         tttCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATc                                      <  1:19573/113‑1 (MQ=255)
                                                                                                                                         tttCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATAtt  >  4:340391/1‑149 (MQ=255)
                                                                                                                                          ttCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTAc          >  2:1164421/1‑140 (MQ=255)
                                                                                                                                          ttCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTAc          <  1:1164421/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATT  >  NC_000913/2171764‑2172049

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