Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,356 T→G 11.8% S259A (TCC→GCC)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,3560T→G11.8% 108.5 ‑0.8 ‑1.4290S259A (TCC→GCC) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (20/14):  ref base (129/126):  total (149/141)
Fisher's exact test strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.61e-02

Only 100 of 296 total aligned reads displayed.

GATGATTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGG  >  NC_007779/3441256‑3441456
                                                                                                    |                                                                                                    
 atgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTc                                                                                                     >  1:8333112/1‑101 (MQ=255)
 atgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTc                                                                                                     >  1:12130960/1‑101 (MQ=255)
  tgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTcc                                                                                                    <  1:11198968/101‑1 (MQ=255)
     ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                 >  1:9973128/1‑101 (MQ=255)
     ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                 <  1:9451475/101‑1 (MQ=255)
      tattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcgc                                                                                                <  1:2737153/101‑1 (MQ=255)
         tGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCt                                                                                             >  1:9599800/1‑101 (MQ=255)
           gccGCAATTTCCTGGTAAAAGTTACCGCCATTACCGGCAACTGGCCGGCCCCTTCTTCCATGGAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTggg                                                                                          <  1:4405255/99‑1 (MQ=255)
           gTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTgg                                                                                           <  1:6607053/101‑1 (MQ=255)
            tCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTggg                                                                                          >  1:2476499/1‑101 (MQ=255)
             cGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTgggg                                                                                         >  1:4344783/1‑101 (MQ=255)
             cGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTgggg                                                                                         <  1:2566281/101‑1 (MQ=255)
               cAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAg                                                                                       >  1:11286503/1‑101 (MQ=255)
                aaTTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGc                                                                                      >  1:1872294/1‑101 (MQ=255)
                aaTTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGc                                                                                      >  1:11816852/1‑101 (MQ=255)
                 aTTTCATGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAATTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCg                                                                                     <  1:4031480/101‑1 (MQ=255)
                  tttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCgg                                                                                    <  1:7312497/101‑1 (MQ=255)
                     ccTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATa                                                                                 >  1:3491845/1‑101 (MQ=255)
                      cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                >  1:6911134/1‑101 (MQ=255)
                       tGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACg                                                                               >  1:9820763/1‑101 (MQ=255)
                        ggTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACgg                                                                              >  1:3129216/1‑101 (MQ=255)
                         gTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGt                                                                             <  1:494457/101‑1 (MQ=255)
                            aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                          >  1:716952/1‑101 (MQ=255)
                            aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                          >  1:11972456/1‑101 (MQ=255)
                            aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                          <  1:14266701/101‑1 (MQ=255)
                             aaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATg                                                                         <  1:4161209/101‑1 (MQ=255)
                             aaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATg                                                                         >  1:9316448/1‑101 (MQ=255)
                             aaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATg                                                                         >  1:12468346/1‑101 (MQ=255)
                                 tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                     >  1:12351320/1‑101 (MQ=255)
                                 tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                     >  1:12768123/1‑101 (MQ=255)
                                 tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                     >  1:11017455/1‑101 (MQ=255)
                                 tAACGCCAATATCGGCAACTCGGCGGTCACCTCTGCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                     <  1:8807045/101‑1 (MQ=255)
                                  aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                    >  1:10455123/1‑101 (MQ=255)
                                  aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                    <  1:7606978/101‑1 (MQ=255)
                                    cGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctct                                                                  >  1:5483074/1‑101 (MQ=255)
                                      ccAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTcc                                                                >  1:9523007/1‑101 (MQ=255)
                                       cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                               <  1:10142593/101‑1 (MQ=255)
                                       cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                               >  1:7855138/1‑101 (MQ=255)
                                        aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                              >  1:4226824/1‑101 (MQ=255)
                                        aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                              >  1:1267546/1‑101 (MQ=255)
                                        aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                              >  1:2225557/1‑101 (MQ=255)
                                           atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                           <  1:13347408/101‑1 (MQ=255)
                                           atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                           >  1:13476074/1‑101 (MQ=255)
                                             cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                         <  1:46575/101‑1 (MQ=255)
                                             cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                         >  1:3812263/1‑101 (MQ=255)
                                             cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCACTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                         >  1:2000289/1‑101 (MQ=255)
                                                 aaCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCta                                                     >  1:684274/1‑101 (MQ=255)
                                                 aaCGCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCta                                                     <  1:129661/101‑1 (MQ=255)
                                                   cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                   <  1:7537819/101‑1 (MQ=255)
                                                   cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                   >  1:6666469/1‑101 (MQ=255)
                                                      ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                <  1:5012446/101‑1 (MQ=255)
                                                      ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                <  1:3880729/101‑1 (MQ=255)
                                                      ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                >  1:7094417/1‑101 (MQ=255)
                                                       gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                               <  1:11570584/101‑1 (MQ=255)
                                                       gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                               <  1:8584691/101‑1 (MQ=255)
                                                        cggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCAc                                              >  1:6246397/1‑101 (MQ=255)
                                                         ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                             >  1:1764209/1‑101 (MQ=255)
                                                         ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                             >  1:10808927/1‑101 (MQ=255)
                                                           tCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaa                                           <  1:12843928/101‑1 (MQ=255)
                                                            cACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaaa                                          <  1:3636979/101‑1 (MQ=255)
                                                               ctctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAccc                                       <  1:8001414/101‑1 (MQ=255)
                                                                tctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcg                                      >  1:922457/1‑101 (MQ=255)
                                                                  ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                    >  1:6169851/1‑101 (MQ=255)
                                                                   tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                   <  1:7726155/101‑1 (MQ=255)
                                                                   tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                   <  1:10877885/101‑1 (MQ=255)
                                                                    ccatctccgaagaggTGGAGAAGATGGTGTGGGCGTCCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAg                                  >  1:6557451/9‑101 (MQ=255)
                                                                    ccATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAg                                  <  1:7266565/101‑1 (MQ=255)
                                                                     catctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAgg                                 >  1:13738428/8‑100 (MQ=255)
                                                                     catctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGt                                 <  1:164767/94‑1 (MQ=255)
                                                                     cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                 <  1:13679159/101‑1 (MQ=255)
                                                                     cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                 <  1:9124762/101‑1 (MQ=255)
                                                                      atctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTg                                <  1:1416718/95‑1 (MQ=255)
                                                                      aTCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTg                                <  1:9027061/101‑1 (MQ=255)
                                                                      aTCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTg                                >  1:6993477/1‑101 (MQ=255)
                                                                       tCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTgg                               >  1:9775658/1‑101 (MQ=255)
                                                                       tCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTgg                               <  1:10394644/101‑1 (MQ=255)
                                                                        ctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGa                              >  1:10765358/5‑101 (MQ=255)
                                                                          ccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATACACGAAACCCGCGAGTGGAtt                            <  1:2848079/99‑1 (MQ=255)
                                                                           agaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAttt                           >  1:13275002/1‑101 (MQ=255)
                                                                            gaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAttct                          >  1:5367111/1‑99 (MQ=255)
                                                                            gaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAttct                          >  1:3572369/1‑99 (MQ=255)
                                                                            gaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtttt                          >  1:8907280/1‑101 (MQ=255)
                                                                             aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                         >  1:4635824/1‑101 (MQ=255)
                                                                             aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                         <  1:11406928/101‑1 (MQ=255)
                                                                              agaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGc                        >  1:5242561/1‑101 (MQ=255)
                                                                              agaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGc                        <  1:5127505/101‑1 (MQ=255)
                                                                               gaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCg                       <  1:1219757/98‑1 (MQ=255)
                                                                                aagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGt                      <  1:333129/101‑1 (MQ=255)
                                                                                    ggAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAAct                  >  1:450494/1‑100 (MQ=255)
                                                                                    ggAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACa                  <  1:10668443/101‑1 (MQ=255)
                                                                                      aaaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGc                >  1:10003788/1‑101 (MQ=255)
                                                                                        aaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGccc              <  1:13975385/101‑1 (MQ=255)
                                                                                         aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg             >  1:5979394/1‑101 (MQ=255)
                                                                                            tGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTg          >  1:12028769/1‑101 (MQ=255)
                                                                                             ggTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGc         <  1:8548026/101‑1 (MQ=255)
                                                                                             ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc         <  1:8835039/101‑1 (MQ=255)
                                                                                             ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGc         >  1:5778246/1‑101 (MQ=255)
                                                                                             ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGATCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc         <  1:10294207/101‑1 (MQ=255)
                                                                                               tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg       >  1:8097868/1‑101 (MQ=255)
                                                                                                    tCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCgg  <  1:1579533/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GATGATTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGG  >  NC_007779/3441256‑3441456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.