Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp pseudogene (915‑916/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.97.5% 326.0 / 1.4 81pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (40/39);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.97.5% 325.2 / 1.0 81pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (40/39);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173222‑2173495
                                                                                                                                              ||                                                                                                                                     
tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                        <  2:1336127/140‑1 (MQ=255)
tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                        <  1:987451/140‑1 (MQ=255)
ttcccACGAGCCGCGGCACCCGCCAGCGCCACGAAGAAGCCGATGGTGGAAAGTTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCTGGC‑TTGCCGATaccac                                                                                                                                        <  2:416948/137‑1 (MQ=255)
 ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATaccacc                                                                                                                                       >  1:494390/1‑140 (MQ=255)
  tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                      <  2:309800/140‑2 (MQ=255)
  tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                      <  2:1055630/140‑2 (MQ=255)
  tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                      <  1:1889264/140‑2 (MQ=255)
         tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATg                                                                                                                             >  1:234007/1‑140 (MQ=255)
                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                   >  2:1801481/1‑140 (MQ=255)
                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                   >  1:809730/1‑140 (MQ=255)
                   ccGCCATCGCCACGAAGAAGCCGATGG‑GGGAAGATCCCCAAACCGGAACCCCTTGTTACCCCGGACACAAACAGGAA‑TTTAAATTG‑GGGGGGGGGTAAAAATTAAGCTTTGCCAAT‑CCACCCC‑CCGAATG‑CCCAacagca                                                                                                                  >  2:59421/1‑140 (MQ=255)
                        aTCGCCGCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                              <  1:374438/140‑1 (MQ=255)
                        aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                              <  1:503763/140‑1 (MQ=255)
                         tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                             <  1:1253665/140‑1 (MQ=255)
                             cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCTGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                         <  2:281323/140‑1 (MQ=255)
                                  agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TT‑AAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                   <  1:1726335/140‑1 (MQ=255)
                                   gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                   <  1:1640475/140‑1 (MQ=255)
                                   gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                   <  1:1697355/140‑1 (MQ=255)
                                     agCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                 >  1:221008/1‑140 (MQ=255)
                                     agCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                 <  1:1541218/140‑1 (MQ=255)
                                      gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                <  2:359643/140‑1 (MQ=255)
                                      gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                <  2:1329156/140‑1 (MQ=255)
                                          gttgtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                            <  1:59421/137‑1 (MQ=255)
                                           tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                           <  2:349971/140‑1 (MQ=255)
                                           tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                           >  1:349971/1‑140 (MQ=255)
                                                  aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGGGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                    <  2:846847/140‑1 (MQ=255)
                                                       cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                               >  1:2001046/1‑140 (MQ=255)
                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                              >  2:1235795/1‑140 (MQ=255)
                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                              <  1:1244987/140‑1 (MQ=255)
                                                         ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                             >  2:788072/1‑140 (MQ=255)
                                                              cGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                        >  2:1710119/1‑140 (MQ=255)
                                                                 cagcaCCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGTAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                     >  1:1378510/1‑140 (MQ=255)
                                                                  agcaCCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                    <  1:764088/140‑1 (MQ=255)
                                                                    caCCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGTAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                  <  2:1378510/140‑1 (MQ=255)
                                                                    caCCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                  <  2:802566/140‑1 (MQ=255)
                                                                     cccTGATTACCCGGCACACAAACAGAAA‑TTAAAATGG‑TGAGGGGGAAAAAAACCAGGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                 <  2:856189/139‑1 (MQ=255)
                                                                         gATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                             >  2:992241/1‑140 (MQ=255)
                                                                          atcaCCCGGCACACAAACACCAA‑TTAAAATGG‑TGAGTGGAATAAAAACCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                            <  2:494390/137‑1 (MQ=255)
                                                                          aTTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                            <  2:2039014/140‑1 (MQ=255)
                                                                          aTTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                            <  1:341627/140‑1 (MQ=255)
                                                                          aTTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                            <  1:1814881/140‑1 (MQ=255)
                                                                               ccGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                       >  2:1417045/1‑140 (MQ=255)
                                                                                 ggCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                     >  2:50763/1‑140 (MQ=255)
                                                                                 ggCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                     >  2:636222/1‑140 (MQ=255)
                                                                                  gCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                    <  1:1656539/140‑1 (MQ=255)
                                                                                   cacacaAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                   >  2:1929901/1‑140 (MQ=255)
                                                                                   cacacaAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                   >  1:1502021/1‑140 (MQ=255)
                                                                                    acacaAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                  <  1:1682830/140‑1 (MQ=255)
                                                                                    acacaAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                  >  2:1130470/1‑140 (MQ=255)
                                                                                    acacaAACAGCAA‑TTAAAATGG‑TGAGGGGGATAAAAATCGGGC‑TCGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                  <  2:392671/140‑1 (MQ=255)
                                                                                    acacaAAC‑GCAATTTAAAGTGG‑GGAGGGGAAAAAAAACCAGCT‑TTCCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                  <  2:2050935/140‑1 (MQ=255)
                                                                                    acaaaAACAGCAA‑TTAAAATGG‑TGAGGGGGATAAAAATCAGGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                  <  2:1326087/140‑1 (MQ=255)
                                                                                         aaCAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                             >  1:1430437/1‑140 (MQ=255)
                                                                                         aaCAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                             >  1:897744/1‑140 (MQ=255)
                                                                                             gCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                         >  2:1113781/1‑140 (MQ=255)
                                                                                               aa‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                       <  1:788072/140‑1 (MQ=255)
                                                                                                 ataaaaaTGG‑TGAGGGGGAAAAAAACCAGGC‑TTGCCAAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                      <  2:1296541/137‑1 (MQ=255)
                                                                                                   tAAAATGGTTGAGTGGAATAAAAGTCAGGC‑TTGCGGAT‑CC‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                     <  2:85688/140‑1 (MQ=255)
                                                                                                   tAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                    >  1:194989/1‑140 (MQ=255)
                                                                                                     aaaTGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                  >  1:380960/1‑140 (MQ=255)
                                                                                                      aaTGG‑GGAGTGGGATAAAATCCGGGC‑TTGCCGATGCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                 <  1:1801481/140‑1 (MQ=255)
                                                                                                       aTGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                >  2:289350/1‑140 (MQ=255)
                                                                                                       aTGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAAGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                >  1:1070410/1‑140 (MQ=255)
                                                                                                         gg‑TGAGTGGGGTAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                              >  2:1554340/1‑140 (MQ=255)
                                                                                                            tGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                            <  1:636222/140‑1 (MQ=255)
                                                                                                             gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                           >  1:1407637/1‑140 (MQ=255)
                                                                                                             gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                           >  2:1230010/1‑140 (MQ=255)
                                                                                                             gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                           >  2:879982/1‑140 (MQ=255)
                                                                                                              aGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                          >  1:804185/1‑140 (MQ=255)
                                                                                                               gTGGGAAAAAAATCAGCC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                         <  2:1407637/140‑1 (MQ=255)
                                                                                                                tGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                        >  1:1439848/1‑140 (MQ=255)
                                                                                                                tGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                        >  2:1797730/1‑140 (MQ=255)
                                                                                                                tGGGATAAAAATCAGGC‑TTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                        <  2:1154071/140‑1 (MQ=255)
                                                                                                                 gggATAAAAATCAGGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                       <  1:171659/140‑1 (MQ=255)
                                                                                                                  ggATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                      >  2:949898/1‑140 (MQ=255)
                                                                                                                     tAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                   <  1:1113781/140‑1 (MQ=255)
                                                                                                                     tAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                   <  1:1974596/140‑1 (MQ=255)
                                                                                                                     aaaaaaTCAGCT‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                   <  2:809730/139‑1 (MQ=255)
                                                                                                                     aaaaaaCCAGGC‑TTCCCATTGCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                   <  1:1154224/139‑1 (MQ=255)
                                                                                                                      aaaaaTCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                  <  1:1797730/140‑1 (MQ=255)
                                                                                                                      aaaaaTCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                  <  1:1929901/140‑1 (MQ=255)
                                                                                                                          aTCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg              >  2:1385149/1‑140 (MQ=255)
                                                                                                                          aTCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg              >  2:534799/1‑140 (MQ=255)
                                                                                                                          aTCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg              >  2:1985984/1‑140 (MQ=255)
                                                                                                                             aGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt           >  2:1824935/1‑140 (MQ=255)
                                                                                                                             aGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt           >  1:1329266/1‑140 (MQ=255)
                                                                                                                              ggC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCAATGATGGGTTTa          >  2:1283421/1‑140 (MQ=255)
                                                                                                                              ggC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa         <  1:457599/140‑1 (MQ=255)
                                                                                                                                c‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt        <  1:1283421/140‑1 (MQ=255)
                                                                                                                                   tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa      >  2:1215050/1‑140 (MQ=255)
                                                                                                                                   tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCAGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa      >  1:1255237/1‑140 (MQ=255)
                                                                                                                                      cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc   <  1:1896261/140‑1 (MQ=255)
                                                                                                                                      cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc   <  1:1235795/140‑1 (MQ=255)
                                                                                                                                      cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc   <  2:2001046/140‑1 (MQ=255)
                                                                                                                                      cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc   <  1:1130470/140‑1 (MQ=255)
                                                                                                                                         gatacccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg  <  2:1005346/135‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                     
TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGG‑TGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173222‑2173495

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