Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,169 T→G D76A (GAT→GCT)  metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,1690TG98.5% 261.2 / ‑3.9 67D76A (GAT→GCT) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (31/35);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

AGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGA  >  NC_000913/4128022‑4128314
                                                                                                                                                   |                                                                                                                                                 
aGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCg                                                                                                                                                 >  2:244636/1‑150 (MQ=255)
 ggCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGc                                                                                                                                               >  1:868537/1‑151 (MQ=255)
   ccTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcat                                                                                                                                             >  2:992443/1‑151 (MQ=255)
     tGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcatc                                                                                                                                            >  1:744465/1‑150 (MQ=255)
     tGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcatc                                                                                                                                            >  1:686783/1‑150 (MQ=255)
     tGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcatc                                                                                                                                            >  1:475712/1‑150 (MQ=255)
        tAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATcc                                                                                                                                        >  2:761686/1‑151 (MQ=255)
          aGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCgg                                                                                                                                      <  1:231957/151‑1 (MQ=255)
          aGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCgg                                                                                                                                      <  2:529007/151‑1 (MQ=255)
           gCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGc                                                                                                                                     <  1:94704/151‑1 (MQ=255)
               aGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAg                                                                                                                                 <  2:1083591/151‑1 (MQ=255)
                   cATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGtt                                                                                                                              <  1:440081/150‑1 (MQ=255)
                   cATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTg                                                                                                                             <  1:716207/151‑1 (MQ=255)
                      cAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGccc                                                                                                                          >  2:746014/1‑151 (MQ=255)
                      cAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGcc                                                                                                                           <  1:532859/150‑1 (MQ=255)
                      cAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGcc                                                                                                                           <  1:319606/150‑1 (MQ=255)
                        ggCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCgg                                                                                                                        >  2:767853/1‑151 (MQ=255)
                         gCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGt                                                                                                                       <  2:400150/151‑1 (MQ=255)
                         gCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGt                                                                                                                       <  1:348361/151‑1 (MQ=255)
                            aTTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTaaa                                                                                                                    <  2:1030887/151‑1 (MQ=255)
                             ttCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTaa                                                                                                                     >  2:533250/1‑149 (MQ=255)
                               ccACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGc                                                                                                                 >  2:958996/1‑151 (MQ=255)
                                    ccGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATg                                                                                                              <  1:918724/149‑1 (MQ=255)
                                    ccGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCa                                                                                                            <  1:81372/151‑1 (MQ=255)
                                    ccGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCa                                                                                                            <  2:886350/151‑1 (MQ=255)
                                           gcTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACg                                                                                                      >  2:11932/1‑150 (MQ=255)
                                           gcTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGc                                                                                                     >  1:868987/1‑151 (MQ=255)
                                                 ttttGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGc                                                                                                <  1:742515/150‑1 (MQ=255)
                                                       tttAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCAcagc                                                                                           >  2:827864/1‑149 (MQ=255)
                                                         tAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACagcagc                                                                                        >  2:946579/1‑150 (MQ=255)
                                                           gTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCg                                                                                     >  2:311588/1‑151 (MQ=255)
                                                             aTTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCt                                                                                   <  2:835315/151‑1 (MQ=255)
                                                              ttCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTACAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTg                                                                                  >  1:749919/1‑151 (MQ=255)
                                                                   aCGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGt                                                                             >  1:619888/1‑151 (MQ=255)
                                                                          cGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTg                                                                       >  2:1048376/1‑150 (MQ=255)
                                                                               tAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCa                                                                  <  2:698262/150‑1 (MQ=255)
                                                                                aaTCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAgg                                                                >  2:46908/1‑151 (MQ=255)
                                                                                           cACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCa                                                         <  2:869032/147‑1 (MQ=255)
                                                                                                aTGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacgac                                                >  2:314707/1‑151 (MQ=255)
                                                                                                     ctctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtg                                           >  1:987098/1‑151 (MQ=255)
                                                                                                       ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcg                                         >  1:175706/1‑151 (MQ=255)
                                                                                                         tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcg                                         <  1:999683/149‑1 (MQ=255)
                                                                                                           tGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTc                                      <  1:437325/150‑1 (MQ=255)
                                                                                                           tGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTc                                      <  2:526827/150‑1 (MQ=255)
                                                                                                            gCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCAt                                    <  1:533216/151‑1 (MQ=255)
                                                                                                            gCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCAt                                    >  1:876205/1‑151 (MQ=255)
                                                                                                            gCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCAt                                    <  1:654580/151‑1 (MQ=255)
                                                                                                               gcctCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCgg                                 <  1:573296/148‑1 (MQ=255)
                                                                                                                   ccGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAg                             <  2:723134/151‑1 (MQ=255)
                                                                                                                   ccGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAg                             <  2:752623/151‑1 (MQ=255)
                                                                                                                    cGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg                            >  1:638157/1‑151 (MQ=255)
                                                                                                                    cGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg                            <  1:903803/151‑1 (MQ=255)
                                                                                                                       gATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtt                         >  1:998082/1‑151 (MQ=255)
                                                                                                                       gATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtt                         >  2:225822/1‑151 (MQ=255)
                                                                                                                        aTTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttt                        >  2:977210/1‑151 (MQ=255)
                                                                                                                         tttCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                       >  1:97769/1‑151 (MQ=255)
                                                                                                                          ttCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                       <  1:339479/150‑1 (MQ=255)
                                                                                                                             gtcgCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAac                   <  2:894471/151‑1 (MQ=255)
                                                                                                                               cgCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAacac                 >  1:928544/1‑151 (MQ=255)
                                                                                                                               cgCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAacac                 <  2:293205/151‑1 (MQ=255)
                                                                                                                                   gcgcTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt              <  2:492922/150‑1 (MQ=255)
                                                                                                                                   gcgcTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCt               <  1:340671/149‑1 (MQ=255)
                                                                                                                                     gcTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag           <  2:991650/151‑1 (MQ=255)
                                                                                                                                        ctTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAaga          <  1:362605/149‑1 (MQ=255)
                                                                                                                                            aCGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAAt     <  2:1026294/150‑1 (MQ=255)
                                                                                                                                            aCGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg    >  2:256291/1‑151 (MQ=255)
                                                                                                                                               cAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa  >  2:808964/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGA  >  NC_000913/4128022‑4128314

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