Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,358 C→G 11.8% S259S (TCC→TCG thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,3580C→G11.8% 98.1 ‑1.0 ‑1.2288S259S (TCC→TCGthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (20/14):  ref base (126/128):  total (146/142)
Fisher's exact test strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.93e-02

Only 100 of 293 total aligned reads displayed.

TGATTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTA  >  NC_007779/3441258‑3441458
                                                                                                    |                                                                                                    
tgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTcc                                                                                                      <  1:1186817/101‑1 (MQ=255)
tgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTcc                                                                                                      <  1:11198968/101‑1 (MQ=255)
 gatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCa                                                                                                     >  1:356302/1‑101 (MQ=255)
  attattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAc                                                                                                    >  1:3886809/1‑101 (MQ=255)
   ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                   >  1:9973128/1‑101 (MQ=255)
   ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                   <  1:7116203/101‑1 (MQ=255)
   ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                   <  1:9057347/101‑1 (MQ=255)
    tattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGGAGAAGTGGAAAAGCTGGTATGGTCCAcgc                                                                                                  <  1:336707/101‑1 (MQ=255)
     attGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcgcg                                                                                                 <  1:4324274/101‑1 (MQ=255)
        ggTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTg                                                                                              >  1:12807015/1‑101 (MQ=255)
         gTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTgg                                                                                             <  1:6607053/101‑1 (MQ=255)
           cGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTgggg                                                                                           >  1:12164271/1‑101 (MQ=255)
             cAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAg                                                                                         >  1:11286503/1‑101 (MQ=255)
              aaTTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGc                                                                                        <  1:2130133/101‑1 (MQ=255)
               aTTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCg                                                                                       <  1:833712/101‑1 (MQ=255)
                tttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCgg                                                                                      <  1:7312497/101‑1 (MQ=255)
                 ttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGa                                                                                     >  1:4191742/1‑101 (MQ=255)
                 ttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGa                                                                                     <  1:469535/101‑1 (MQ=255)
                 ttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGa                                                                                     >  1:8732792/1‑101 (MQ=255)
                    cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                  <  1:4584514/101‑1 (MQ=255)
                     tGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACg                                                                                 >  1:9820763/1‑101 (MQ=255)
                      ggTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATCCgg                                                                                >  1:8577065/1‑101 (MQ=255)
                      ggTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACgg                                                                                >  1:9650010/1‑101 (MQ=255)
                       gTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGt                                                                               >  1:12693232/1‑101 (MQ=255)
                       gTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGt                                                                               <  1:494457/101‑1 (MQ=255)
                         aaaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGa                                                                             <  1:6956263/101‑1 (MQ=255)
                          aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                            >  1:11972456/1‑101 (MQ=255)
                          aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                            <  1:7026674/101‑1 (MQ=255)
                          aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                            >  1:3114969/1‑101 (MQ=255)
                          aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                            >  1:9192328/1‑101 (MQ=255)
                            aGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGCGCGGATACGGTGATgg                                                                          <  1:9591147/101‑1 (MQ=255)
                              ttAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAt                                                                        >  1:2730682/1‑101 (MQ=255)
                              ttAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAt                                                                        <  1:8821340/101‑1 (MQ=255)
                               tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                       <  1:3148032/101‑1 (MQ=255)
                               tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                       >  1:11017455/1‑101 (MQ=255)
                                aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                      <  1:7606978/101‑1 (MQ=255)
                                aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                      >  1:10455123/1‑101 (MQ=255)
                                  cGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctct                                                                    >  1:5483074/1‑101 (MQ=255)
                                    ccAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTcc                                                                  >  1:4482548/1‑101 (MQ=255)
                                     cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                                 <  1:10142593/101‑1 (MQ=255)
                                     cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                                 >  1:9433634/1‑101 (MQ=255)
                                     cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                                 >  1:9241194/1‑101 (MQ=255)
                                     cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                                 >  1:7855138/1‑101 (MQ=255)
                                     cAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCa                                                                 <  1:6901879/101‑1 (MQ=255)
                                      aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                                <  1:12459816/101‑1 (MQ=255)
                                       atatCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAcc                                                               <  1:1419811/101‑1 (MQ=255)
                                         atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                             >  1:13476074/1‑101 (MQ=255)
                                         atCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCgg                                                             <  1:13347408/101‑1 (MQ=255)
                                            ggCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCg                                                          <  1:899729/101‑1 (MQ=255)
                                             gCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGGTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGc                                                         >  1:12101124/1‑101 (MQ=255)
                                              cAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCt                                                        <  1:2803231/101‑1 (MQ=255)
                                               aaCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCta                                                       >  1:684274/1‑101 (MQ=255)
                                                 cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                     >  1:6666469/1‑101 (MQ=255)
                                                 cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                     <  1:4638899/101‑1 (MQ=255)
                                                 cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                     <  1:7537819/101‑1 (MQ=255)
                                                 cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                     <  1:1685407/101‑1 (MQ=255)
                                                    ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                  <  1:5012446/101‑1 (MQ=255)
                                                     gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                                 >  1:11450826/1‑101 (MQ=255)
                                                     gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                                 <  1:8584691/101‑1 (MQ=255)
                                                      cggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCAc                                                >  1:6246397/1‑101 (MQ=255)
                                                       ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                               >  1:10808927/1‑101 (MQ=255)
                                                        gTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtagtgctgaag                                              <  1:4140997/101‑10 (MQ=255)
                                                            ccTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAcc                                          <  1:8897901/101‑1 (MQ=255)
                                                             ctctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAccc                                         <  1:8001414/101‑1 (MQ=255)
                                                                ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                      >  1:6169851/1‑101 (MQ=255)
                                                                 tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                     <  1:10877885/101‑1 (MQ=255)
                                                                 tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                     <  1:7726155/101‑1 (MQ=255)
                                                                 tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTACACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                     <  1:1811766/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)
                                                                   cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                   <  1:1528781/101‑1 (MQ=255)
                                                                   cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                   <  1:13679159/101‑1 (MQ=255)
                                                                     tCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTgg                                 <  1:10394644/101‑1 (MQ=255)
                                                                      ctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGa                                >  1:10765358/5‑101 (MQ=255)
                                                                        aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                              >  1:12882687/1‑101 (MQ=255)
                                                                        aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                              >  1:3586474/1‑101 (MQ=255)
                                                                        aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                              >  1:3812737/1‑101 (MQ=255)
                                                                        aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                              >  1:13769789/1‑101 (MQ=255)
                                                                        aagaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAtt                              <  1:11754605/101‑1 (MQ=255)
                                                                         agaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAttt                             <  1:950784/101‑1 (MQ=255)
                                                                          gaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAttct                            >  1:5367111/1‑99 (MQ=255)
                                                                            agaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGc                          >  1:161751/1‑101 (MQ=255)
                                                                            agaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGc                          >  1:5242561/1‑101 (MQ=255)
                                                                            agaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGc                          >  1:12550274/1‑101 (MQ=255)
                                                                            agaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGc                          <  1:5127505/101‑1 (MQ=255)
                                                                             gaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCg                         >  1:2434578/4‑101 (MQ=255)
                                                                               ggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTa                       >  1:1998057/2‑101 (MQ=255)
                                                                               ggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTa                       <  1:4356787/100‑1 (MQ=255)
                                                                                  ggAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACa                    <  1:3766939/101‑1 (MQ=255)
                                                                                  ggAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACa                    <  1:10668443/101‑1 (MQ=255)
                                                                                    aaaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGc                  >  1:10003788/1‑101 (MQ=255)
                                                                                       agatGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCg               <  1:13064131/98‑1 (MQ=255)
                                                                                       aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg               >  1:5979394/1‑101 (MQ=255)
                                                                                           ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGc           >  1:5778246/1‑101 (MQ=255)
                                                                                           ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGATCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc           <  1:10294207/101‑1 (MQ=255)
                                                                                             tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg         >  1:8097868/1‑101 (MQ=255)
                                                                                             tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg         >  1:1408884/1‑101 (MQ=255)
                                                                                              aTGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGa        >  1:12334138/1‑101 (MQ=255)
                                                                                                  tCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCgg    <  1:2528705/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGATTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTA  >  NC_007779/3441258‑3441458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ 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.