Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,044 T→C 5.3% F488S (TTT→TCT)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,0440T→C 5.3% 24.0 ‑2.0 ‑1.1225F488S (TTT→TCT) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (8/4):  ref base (90/123):  total (98/127)
Fisher's exact test strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.87e-02

Only 100 of 244 total aligned reads displayed.

CACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGC  >  NC_007779/3441944‑3442143
                                                                                                     |                                                                                                   
cACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGtt                                                                                                     <  1:6453886/101‑1 (MQ=255)
  cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                                                                                   >  1:8909127/1‑101 (MQ=255)
     aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                                                                                <  1:4235961/101‑1 (MQ=255)
      gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                                                                                               <  1:12116639/101‑1 (MQ=255)
       cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTTGCTgc                                                                                              >  1:12675133/1‑101 (MQ=255)
        gcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcg                                                                                             <  1:3060790/101‑1 (MQ=255)
         ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCCCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                                            <  1:11992107/101‑1 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           >  1:6261907/1‑101 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           >  1:8278608/1‑101 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           <  1:8385815/101‑1 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           >  1:5195194/1‑101 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           >  1:5061392/1‑101 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           <  1:10712738/101‑1 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           >  1:9025261/1‑101 (MQ=255)
          cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                           >  1:3307960/1‑101 (MQ=255)
            ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                                                         >  1:2420846/1‑101 (MQ=255)
            ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGAt                                                                                         >  1:1358037/1‑101 (MQ=255)
             tttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATg                                                                                        >  1:1499119/1‑101 (MQ=255)
               tACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGct                                                                                      <  1:79210/101‑1 (MQ=255)
                acacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctc                                                                                     <  1:5687325/101‑1 (MQ=255)
                 cacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGGTGCGCGATGctct                                                                                    >  1:3744810/1‑101 (MQ=255)
                 cacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctct                                                                                    >  1:1325540/1‑101 (MQ=255)
                     ctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTa                                                                                <  1:5765864/101‑1 (MQ=255)
                       ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                                              >  1:13308615/1‑101 (MQ=255)
                       ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                                              <  1:13936531/101‑1 (MQ=255)
                       ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                                              >  1:9255151/1‑101 (MQ=255)
                       agggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                                              >  1:10058201/2‑101 (MQ=255)
                          gCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTaa                                                                           <  1:7123666/101‑1 (MQ=255)
                          gCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTaa                                                                           >  1:11800486/1‑101 (MQ=255)
                           ccGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAAc                                                                          >  1:8986606/1‑101 (MQ=255)
                            cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                                                         <  1:9736532/101‑1 (MQ=255)
                            cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACg                                                                         >  1:7682134/1‑101 (MQ=255)
                              cTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGcc                                                                       <  1:760251/101‑1 (MQ=255)
                              cTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGcc                                                                       <  1:163534/101‑1 (MQ=255)
                               tAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCa                                                                      >  1:3644447/1‑101 (MQ=255)
                                aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                     >  1:10644174/1‑101 (MQ=255)
                                aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                     <  1:6938045/101‑1 (MQ=255)
                                aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCCGCGCGATGCTCTGTTACGTAACGCCaa                                                                     >  1:4261282/1‑101 (MQ=255)
                                  cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                                                   <  1:4804565/101‑1 (MQ=255)
                                  cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                                                   <  1:5378619/101‑1 (MQ=255)
                                   tACCGATATTGCGCCGGGCTATGTCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                                  <  1:8721583/101‑1 (MQ=255)
                                   tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                                  >  1:10324378/1‑101 (MQ=255)
                                   tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                                  <  1:7395214/101‑1 (MQ=255)
                                   aaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCgaagg                                                                  >  1:4615764/2‑99 (MQ=255)
                                   aaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCgaagg                                                                  >  1:12275970/2‑99 (MQ=255)
                                    aCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAaga                                                                 <  1:765434/101‑1 (MQ=255)
                                     ccGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                                >  1:5298895/1‑101 (MQ=255)
                                       gacatcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCa                                                              >  1:12888273/7‑101 (MQ=255)
                                         tatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATc                                                            >  1:2825191/1‑101 (MQ=255)
                                          atcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCt                                                           <  1:5969454/98‑1 (MQ=255)
                                          atTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCt                                                           <  1:10447992/101‑1 (MQ=255)
                                           ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg                                                          <  1:13141453/101‑1 (MQ=255)
                                            tGCGCCGGGCTATGACCACTTCACGTCGGCGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg                                                         >  1:8142239/1‑101 (MQ=255)
                                              cgcCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGt                                                       <  1:11730317/101‑1 (MQ=255)
                                               gcCGGGCTATGACCACTTCACGTCGGGGATTGGAGCGGCGTTGTTTGGCTGGTTTGGCTGCGCGACGCTCTGTTTCGTAACGCCAAAAGAGCATCTgggcc                                                      >  1:1385426/1‑99 (MQ=255)
                                                 cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                                                    >  1:8530109/1‑101 (MQ=255)
                                                 cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                                                    >  1:7014789/1‑101 (MQ=255)
                                                  gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                                   <  1:5817937/101‑1 (MQ=255)
                                                   ggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGcc                                                  <  1:4103267/101‑1 (MQ=255)
                                                        tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                             <  1:1049386/101‑1 (MQ=255)
                                                         gACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa                                            >  1:2204692/1‑101 (MQ=255)
                                                          aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaaa                                           <  1:8186025/101‑1 (MQ=255)
                                                            cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGGGATGCTCTGTTACGGAACGCCAAAAGAGCATAGGGGGTTGGCCCATAaaga                                         >  1:3439701/1‑101 (MQ=255)
                                                            cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                                         >  1:11290704/1‑101 (MQ=255)
                                                               ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGTCCAATAAagaaga                                      <  1:1947835/101‑1 (MQ=255)
                                                               ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                                      <  1:8591628/101‑1 (MQ=255)
                                                                tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                     <  1:13559411/101‑1 (MQ=255)
                                                                 cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                                    <  1:13396758/101‑1 (MQ=255)
                                                                 cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                                    >  1:7726340/1‑101 (MQ=255)
                                                                 cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGGGAGGCTCTGTTACGGAAAGCCAAAAGAGCATCCGGGTCTGCCCAATGAAGAAGATg                                    >  1:319961/1‑101 (MQ=255)
                                                                  aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt                                   <  1:5627652/101‑1 (MQ=255)
                                                                   cGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                                  >  1:6715266/1‑101 (MQ=255)
                                                                     tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa                                <  1:4894091/101‑1 (MQ=255)
                                                                      cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGTCCAATAAAGAAGATGTTAAg                               <  1:1785350/101‑1 (MQ=255)
                                                                      cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg                               >  1:802055/1‑101 (MQ=255)
                                                                       ggggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc                              <  1:9592434/101‑1 (MQ=255)
                                                                        gggCTTGGTGCGCCGATAATTGGCTGGTTTGGCTGCGCGAAGCGCTGTTGCGTAACGCCAAAAGAGGAACTGGGTCTGACCAATAAAGAAGATGTTAAGCa                             >  1:9032303/1‑101 (MQ=255)
                                                                        gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa                             >  1:12786883/1‑101 (MQ=255)
                                                                        gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa                             <  1:2871870/101‑1 (MQ=255)
                                                                          gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg                           <  1:5165789/101‑1 (MQ=255)
                                                                          gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg                           <  1:497769/101‑1 (MQ=255)
                                                                           aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAggg                          <  1:13659244/101‑1 (MQ=255)
                                                                            ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg                         <  1:10026019/101‑1 (MQ=255)
                                                                              ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt                       >  1:1719339/1‑101 (MQ=255)
                                                                               gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt                      <  1:8646174/101‑1 (MQ=255)
                                                                               gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt                      <  1:12907949/101‑1 (MQ=255)
                                                                                 gcggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTAt                    >  1:14254548/1‑101 (MQ=255)
                                                                                  cggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATc                   <  1:2556964/101‑1 (MQ=255)
                                                                                   ggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa                  >  1:10929483/1‑101 (MQ=255)
                                                                                     cgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCAcc                >  1:5568237/1‑101 (MQ=255)
                                                                                      gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt               <  1:8877445/101‑1 (MQ=255)
                                                                                      gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt               <  1:12148188/101‑1 (MQ=255)
                                                                                      gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt               <  1:9730579/101‑1 (MQ=255)
                                                                                         gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtata            >  1:9610361/1‑101 (MQ=255)
                                                                                         gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtata            <  1:8359349/101‑1 (MQ=255)
                                                                                            tGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGa         >  1:7224123/1‑101 (MQ=255)
                                                                                             ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt        >  1:9388689/1‑101 (MQ=255)
                                                                                                  ggTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtgct    >  1:12551948/1‑101 (MQ=255)
                                                                                                   gTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtgctg   >  1:12198544/1‑101 (MQ=255)
                                                                                                    tttGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtgctgc  >  1:6556076/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
CACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGC  >  NC_007779/3441944‑3442143

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