Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,731 T→G 26.6% N1236H (AAC→CAC)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,7310T→G26.6% 273.5 ‑0.5 ‑0.0301N1236H (AAC→CAC) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (45/35):  ref base (115/106):  total (160/141)
Fisher's exact test strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.33e-01

Only 100 of 313 total aligned reads displayed.

CGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCG  >  NC_007779/3451631‑3451831
                                                                                                    |                                                                                                    
  cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                                                                    >  1:10191905/1‑101 (MQ=255)
  cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                                                                    <  1:563206/101‑1 (MQ=255)
  cTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCa                                                                                                    <  1:11821085/101‑1 (MQ=255)
   tGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAg                                                                                                   <  1:7368846/101‑1 (MQ=255)
   tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                                                   <  1:485599/101‑1 (MQ=255)
   tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                                                   <  1:3661027/101‑1 (MQ=255)
   tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                                                   <  1:14594660/101‑1 (MQ=255)
    gACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGt                                                                                                  >  1:14137705/1‑101 (MQ=255)
     accaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGtt                                                                                                 >  1:1868466/1‑101 (MQ=255)
     accaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGtt                                                                                                 >  1:8586959/1‑101 (MQ=255)
      ccaccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGctt                                                                                                <  1:6443706/101‑1 (MQ=255)
       caccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTc                                                                                               <  1:5668249/101‑1 (MQ=255)
       caccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTc                                                                                               <  1:6381281/101‑1 (MQ=255)
        accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                                                                              >  1:8320416/1‑101 (MQ=255)
           gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                                           <  1:13156732/101‑1 (MQ=255)
           gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                                           >  1:8392202/1‑101 (MQ=255)
           gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                                           >  1:6329484/1‑101 (MQ=255)
            aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                                                                          <  1:5052691/101‑1 (MQ=255)
             aCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatg                                                                                         <  1:14824831/101‑1 (MQ=255)
              cTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgt                                                                                        <  1:1222986/101‑1 (MQ=255)
              cTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgt                                                                                        >  1:6201604/1‑101 (MQ=255)
               tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGcatgta                                                                                       >  1:13052361/1‑101 (MQ=255)
               tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgta                                                                                       >  1:9433069/1‑101 (MQ=255)
               tgggCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgta                                                                                       >  1:9616076/4‑101 (MQ=255)
                  gCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacat                                                                                    <  1:13649281/101‑1 (MQ=255)
                   ccTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatg                                                                                   >  1:5092396/1‑101 (MQ=255)
                    cTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgt                                                                                  >  1:1099812/1‑101 (MQ=255)
                    cTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatgt                                                                                  <  1:3264961/101‑1 (MQ=255)
                    cTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatgt                                                                                  <  1:14305986/101‑1 (MQ=255)
                      tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                                                >  1:13573798/1‑101 (MQ=255)
                      tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                                                <  1:15080073/101‑1 (MQ=255)
                      tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                                                <  1:7031146/101‑1 (MQ=255)
                      tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                                                >  1:2061296/1‑101 (MQ=255)
                         caccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccc                                                                             >  1:7676872/1‑100 (MQ=255)
                          accCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccca                                                                            <  1:884714/101‑3 (MQ=255)
                           cccAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAc                                                                           >  1:10774932/1‑101 (MQ=255)
                           cccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAAc                                                                           <  1:10837943/101‑1 (MQ=255)
                            ccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACg                                                                          >  1:1285489/1‑101 (MQ=255)
                              aGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccga                                                                        >  1:1265599/1‑98 (MQ=255)
                              aGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGt                                                                        <  1:12095724/101‑1 (MQ=255)
                               gCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATATTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGtt                                                                       <  1:4360861/101‑1 (MQ=255)
                                gggCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaa                                                                      <  1:15146871/100‑6 (MQ=255)
                                cGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaa                                                                      <  1:2503915/101‑6 (MQ=255)
                                cGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaa                                                                      <  1:1706622/101‑6 (MQ=255)
                                cGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTa                                                                      <  1:14703304/101‑1 (MQ=255)
                                 ggCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaac                                                                     <  1:13196794/101‑7 (MQ=255)
                                 ggCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAc                                                                     >  1:4037800/1‑101 (MQ=255)
                                  gctgctGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAcc                                                                    <  1:2443572/101‑1 (MQ=255)
                                    tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg                                                                  <  1:9770445/101‑1 (MQ=255)
                                    tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg                                                                  >  1:1578919/1‑101 (MQ=255)
                                       tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg                                                               >  1:10421470/1‑101 (MQ=255)
                                       tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg                                                               >  1:6592507/1‑101 (MQ=255)
                                          gTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTc                                                            <  1:982089/101‑1 (MQ=255)
                                          gTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTc                                                            >  1:1192616/1‑101 (MQ=255)
                                           tAACCATGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                                                           >  1:8727979/1‑101 (MQ=255)
                                            aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                                                          <  1:12786146/101‑1 (MQ=255)
                                            aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                                                          >  1:809216/1‑101 (MQ=255)
                                             aCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGaa                                                         <  1:2726975/101‑1 (MQ=255)
                                              ccAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAAc                                                        <  1:4317000/101‑1 (MQ=255)
                                               cAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACt                                                       >  1:341333/1‑101 (MQ=255)
                                                   cTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTc                                                   >  1:6084498/1‑101 (MQ=255)
                                                    tGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCa                                                  >  1:1610701/1‑101 (MQ=255)
                                                      tAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCAcc                                                <  1:7160175/101‑1 (MQ=255)
                                                       aaGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCa                                               >  1:13699774/1‑101 (MQ=255)
                                                       aaGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCa                                               <  1:5951634/101‑1 (MQ=255)
                                                         gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                                             >  1:3884852/1‑101 (MQ=255)
                                                         gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                                             >  1:2205089/1‑101 (MQ=255)
                                                           aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                           >  1:11612248/1‑101 (MQ=255)
                                                           aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                           <  1:8480102/101‑1 (MQ=255)
                                                             cGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGa                                         >  1:7247919/1‑101 (MQ=255)
                                                              ggTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGAc                                        >  1:2941227/1‑101 (MQ=255)
                                                               gTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGAcc                                       <  1:12984496/101‑1 (MQ=255)
                                                                tGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccg                                      >  1:9937658/1‑100 (MQ=255)
                                                                tGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccg                                      >  1:4212530/1‑100 (MQ=255)
                                                                 ggAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCAt                                     <  1:7426083/101‑1 (MQ=255)
                                                                 ggAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCAt                                     <  1:13936059/101‑1 (MQ=255)
                                                                  gAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCGTc                                    <  1:10740883/101‑1 (MQ=255)
                                                                  gAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATc                                    >  1:1897187/1‑101 (MQ=255)
                                                                   aaCGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCg                                   >  1:7841991/1‑101 (MQ=255)
                                                                    aCGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCga                                  >  1:4130545/1‑100 (MQ=255)
                                                                      gcgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTAc                                >  1:12386711/1‑101 (MQ=255)
                                                                       agcgTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACa                               >  1:4887033/2‑101 (MQ=255)
                                                                         cgTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgc                             >  1:14895327/1‑100 (MQ=255)
                                                                          ttGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGc                            >  1:4599958/2‑101 (MQ=255)
                                                                           tgaaTCTTCTCGTCGACCAGGTGGTGCATCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgcgt                           <  1:10246661/98‑4 (MQ=255)
                                                                            gCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCgg                          <  1:813583/101‑1 (MQ=255)
                                                                             aaTCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgcgtgc                         >  1:5190826/2‑96 (MQ=255)
                                                                                cTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgcgtggcct                      >  1:8529704/1‑93 (MQ=255)
                                                                                cTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgcgtggcct                      >  1:9267102/1‑93 (MQ=255)
                                                                                   gtcgtcGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGAc                   <  1:6949574/101‑1 (MQ=255)
                                                                                     cgtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCg                 >  1:4701792/1‑101 (MQ=255)
                                                                                      gtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCgg                >  1:5259181/1‑101 (MQ=255)
                                                                                      gtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCgg                >  1:10536403/1‑101 (MQ=255)
                                                                                           ccAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTc           >  1:11381051/1‑101 (MQ=255)
                                                                                            cAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCg          <  1:11181214/101‑1 (MQ=255)
                                                                                               gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa       >  1:13334516/1‑101 (MQ=255)
                                                                                               gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa       <  1:1940705/101‑1 (MQ=255)
                                                                                               gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa       >  1:6269678/1‑101 (MQ=255)
                                                                                                    tCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCg  >  1:314407/1‑101 (MQ=255)
                                                                                                    tCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCg  <  1:989091/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCG  >  NC_007779/3451631‑3451831

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.