Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,432,583 T→G F58C (TTC→TGC)  tfaR → Rac prophage; putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,432,5830TG100.0% 105.0 / NA 57F58C (TTC→TGC) tfaRRac prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (32/25);  total (32/25)

GTGAACAACCACGGACCATAAAAATTTATAATCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAAACA  >  NC_000913/1432427‑1432813
                                                                                                                                                            |                                                                                                                                                                                                                                      
gTGAACAACCACGGACCATAAAAATTTATAATCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTACAACGGGTTATGACGTGGCTTCCGGCaa                                                                                                                                                             >  1:691076/1‑232 (MQ=11)
     caacCACGGACCATAAAAATTTATAATCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGtaatt                                                                                                                                                    >  1:308147/1‑236 (MQ=14)
                  tAAAAATTTATAATCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTg                                                                                                                                            >  1:655105/1‑231 (MQ=14)
                  tAAAAATTTATAATCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACt                                                                                                                                        <  1:342968/235‑1 (MQ=14)
                               tCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCgg                                                                                                                          >  2:614055/1‑236 (MQ=14)
                                 tgctgGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGaa                                                                                                                        <  1:741844/236‑1 (MQ=14)
                                      gCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATg                                                                                                                    >  2:282694/1‑235 (MQ=14)
                                       ccGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATGCAGGTTGGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGaaaa                                                                                                                      <  1:678051/232‑1 (MQ=9)
                                                     tttATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGcc                                                                                                    >  2:419031/1‑236 (MQ=14)
                                                          tGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAg                                                                                               >  2:302012/1‑236 (MQ=14)
                                                             tGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgggg                                                                                            <  1:600868/236‑1 (MQ=14)
                                                             tGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgggg                                                                                            <  2:689116/236‑1 (MQ=14)
                                                               aaGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAg                                                                                          >  2:71811/1‑236 (MQ=14)
                                                               aaGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATAGCCGGAAgggta                                                                                           >  1:730170/1‑233 (MQ=14)
                                                                      tGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAgc                                                                                   >  1:617185/1‑235 (MQ=14)
                                                                         atttatTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATAGTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                                                                                 <  1:711632/232‑1 (MQ=14)
                                                                         atatatTCCGCCTCATCCAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                                                                                 <  1:528560/235‑1 (MQ=14)
                                                                            tatTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGGTATAGCCGGAAGGGGAGTTTCTGTAGTGGa                                                                             >  2:342560/1‑236 (MQ=9)
                                                                                 cGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCCGGAGTGGAACTGc                                                                        >  2:397745/1‑236 (MQ=9)
                                                                                 cGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGc                                                                        <  1:545481/236‑1 (MQ=14)
                                                                                 cGCCTAATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGATTCATCGTGACATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGc                                                                        <  1:13170/236‑1 (MQ=9)
                                                                                  gCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTTTCGCCGGAAGGGGAGTTTTCGAAGTGGAACGGc                                                                        >  2:16152/1‑235 (MQ=9)
                                                                                  gCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGca                                                                       <  1:419031/236‑1 (MQ=14)
                                                                                   ccTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGTTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGcac                                                                      >  2:702973/1‑236 (MQ=14)
                                                                                   ccTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTGACCGGAAAATGTTACCTGGTTATCGCCGGCAGGGGAGTTTCAGAAGTGGAAAGGcaa                                                                      >  2:435062/1‑235 (MQ=9)
                                                                                      cATGCAGGTCTGCCAGGAAGCAGTACCGATATTGCGCCGCCAGATTGTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGc                                                                   <  1:87841/236‑1 (MQ=255)
                                                                                      cATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGc                                                                   <  1:16152/236‑1 (MQ=14)
                                                                                                tGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAg                                                         >  2:271912/1‑236 (MQ=14)
                                                                                                 gCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGACAGGGAGTTTCAGCAGTGGAACGGCACAGCCTCGCTGAAgg                                                        >  2:595995/1‑236 (MQ=9)
                                                                                                 gCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTgg                                                               <  1:302012/229‑1 (MQ=14)
                                                                                                 gCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTg                                                            <  1:614055/232‑1 (MQ=14)
                                                                                                  ccAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGa                                                       >  2:656508/1‑236 (MQ=14)
                                                                                                  ccAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGa                                                       >  1:266312/1‑236 (MQ=14)
                                                                                                  ccAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGa                                                       >  1:363338/1‑236 (MQ=14)
                                                                                                  ccAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCCTGGTGAAGGa                                                       >  2:286913/1‑236 (MQ=14)
                                                                                                       aaaCAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAAcgcg                                                                                                                                                         <  1:723045/133‑1 (MQ=17)
                                                                                                       aaaCAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAAcgcg                                                                                                                                                         >  2:723045/1‑133 (MQ=17)
                                                                                                           tgTACCGTTATTGCACCGCCAGATATGCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGc                                              <  2:308147/235‑1 (MQ=14)
                                                                                                            gTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCa                                             >  1:195959/1‑236 (MQ=14)
                                                                                                            gTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGACATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACCTGGCTTCCGGCAACGCGTTATTTATTTTTGAACTCGGTCCGTTACCGGAACATGTTACCTGGCTATCGTCGCACGGGGAGTTTCTGAAGTGGAACGGCACAGCCTGGGTGATGGAGACGGAAgat                                             >  2:295196/1‑234 (MQ=11)
                                                                                                                  atatTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGa                                           >  1:14620/1‑232 (MQ=14)
                                                                                                                       gCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGtt                                  >  1:369064/1‑236 (MQ=14)
                                                                                                                       gCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGtt                                  <  1:71811/236‑1 (MQ=14)
                                                                                                                        cACCGCCAGTTTTGCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGGTCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTc                                 <  2:655105/236‑1 (MQ=9)
                                                                                                                         aCCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACt                                                                                                                                        <  2:214127/132‑1 (MQ=17)
                                                                                                                         aCCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACt                                                                                                                                        >  1:214127/1‑132 (MQ=17)
                                                                                                                          ccGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCg                               <  2:363338/236‑1 (MQ=14)
                                                                                                                           cGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgg                                                                                              >  2:764729/1‑172 (MQ=14)
                                                                                                                           cGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgg                                                                                              <  1:764729/172‑1 (MQ=14)
                                                                                                                             ccAGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCAAAGCCTGGGTGAAGGATACGGAAGCGGAAAAACGGTTCGGGAt                            >  2:102855/1‑236 (MQ=9)
                                                                                                                               aGATATTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATcc                          <  2:651918/236‑1 (MQ=14)
                                                                                                                                 atatTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTg                                                            <  2:18028/200‑1 (MQ=14)
                                                                                                                                 atatTCCGGCTGGCTTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTg                                                            >  1:18028/1‑200 (MQ=14)
                                                                                                                                          ctggctTCGTGGCTGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGgaa               >  1:534912/1‑236 (MQ=14)
                                                                                                                                               ttCGTGGCTGTTTGCAACATTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAAc          <  1:173081/236‑1 (MQ=14)
                                                                                                                                                tCGTGGGTGTTTGCAACAGTTATGAGTAAGCGTAGCATCTCGTTGAAGATCATCGGGGTAAAGCGGTTTAGGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCAGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAg                                               <  1:201453/198‑1 (MQ=11)
                                                                                                                                                       tGTTTGCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAaaca  >  1:738213/1‑236 (MQ=14)
                                                                                                                                                            |                                                                                                                                                                                                                                      
GTGAACAACCACGGACCATAAAAATTTATAATCTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGATGCATATATTCCGCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAAACA  >  NC_000913/1432427‑1432813

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