Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,356 T→G 12.4% 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→G12.4% 100.9 ‑1.6 ‑0.3264S259A (TCC→GCC) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (13/20):  ref base (120/111):  total (133/131)
Fisher's exact test strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.36e-01

Only 100 of 272 total aligned reads displayed.

GATGATTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGG  >  NC_007779/3441256‑3441456
                                                                                                    |                                                                                                    
gatgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGt                                                                                                      >  1:6448858/1‑101 (MQ=255)
 atgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTc                                                                                                     >  1:12560880/1‑101 (MQ=255)
 atgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTc                                                                                                     >  1:13799180/1‑101 (MQ=255)
 atgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTc                                                                                                     >  1:9093593/1‑101 (MQ=255)
  tgatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTcc                                                                                                    >  1:15010707/1‑101 (MQ=255)
   gatTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCa                                                                                                   <  1:14651841/101‑1 (MQ=255)
   gatTAGTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCa                                                                                                   <  1:432207/101‑1 (MQ=255)
     ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                 <  1:462335/101‑1 (MQ=255)
     ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                 <  1:14205474/101‑1 (MQ=255)
     ttattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcg                                                                                                 >  1:158234/1‑101 (MQ=255)
      tattGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCAcgc                                                                                                >  1:1065693/1‑101 (MQ=255)
         tGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCt                                                                                             >  1:8556367/1‑101 (MQ=255)
          ggTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTg                                                                                            <  1:606361/101‑1 (MQ=255)
            tCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTggg                                                                                          >  1:658858/1‑101 (MQ=255)
              gCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGa                                                                                        >  1:7407323/1‑101 (MQ=255)
                aaTTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGc                                                                                      >  1:1631621/1‑101 (MQ=255)
                   ttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGa                                                                                   >  1:3677357/1‑101 (MQ=255)
                   ttCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGa                                                                                   >  1:8093828/1‑101 (MQ=255)
                    tCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGAt                                                                                  <  1:12216897/101‑1 (MQ=255)
                    tCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGAt                                                                                  <  1:3837291/101‑1 (MQ=255)
                    tCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGAt                                                                                  <  1:9125040/101‑1 (MQ=255)
                      cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                <  1:6430461/101‑1 (MQ=255)
                      cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                <  1:7015799/101‑1 (MQ=255)
                      cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                >  1:4768258/1‑101 (MQ=255)
                      cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                <  1:9569650/101‑1 (MQ=255)
                      cTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATAc                                                                                <  1:10479161/101‑1 (MQ=255)
                           aaaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGa                                                                           >  1:5412990/1‑101 (MQ=255)
                            aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                          <  1:8456681/101‑1 (MQ=255)
                            aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                          >  1:14462750/1‑101 (MQ=255)
                            aaaGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGAt                                                                          >  1:3369951/1‑101 (MQ=255)
                               gTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGa                                                                       <  1:5111644/101‑1 (MQ=255)
                               gTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGa                                                                       <  1:8910/101‑1 (MQ=255)
                                ttAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAt                                                                      <  1:1325937/101‑1 (MQ=255)
                                 tAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtc                                                                     >  1:2571231/1‑101 (MQ=255)
                                  aaCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtct                                                                    >  1:11188752/1‑101 (MQ=255)
                                   aCGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctc                                                                   >  1:148349/1‑101 (MQ=255)
                                     gCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGAtctctc                                                                 >  1:2850577/1‑101 (MQ=255)
                                        aaTATCGGCAACTCGGCGGTCGCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                              >  1:10854029/1‑101 (MQ=255)
                                        aaTATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAc                                                              >  1:5891866/1‑101 (MQ=255)
                                         atatCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAcc                                                             <  1:3497712/101‑1 (MQ=255)
                                         atatCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCAcc                                                             <  1:1088106/101‑1 (MQ=255)
                                            tCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGt                                                          >  1:1981544/1‑101 (MQ=255)
                                            tCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGt                                                          <  1:1717641/101‑1 (MQ=255)
                                             cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                         >  1:10984352/1‑101 (MQ=255)
                                             cGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTc                                                         >  1:14858600/1‑101 (MQ=255)
                                               gCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGc                                                       >  1:8752032/1‑101 (MQ=255)
                                               gCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGc                                                       >  1:2888884/1‑101 (MQ=255)
                                                cAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCt                                                      >  1:12131150/1‑101 (MQ=255)
                                                 aaCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCta                                                     >  1:7574659/1‑101 (MQ=255)
                                                 aaCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCta                                                     >  1:2447416/1‑101 (MQ=255)
                                                  aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                    <  1:11432361/101‑1 (MQ=255)
                                                    tCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTatat                                                  >  1:12382036/1‑101 (MQ=255)
                                                       gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                               >  1:2294647/1‑101 (MQ=255)
                                                       gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                               <  1:6705030/101‑1 (MQ=255)
                                                       gcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCa                                               <  1:7820595/101‑1 (MQ=255)
                                                         ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                             <  1:10593037/101‑1 (MQ=255)
                                                           tCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaa                                           >  1:5009494/1‑101 (MQ=255)
                                                              ccTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAcc                                        >  1:9447939/1‑101 (MQ=255)
                                                                tctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcg                                      >  1:13224659/1‑101 (MQ=255)
                                                                  ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                    <  1:4481343/101‑1 (MQ=255)
                                                                  ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                    >  1:3947341/1‑101 (MQ=255)
                                                                   tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                   <  1:6993111/101‑1 (MQ=255)
                                                                    ccATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAg                                  >  1:1196915/1‑101 (MQ=255)
                                                                     cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                 >  1:7080633/1‑101 (MQ=255)
                                                                     cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                 >  1:2167253/1‑101 (MQ=255)
                                                                     cATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGt                                 <  1:5808465/101‑1 (MQ=255)
                                                                      aTCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTg                                <  1:13645554/101‑1 (MQ=255)
                                                                      aTCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGCCGCTATATTCACGAAACCCGCGAggg                                >  1:10115632/1‑99 (MQ=255)
                                                                        cGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGa                              <  1:1293203/101‑1 (MQ=255)
                                                                         tccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAt                             <  1:980838/98‑1 (MQ=255)
                                                                         gaagaagaaGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAt                             >  1:7699005/1‑101 (MQ=255)
                                                                          ccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAtt                            >  1:8336660/3‑101 (MQ=255)
                                                                           agaagacgTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAttt                           >  1:13407451/1‑101 (MQ=255)
                                                                            gaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAttct                          >  1:7241428/1‑99 (MQ=255)
                                                                              agaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGc                        <  1:15139730/101‑1 (MQ=255)
                                                                               gaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCg                       >  1:12672966/1‑101 (MQ=255)
                                                                                aagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGt                      <  1:11102296/101‑1 (MQ=255)
                                                                                aagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGt                      <  1:13508217/101‑1 (MQ=255)
                                                                                 ggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTa                     >  1:8235383/2‑101 (MQ=255)
                                                                                 agTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTa                     <  1:11544296/101‑1 (MQ=255)
                                                                                 agTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTa                     <  1:13277462/101‑1 (MQ=255)
                                                                                   tGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAAc                   >  1:12025035/1‑101 (MQ=255)
                                                                                    ggAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACa                  <  1:6648285/101‑1 (MQ=255)
                                                                                     gagaaGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAActc                 <  1:5300940/98‑3 (MQ=255)
                                                                                       aaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGcc               <  1:10575175/101‑1 (MQ=255)
                                                                                        aaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGccc              <  1:9688304/101‑1 (MQ=255)
                                                                                         agatGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCg             <  1:13090043/98‑1 (MQ=255)
                                                                                         aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg             <  1:2096677/101‑1 (MQ=255)
                                                                                          gCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCgg            <  1:6760159/101‑1 (MQ=255)
                                                                                          gCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTAACGAAACCCGCGAGGGGCTTTTGCGTAACAGCCCgg            >  1:9563783/1‑101 (MQ=255)
                                                                                             ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc         >  1:12856641/1‑101 (MQ=255)
                                                                                             ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc         <  1:5996975/101‑1 (MQ=255)
                                                                                              gTATAGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCTGTGcc        <  1:15086925/101‑1 (MQ=255)
                                                                                               tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg       <  1:2997684/101‑1 (MQ=255)
                                                                                                aTGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGa      >  1:11783796/1‑101 (MQ=255)
                                                                                                aTGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGa      >  1:11677590/1‑101 (MQ=255)
                                                                                                aTGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGa      >  1:13998493/1‑101 (MQ=255)
                                                                                                 tGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGAt     <  1:3108688/101‑1 (MQ=255)
                                                                                                  ggTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATc    <  1:10253501/101‑1 (MQ=255)
                                                                                                    tCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCgg  >  1:9396525/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GATGATTATTGGTCGCAATTTCCTGGTAAAAGTTAACGCCAATATCGGCAACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGG  >  NC_007779/3441256‑3441456

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