Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 366,409 C→T intergenic (‑104/+19) lacZ ← / ← lacI beta‑galactosidase/DNA‑binding transcriptional repressor LacI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,4090CT100.0% 128.8 / NA 40intergenic (‑104/+19)lacZ/lacIbeta‑galactosidase/DNA‑binding transcriptional repressor LacI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/22);  total (18/22)

GGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGT  >  NC_000913/366163‑366652
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gggCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACAt                                                                                                                                                                                                                                                    >  1:371515/1‑251 (MQ=255)
    gATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAAt                                                                                                                                                                                                                                                >  2:157862/1‑251 (MQ=255)
          tGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGGGTAAAGCCTGGGGGGCCTAATGAGTGAGCTAACTTACATTAAttgcgtt                                                                                                                                                                                                                                          >  1:104631/1‑251 (MQ=255)
                ccgctTCGCTATTACGCCAGCTGGCGAAAGGGTGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTc                                                                                                                                                                                                                                    <  2:117468/248‑1 (MQ=255)
                 ctctTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTc                                                                                                                                                                                                                                    >  1:9831/1‑250 (MQ=255)
                  tctTCGCAATTACGCCAGCCGGCGAAAGGAGGATGCGCTGCACGGCGATTCAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAACTTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGATGGAGCTACCTTACATTATTTGGGTGGCGCTCAc                                                                                                                                                                                                                                  >  2:333143/1‑251 (MQ=255)
                             aCGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCttt                                                                                                                                                                                                                        <  2:381844/250‑1 (MQ=255)
                                          aaaGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAAcc                                                                                                                                                                                                          >  2:461924/1‑251 (MQ=255)
                                                         gcAAGGCGATTAAGTGGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCGGTCGTGCCAGCTg                                                                                                                                                                                           <  1:69762/250‑1 (MQ=255)
                                                           aaGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCAt                                                                                                                                                                                        >  1:157481/1‑251 (MQ=255)
                                                               cGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAAt                                                                                                                                                                                    <  1:429747/251‑1 (MQ=255)
                                                                gATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGAGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATg                                                                                                                                                                                   >  1:487507/1‑251 (MQ=255)
                                                                 aTTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGTCGTTGTAAAACGACGGCCAGTGAATCCGTCATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTACTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGAGCCAGCTGCATTAATgc                                                                                                                                                                                  >  2:359213/1‑250 (MQ=255)
                                                                   tAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGCCGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAt                                                                                                                                                                                >  1:140668/1‑251 (MQ=255)
                                                                        tGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGc                                                                                                                                                                            <  1:500948/250‑1 (MQ=255)
                                                                             aaCGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAAcgc                                                                                                                                                                      >  1:104978/1‑251 (MQ=255)
                                                                               cGCCAGGGTTTTCCCAGTCACGACGTTGTAAACCGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCTGCTGCATTAATGAATCGGCCAAcgcgc                                                                                                                                                                    >  1:175199/1‑251 (MQ=255)
                                                                                     ggTTTTCCCAGGCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGGTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGgag                                                                                                                                                              <  2:465551/251‑1 (MQ=255)
                                                                                                 cACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTTCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGc                                                                                                                                                   <  1:304705/250‑1 (MQ=255)
                                                                                                     acgTTGAAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGAATAAAGTGTAAAGCCTGGGGAGCCTCATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCCTTAATGAATCGGCCAACGCGCGGGGAGAAGCGGTTTGCGTAtt                                                                                                                                              >  1:450243/1‑251 (MQ=255)
                                                                                                       gTTGTAGAACGACGGCCAGAGAATCCGGAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAATATACGAGACGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGtgtttg                                                                                                                                            >  2:410746/1‑246 (MQ=255)
                                                                                                         agTAAAACGCCGGCCAGTGAAGCCGTAAGCATGGGCATAGCTTTTTCCTGTGTGAAATTGTGATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCTCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgc                                                                                                                                          <  1:266716/250‑1 (MQ=255)
                                                                                                          gTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgcg                                                                                                                                         <  1:279070/251‑1 (MQ=255)
                                                                                                                           gAATCCTTAATCATGGTCAGAGCTTGTTCCTGTGTGAAATTGTTATCCGCTCACAAGTCCCCACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGGCGGGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCTCCAGGGTGGTttttct                                                                                                                        <  2:488517/250‑1 (MQ=255)
                                                                                                                               ccGTAATCATGGGCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACTCAACATACTAGCCGGAAGCATAAAGTGTAAAGCGTGGGGTGCCTAATGAGTTAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCTTTCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                                                                                   <  1:54025/251‑1 (MQ=255)
                                                                                                                                  tAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGAGGGGAGAGGCGGTTTTCGTATTGGGCGCCCGGGTGGTGTTTTTTTTTAcc                                                                                                                 >  1:472677/1‑250 (MQ=255)
                                                                                                                                   aaTCATGGTCATAGCTGTTTCCTGTGAGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAg                                                                                                               <  1:279381/251‑1 (MQ=255)
                                                                                                                                          gTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCCCAATTCCACACAACATACGAGGCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                        <  1:333143/251‑1 (MQ=255)
                                                                                                                                            cATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                       <  1:461924/250‑1 (MQ=255)
                                                                                                                                                   gTTTCCTGAGTGAAATTGTTATCCTCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                <  1:421389/250‑1 (MQ=255)
                                                                                                                                                          gtgtgAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                        <  1:442865/251‑1 (MQ=255)
                                                                                                                                                                   tGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGTGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACATCTGATTGCCCTTCAcc                                                                                <  2:307551/250‑1 (MQ=255)
                                                                                                                                                                                            aaatcagcgCCGGAAGCATAAAGTGGAAAGCCTGGGGTGCCTAATTAGTGAGCTAACTTACATTAATTGCGTTGCGCACACTGCCCGCTTTCCAGTCGGGAAACGTGTCGTGCCAGCTGCTTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCCGGGTGGTTTATCTGTTCACCAGTGGTACGGGTAACAGCTGATTTCCCTTCACCGCCTGGCCCTCAGAGGTTTGCAGCaa                                                      <  2:223361/243‑1 (MQ=255)
                                                                                                                                                                                                     ccGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCAc                                             <  1:457773/251‑1 (MQ=255)
                                                                                                                                                                                                               aaaGTGTAAAGCCTGGGGAGCCTAATGAGTGAGCTAACTTAAATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGGGGCCGGGGGGGTTTTTTTTTTAACCGTGAGACGGGGAACAGCTGATTGGCCCTCACCACCCGGCCCCGAGGGAGTTGCAGCaaagggtcccccccgggttgc                                   >  1:400774/1‑232 (MQ=255)
                                                                                                                                                                                                                              gggTGCCTAATGAGTGAGCTAACTTACATTAATTGCGGTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAAGTAATCGGCCAACGCGCGGGGAGAGGCGTTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTGCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGTCCCTGAGAGAGGTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                   <  2:467318/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTAGGATTTGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGATAAGGGGCACAGCTGTGTGCCCTTCACCGCCTCGCCCTGTGAGAATTGCAGCCAACGGGCCCCGCTGGTTTGCCCCAGCAGGCGTAAATACTgat            >  1:16123/1‑249 (MQ=255)
                                                                                                                                                                                                                                      aaTGAGTGAGCTAACGTACATTAATGGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCGTGCGTGCCAGCTGCTTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCTGGGCCTTAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt           <  2:40208/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg          <  1:293741/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGCGGCAGGGGGGGTTTTCTTTTTACCAGTGAGGCGGGGAACAGATGATTTCCCTTCCCCGCCTGGCCCTCAGGGGGATGCAGCAAACGGGCCACCCCGGTTTGGCCCCACAGGCggcaatccctgtttattgtggg  >  1:385129/1‑230 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGT  >  NC_000913/366163‑366652

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