Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,139,639 T→G 6.2% intergenic (+46/+77) yghU → / ← hybG putative S‑transferase/hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,139,6390TG 6.2% 373.2 / 11.5 119intergenic (+46/+77)yghU/hybGputative S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base T (58/53);  new base G (3/5);  total (61/58)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 121 total aligned reads displayed.

GCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCG  >  NC_000913/3139491‑3139785
                                                                                                                                                    |                                                                                                                                                  
gCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGtt                                                                                                                                                    <  3:142252/149‑1 (MQ=255)
      tGTTAACCGCACCAACGGACCGCTGAATCAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTAGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                              >  2:145234/1‑149 (MQ=255)
        ttAACCGCACCAACGGCCCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCCGCAAGCCATCCAGGCCGGATAAGGCGTTCacgacg                                                                                                                                             >  4:96633/1‑148 (MQ=255)
           aCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                          <  7:397784/148‑5 (MQ=255)
           aCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                          <  5:245599/148‑1 (MQ=255)
           aCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGGTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                         <  1:286330/149‑1 (MQ=255)
           aCCGCACCAACGGACCGCTGAATGAGCAGGTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                          <  1:146973/148‑5 (MQ=255)
              gCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                      <  7:362847/149‑1 (MQ=255)
              gCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGCTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                      <  7:151567/149‑1 (MQ=255)
                   aaCGGACCGCAGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACTAATACGGAAGATAACCGTCAGGGGTAAGGGGTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATc                                                                                                                                 >  4:171830/1‑149 (MQ=255)
                    aCGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                              <  6:77233/135‑1 (MQ=255)
                    aCGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                              >  5:77233/1‑135 (MQ=255)
                     cGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGt                                                                                                                               <  8:260859/149‑1 (MQ=255)
                       gACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTg                                                                                                                              >  5:319540/1‑148 (MQ=255)
                            cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAAt                                                                                                                        <  8:306980/149‑1 (MQ=255)
                                aTGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                                                                                    >  2:96634/1‑149 (MQ=255)
                                    gCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGAt                                                                                                                 <  7:255716/148‑1 (MQ=255)
                                     cAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGc                                                                                                               >  8:415409/1‑149 (MQ=255)
                                     cAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGc                                                                                                               >  2:43946/1‑149 (MQ=255)
                                          gCATGAGCTCCATGACGCCAGTGATTTCTAGACGAATACGGAAGATCAGCGTCATGTGTAAGGGTTTTTGTTCTTCTCAGCAAGCCATCCAGTCTGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACg                                                                                                           <  6:339374/148‑1 (MQ=37)
                                            aTGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                                                        >  4:250133/1‑149 (MQ=255)
                                             tGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                                                                                                  <  8:258331/122‑1 (MQ=255)
                                               aGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCTTCTCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACATTGCCTGTTGCGACGCTGcc                                                                                                      <  4:282415/148‑1 (MQ=255)
                                               aGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                                                                                                      <  2:155831/148‑1 (MQ=255)
                                                    cATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCt                                                                                                >  7:16803/1‑149 (MQ=255)
                                                      tGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                              <  2:81456/149‑1 (MQ=255)
                                                             aGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGcc                                                                                       <  3:159290/149‑1 (MQ=255)
                                                                 aTTTCGAGACGAATCCGGAAGATAAGCGTCAGGGGGAAGGGGGGGGGGGCGTCGCAGCAAGCCATCCAGGCCGGAGAGGGCGTGCACGCCGCAGGCGGCAATCGTGCCGCTTGCCTGATGCGCCGCTGCGGCGTGTTTTCAGGCCtgc                                                                                    >  8:305817/1‑146 (MQ=255)
                                                                  tttCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                        >  2:366099/1‑95 (MQ=255)
                                                                  tttCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                        <  1:366099/95‑1 (MQ=255)
                                                                  tttCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCCGGCCTACg                                                                                   >  4:163499/1‑148 (MQ=255)
                                                                  tttCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGGTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                        <  1:372063/95‑1 (MQ=255)
                                                                         aCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGATGCCGCGGCTTATCAGGCCTACGAAATCTaa                                                                           >  8:334612/1‑149 (MQ=255)
                                                                           gAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                          >  3:353999/1‑148 (MQ=255)
                                                                                 ggAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcg                                                                   >  3:361438/1‑149 (MQ=255)
                                                                                   aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctg                                                                                                                                  <  1:15640/84‑4 (MQ=21)
                                                                                   aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                                                                                    <  5:236201/98‑1 (MQ=255)
                                                                                   aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                                                                                    >  6:236201/1‑98 (MQ=255)
                                                                                   aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGAAAGCCATCCAGTCCGGATAAGGCGTTCACGCCGCATACGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgcg                                                                 >  1:310236/1‑149 (MQ=255)
                                                                                     gATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt                                                               <  1:43946/149‑1 (MQ=255)
                                                                                       tAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                                                             <  1:99146/149‑1 (MQ=255)
                                                                                         aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                              <  5:7767/114‑1 (MQ=255)
                                                                                         aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                              >  6:7767/1‑114 (MQ=255)
                                                                                          gCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                                                          >  2:237518/1‑149 (MQ=255)
                                                                                           cGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAAt                                                         <  4:353999/149‑1 (MQ=255)
                                                                                            gTCAGGGGTAAGGGTTGGTGGTCGTCGCAGCAAGCCATCCAGGCCGGTTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                                                        <  7:126102/149‑1 (MQ=255)
                                                                                             tCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                                                       >  3:21380/1‑149 (MQ=255)
                                                                                               aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTc                                                      <  3:250133/148‑1 (MQ=255)
                                                                                               aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTc                                                      <  3:98558/148‑1 (MQ=255)
                                                                                                ggggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGt                                                    <  6:247898/149‑1 (MQ=255)
                                                                                                  ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcatc                                                                      >  7:373048/1‑129 (MQ=255)
                                                                                                  ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcatc                                                                      <  8:373048/129‑1 (MQ=255)
                                                                                                  ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAc                                                                    <  5:252516/131‑1 (MQ=255)
                                                                                                  ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAc                                                                    >  6:252516/1‑131 (MQ=255)
                                                                                                  ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTaa                                                  <  8:129268/149‑1 (MQ=255)
                                                                                                         gtgggTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCATTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAttt                                           <  8:186890/146‑1 (MQ=255)
                                                                                                          ttGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTg                                          <  7:328379/149‑1 (MQ=255)
                                                                                                          ttGGTGTTCGTCACAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTg                                          <  3:170230/149‑1 (MQ=255)
                                                                                                               gtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcg                                       >  4:24769/1‑147 (MQ=255)
                                                                                                                ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTAc                                                                                    <  2:203418/101‑1 (MQ=255)
                                                                                                                 tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGt                                   >  6:94810/1‑149 (MQ=255)
                                                                                                                  cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGt                                   <  5:133524/148‑1 (MQ=255)
                                                                                                                      gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCTGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                              >  2:272406/1‑149 (MQ=255)
                                                                                                                      gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                              >  1:151275/1‑149 (MQ=255)
                                                                                                                      gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                              <  1:237518/149‑1 (MQ=255)
                                                                                                                      gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                              >  3:322445/1‑149 (MQ=255)
                                                                                                                      gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCGGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                              <  3:171830/149‑1 (MQ=255)
                                                                                                                        agcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                                                     <  3:43535/124‑1 (MQ=255)
                                                                                                                        agcaAGCCATCCACGCCGGATAAGGCGTTCACCCCGCATCCGGCAAACCTGAACAATGCCTGATGCCACGCTGCCGCGACTTATCAGGCCTACGAAAGCTAATCATCACGCCCAGGTAATGTCg                                                     >  4:43535/1‑124 (MQ=38)
                                                                                                                          caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAt                          <  2:310236/149‑1 (MQ=255)
                                                                                                                          caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAt                          <  6:61666/149‑1 (MQ=255)
                                                                                                                           aaGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAtg                         <  4:296602/149‑1 (MQ=255)
                                                                                                                             gCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                                                                >  5:55359/1‑108 (MQ=255)
                                                                                                                             gCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAtgtg                       >  6:351414/1‑149 (MQ=255)
                                                                                                                              ccATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                      >  1:156054/1‑149 (MQ=255)
                                                                                                                              ccATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                      >  7:196284/1‑149 (MQ=255)
                                                                                                                               cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                              <  3:154602/76‑1 (MQ=255)
                                                                                                                               cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                              >  4:154602/1‑76 (MQ=255)
                                                                                                                                aTCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggct                    >  4:340374/1‑149 (MQ=255)
                                                                                                                                aTCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggct                    >  3:4190/1‑149 (MQ=255)
                                                                                                                                  gattgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                      >  8:275721/5‑145 (MQ=38)
                                                                                                                                  gattgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                      <  7:275721/141‑1 (MQ=38)
                                                                                                                                  ccAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATcc                                               >  3:335160/1‑120 (MQ=255)
                                                                                                                                  ccAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATcc                                               <  4:335160/120‑1 (MQ=255)
                                                                                                                                   cAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggcttt                  <  2:151275/148‑1 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttgg                >  6:328056/1‑149 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttgg                >  7:297479/1‑149 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggct              <  4:322445/149‑1 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggct              >  7:396787/1‑149 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggctt             <  3:221459/149‑1 (MQ=255)
                                                                                                                                        cGGATAAGGCGTGCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCCTCAAGCGCTGGTAATTTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTc            <  6:344035/149‑1 (MQ=255)
                                                                                                                                          gATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCAt          >  5:171086/1‑149 (MQ=255)
                                                                                                                                           aTAAGGCGTTCTCGCGGCATCCGGCAATCGTGCACAATGCGTGATGCGACGCTGCCGCGTCTTATAAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATc         <  1:145234/149‑1 (MQ=255)
                                                                                                                                             aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCtt       >  2:55105/1‑149 (MQ=255)
                                                                                                                                              aGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtc      <  8:16803/149‑1 (MQ=255)
                                                                                                                                                gCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtcgt    >  8:340671/1‑149 (MQ=255)
                                                                                                                                                 cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtcgtc   <  3:24769/149‑1 (MQ=255)
                                                                                                                                                 agTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCCCGCTCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGGCTAATGTGGCTTTGGCTTCATCTtcgtc   <  4:42369/148‑1 (MQ=255)
                                                                                                                                                  gTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtcgtc   >  6:373493/1‑148 (MQ=255)
                                                                                                                                                    tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTtcgtcg  <  5:328056/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCG  >  NC_000913/3139491‑3139785

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: