Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,139,632 A→G 5.2% intergenic (+39/+84) 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,6320AG 5.2% 365.7 / 7.6 115intergenic (+39/+84)yghU/hybGputative S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base A (56/53);  new base G (2/4);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 117 total aligned reads displayed.

ACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTT  >  NC_000913/3139485‑3139780
                                                                                                                                                   |                                                                                                                                                    
aCGGGGGCGTATTGTTAACCGCACCAACGGACCGCTGACTGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATaa                                                                                                                                                     <  4:267235/149‑1 (MQ=255)
  gTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAgg                                                                                                                                                   >  8:310049/1‑149 (MQ=255)
      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)
                                          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)
                                                gCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGc                                                                                                     >  2:99146/1‑149 (MQ=255)
                                                 cATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTTGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                                                                                    >  4:51764/1‑149 (MQ=255)
                                                  aTGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                                                   >  4:332721/1‑149 (MQ=255)
                                                  aTGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                                                   >  4:250133/1‑149 (MQ=255)
                                                   tGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                                                                                             <  8:258331/122‑1 (MQ=255)
                                                   tGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                                                                                             >  7:258331/1‑122 (MQ=255)
                                                     aGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCTTCTCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACATTGCCTGTTGCGACGCTGcc                                                                                                 <  4:282415/148‑1 (MQ=255)
                                                     aGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                                                                                                 <  2:155831/148‑1 (MQ=255)
                                                         ccATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTc                                                                                            >  8:328379/1‑149 (MQ=255)
                                                          cATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCt                                                                                           >  7:16803/1‑149 (MQ=255)
                                                            tGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                         <  2:81456/149‑1 (MQ=255)
                                                                gCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAg                                                                                     >  5:61666/1‑149 (MQ=255)
                                                                gCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAg                                                                                     <  6:319540/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)
                                                                        tttCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGGTGGTGGTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                   >  2:372063/1‑95 (MQ=255)
                                                                               aCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGATGCCGCGGCTTATCAGGCCTACGAAATCTaa                                                                      >  8:334612/1‑149 (MQ=255)
                                                                                 gAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                     >  3:353999/1‑148 (MQ=255)
                                                                                 gAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                     >  1:430208/1‑148 (MQ=255)
                                                                                         aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctg                                                                                                                             >  2:15640/1‑81 (MQ=21)
                                                                                         aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctg                                                                                                                             <  1:15640/84‑4 (MQ=21)
                                                                                         aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                                                                               <  5:236201/98‑1 (MQ=255)
                                                                                         aaGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGAAAGCCATCCAGTCCGGATAAGGCGTTCACGCCGCATACGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgcg                                                            >  1:310236/1‑149 (MQ=255)
                                                                                           gATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt                                                          <  1:43946/149‑1 (MQ=255)
                                                                                               aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                         >  6:7767/1‑114 (MQ=255)
                                                                                               aGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                         <  5:7767/114‑1 (MQ=255)
                                                                                                gCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                                                     >  2:237518/1‑149 (MQ=255)
                                                                                                 cGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAAt                                                    <  4:353999/149‑1 (MQ=255)
                                                                                                  gTCAGGGGTAAGGGTTGGTGGTCGTCGCAGCAAGCCATCCAGGCCGGTTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                                                   <  7:126102/149‑1 (MQ=255)
                                                                                                     aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTc                                                 <  3:98558/148‑1 (MQ=255)
                                                                                                     aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTc                                                 <  3:250133/148‑1 (MQ=255)
                                                                                                      ggggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGt                                               <  6:247898/149‑1 (MQ=255)
                                                                                                        ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcatc                                                                 <  8:373048/129‑1 (MQ=255)
                                                                                                        ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcatc                                                                 >  7:373048/1‑129 (MQ=255)
                                                                                                        ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAc                                                               >  6:252516/1‑131 (MQ=255)
                                                                                                        ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAc                                                               <  5:252516/131‑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)
                                                                                                                      ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTAc                                                                               >  1:203418/1‑101 (MQ=255)
                                                                                                                       tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGt                              >  6:94810/1‑149 (MQ=255)
                                                                                                                        cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGt                              <  5:133524/148‑1 (MQ=255)
                                                                                                                            gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         >  3:322445/1‑149 (MQ=255)
                                                                                                                            gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         >  1:151275/1‑149 (MQ=255)
                                                                                                                            gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCGGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         <  3:171830/149‑1 (MQ=255)
                                                                                                                              agcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                                                <  3:43535/124‑1 (MQ=255)
                                                                                                                                caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAt                     <  6:61666/149‑1 (MQ=255)
                                                                                                                                caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAt                     <  3:332721/149‑1 (MQ=255)
                                                                                                                                caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAt                     <  2:310236/149‑1 (MQ=255)
                                                                                                                                  aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAtgt                   <  4:21380/149‑1 (MQ=255)
                                                                                                                                  aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAtgt                   >  8:247639/1‑149 (MQ=255)
                                                                                                                                   gCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                                                           >  5:55359/1‑108 (MQ=255)
                                                                                                                                   gCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                                                           <  6:55359/108‑1 (MQ=255)
                                                                                                                                    ccATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                 >  1:156054/1‑149 (MQ=255)
                                                                                                                                    ccATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                 >  7:196284/1‑149 (MQ=255)
                                                                                                                                     cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                         >  4:154602/1‑76 (MQ=255)
                                                                                                                                     cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                         <  3:154602/76‑1 (MQ=255)
                                                                                                                                      aTCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggct               >  3:4190/1‑149 (MQ=255)
                                                                                                                                      aTCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggct               >  4:340374/1‑149 (MQ=255)
                                                                                                                                        gattgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                 <  7:275721/141‑1 (MQ=38)
                                                                                                                                        gattgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtgg                 >  8:275721/5‑145 (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)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggct         >  7:396787/1‑149 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggct         >  3:6918/1‑149 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggct         <  4:322445/149‑1 (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)
                                                                                                                                                  tAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCt   <  7:393658/149‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCtt  >  2:55105/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTT  >  NC_000913/3139485‑3139780

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: