Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 284,450 A→G E417G (GAG→GGG)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913284,4500AG100.0% 393.3 / NA 106E417G (GAG→GGG) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (61/45);  total (61/45)

Only 100 of 106 total aligned reads displayed.

GTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGT  >  NC_000913/284302‑284597
                                                                                                                                                    |                                                                                                                                                   
gTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGc                                                                                                                                                   <  2:461444/151‑1 (MQ=255)
gTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGc                                                                                                                                                   <  2:409750/151‑1 (MQ=255)
  gATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGcc                                                                                                                                                  <  1:598224/150‑1 (MQ=255)
     gCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGAt                                                                                                                                              <  2:829503/151‑1 (MQ=255)
      cTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGAt                                                                                                                                              <  1:1193693/150‑1 (MQ=255)
       tCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATg                                                                                                                                             >  2:995477/1‑150 (MQ=255)
        ccACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGAc                                                                                                                                           >  1:567775/1‑151 (MQ=255)
            cTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTg                                                                                                                                        >  1:1100774/1‑150 (MQ=255)
              gcgcGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGa                                                                                                                                       <  2:382927/149‑1 (MQ=255)
              gcgcGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGa                                                                                                                                       <  4:436885/149‑1 (MQ=255)
              gcgcGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGGAGGGCCGGATGACGTGATc                                                                                                                                     >  4:187733/1‑151 (MQ=255)
                 cgACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTg                                                                                                                                  <  2:1089747/151‑1 (MQ=255)
                  gACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTgc                                                                                                                                 >  2:1132236/1‑151 (MQ=255)
                       cGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTgccga                                                                                                                              <  2:81172/149‑2 (MQ=255)
                        ggCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                           >  3:175094/1‑151 (MQ=255)
                         gCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGga                                                                                                                          <  4:375196/151‑1 (MQ=255)
                          ccTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                           >  2:169758/1‑149 (MQ=255)
                           ctgctgCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                           >  2:862349/1‑148 (MQ=255)
                           ctgctgCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                           >  1:908223/1‑148 (MQ=255)
                            tgctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGgaga                                                                                                                        >  4:189632/1‑150 (MQ=255)
                              ctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAgg                                                                                                                     >  2:865114/1‑151 (MQ=255)
                                 cATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCaa                                                                                                                  <  4:394413/151‑1 (MQ=255)
                                 cATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCATAGGGCCGGCTGACGTGATCCTGCCGCCGGAGACGGCa                                                                                                                   >  3:144753/1‑150 (MQ=255)
                                   tCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCaaa                                                                                                                 >  3:354244/1‑150 (MQ=255)
                                     tGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCaaaa                                                                                                                >  1:278318/1‑149 (MQ=255)
                                      ggACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGc                                                                                                              <  1:651723/150‑1 (MQ=255)
                                        aCGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGa                                                                                                            >  2:305944/1‑150 (MQ=255)
                                         cGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGaaa                                                                                                          <  1:929589/151‑1 (MQ=255)
                                         cGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGaa                                                                                                           >  1:661832/1‑150 (MQ=255)
                                          gCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGaa                                                                                                           >  1:1261117/1‑149 (MQ=255)
                                          gCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAg                                                                                                         <  1:1043167/151‑1 (MQ=255)
                                            cATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAggg                                                                                                       <  2:289950/151‑1 (MQ=255)
                                               gaccgtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTg                                                                                                    >  2:1151674/1‑151 (MQ=255)
                                                     gaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCg                                                                                              >  2:841046/1‑151 (MQ=255)
                                                     gaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCg                                                                                              >  2:1094557/1‑151 (MQ=255)
                                                     gaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCg                                                                                              >  1:182367/1‑151 (MQ=255)
                                                     gaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCAAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCg                                                                                              <  4:419471/151‑1 (MQ=255)
                                                       ccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGAc                                                                                            >  2:921402/1‑151 (MQ=255)
                                                          gCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACgg                                                                                          >  2:1058638/1‑150 (MQ=255)
                                                           ccAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTc                                                                                        >  3:16101/1‑151 (MQ=255)
                                                           ccAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTc                                                                                        <  4:403860/151‑1 (MQ=255)
                                                                cGGTGGGCGCGAAACTTGAATGGTGGCAGGCGTACGAGCGCCGGGCGCGCTTCCGCCAGAGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAACGGGCTGACCTCGTCGGTCTGCtt                                                                                   >  2:887273/1‑151 (MQ=255)
                                                                cGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCtt                                                                                   <  1:435101/151‑1 (MQ=255)
                                                                 ggTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCtt                                                                                   <  2:793161/150‑1 (MQ=255)
                                                                  gTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTcc                                                                                 >  1:891223/1‑151 (MQ=255)
                                                                     ggCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGa                                                                              <  2:195427/151‑1 (MQ=255)
                                                                       cGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACg                                                                            >  3:445741/1‑151 (MQ=255)
                                                                         agaACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACggg                                                                          >  1:132458/1‑151 (MQ=255)
                                                                          gaACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACggg                                                                          <  2:1141687/150‑1 (MQ=255)
                                                                           aaCCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGc                                                                         >  2:13338/1‑150 (MQ=255)
                                                                             ccTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAAc                                                                      >  2:1112342/1‑151 (MQ=255)
                                                                             ccTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAAc                                                                      >  3:378587/1‑151 (MQ=255)
                                                                               ttGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACAt                                                                    >  4:458088/1‑151 (MQ=255)
                                                                               ttGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACAt                                                                    >  1:778112/1‑151 (MQ=255)
                                                                                 gAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATc                                                                   <  1:1066262/150‑1 (MQ=255)
                                                                                  aaTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCCGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATc                                                                   >  1:664524/1‑149 (MQ=255)
                                                                                  aaTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCg                                                                  <  2:660296/150‑1 (MQ=255)
                                                                                  aaTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGc                                                                 <  1:640031/151‑1 (MQ=255)
                                                                                   aTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCt                                                                <  2:282008/151‑1 (MQ=255)
                                                                                     ggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTcc                                                              <  2:209009/151‑1 (MQ=255)
                                                                                      gtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCg                                                             >  1:212563/1‑151 (MQ=255)
                                                                                      gtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCg                                                             <  1:430348/151‑1 (MQ=255)
                                                                                               gTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCg                                                    >  1:583399/1‑151 (MQ=255)
                                                                                                 ccGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCgg                                                   >  2:1100413/1‑150 (MQ=255)
                                                                                                 ccGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGt                                                  <  2:148983/151‑1 (MQ=255)
                                                                                                  tgAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                                 >  2:1008706/2‑151 (MQ=255)
                                                                                                  cGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                                 <  1:774721/151‑1 (MQ=255)
                                                                                                     gcgcCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATc                                              >  4:358924/1‑151 (MQ=255)
                                                                                                      cgcCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCa                                             <  2:527561/151‑1 (MQ=255)
                                                                                                       gcCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCa                                             <  2:853553/150‑1 (MQ=255)
                                                                                                        ccGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCa                                             >  3:284205/1‑149 (MQ=255)
                                                                                                             cgcgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCAc                                      <  1:370017/151‑1 (MQ=255)
                                                                                                               cgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACgg                                    <  2:198128/151‑1 (MQ=255)
                                                                                                               cgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACgg                                    <  1:793702/151‑1 (MQ=255)
                                                                                                               cgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACg                                     <  1:952931/150‑1 (MQ=255)
                                                                                                                 cTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGc                                   <  3:50547/150‑1 (MQ=255)
                                                                                                                  ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                  >  2:224481/1‑150 (MQ=255)
                                                                                                                  ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                  >  2:1038108/1‑150 (MQ=255)
                                                                                                                    ccGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATc                               >  2:1282236/1‑151 (MQ=255)
                                                                                                                       ccAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAc                            >  1:226562/1‑151 (MQ=255)
                                                                                                                         aGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAccc                          <  1:180759/151‑1 (MQ=255)
                                                                                                                           ttccagcgcGAGCAGGCCGCCGTAGGGCCGGATGACGTAATCCTGCCGCCGGAGTAGGCAAAACCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGACCCGGGAGGCTCGGTGATAAAGGCCACGGCGATCgtgcc                          >  3:148740/6‑145 (MQ=255)
                                                                                                                            gCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTc                       <  2:278278/151‑1 (MQ=255)
                                                                                                                              cTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCgg                     >  2:1249268/1‑151 (MQ=255)
                                                                                                                                gcgcGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCgg                     <  4:287923/149‑1 (MQ=255)
                                                                                                                                  gcgAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtggt                 >  1:870808/1‑151 (MQ=255)
                                                                                                                                  gcgAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtggt                 <  3:438318/151‑1 (MQ=255)
                                                                                                                                   cgAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGgtggtg                >  1:1209827/1‑151 (MQ=255)
                                                                                                                                   cgAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGgtggtg                >  1:814327/1‑151 (MQ=255)
                                                                                                                                       cAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCg            >  2:119154/1‑151 (MQ=255)
                                                                                                                                        aGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa           <  4:319615/151‑1 (MQ=255)
                                                                                                                                         ggACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa           <  2:279443/150‑1 (MQ=255)
                                                                                                                                          gACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa           >  4:29805/1‑149 (MQ=255)
                                                                                                                                           aCGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa           >  2:1061347/1‑148 (MQ=255)
                                                                                                                                            cGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa           >  3:6217/1‑147 (MQ=255)
                                                                                                                                             ggCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATg      >  1:595728/1‑151 (MQ=255)
                                                                                                                                               cGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGc    >  2:186169/1‑151 (MQ=255)
                                                                                                                                               cGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGc    >  2:201970/1‑151 (MQ=255)
                                                                                                                                               cGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGc    >  1:14671/1‑151 (MQ=255)
                                                                                                                                                 tAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGt  >  1:145087/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGT  >  NC_000913/284302‑284597

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