Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,200,197 A→G 14.9% K271E (AAA→GAA)  rpoA → RNA polymerase subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,200,1970A→G14.9% 145.5 ‑1.0 ‑0.2323K271E (AAA→GAA) rpoARNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  new base (22/26):  ref base (147/128):  total (169/154)
Fisher's exact test strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.92e-01

Only 100 of 345 total aligned reads displayed.

ACGTCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTACTGAGATTAA  >  NC_007779/4200097‑4200297
                                                                                                    |                                                                                                    
aCGTCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTa                                                                                                      <  1:2609682/101‑1 (MQ=255)
 cGTCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTaa                                                                                                     >  1:10265171/1‑101 (MQ=255)
  gTCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTgaa                                                                                                    >  1:8483938/1‑101 (MQ=255)
  gTCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTaaa                                                                                                    >  1:10838356/1‑101 (MQ=255)
   tCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAg                                                                                                   >  1:6744460/1‑101 (MQ=255)
   tCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAg                                                                                                   <  1:8743326/101‑1 (MQ=255)
        cTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGaa                                                                                              >  1:13738863/1‑101 (MQ=255)
         tGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAg                                                                                             >  1:4187437/1‑101 (MQ=255)
         tGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAg                                                                                             <  1:12015840/101‑1 (MQ=255)
            agagAAAGAAGAGAAACCAGAGTTCGATCCGATGCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTCAAGCAGAAGCTa                                                                                          >  1:9742977/4‑101 (MQ=255)
            aGTGAAAGATGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTa                                                                                          >  1:12515680/1‑101 (MQ=255)
            aGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTa                                                                                          >  1:9249185/1‑101 (MQ=255)
              tGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATc                                                                                        <  1:3033883/101‑1 (MQ=255)
              tGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATc                                                                                        >  1:8331239/1‑101 (MQ=255)
              tGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATc                                                                                        <  1:5181557/101‑1 (MQ=255)
               gAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATcc                                                                                       >  1:3972807/1‑101 (MQ=255)
                aaaGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCa                                                                                      <  1:8552803/101‑1 (MQ=255)
                aaaGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCa                                                                                      >  1:6890708/1‑101 (MQ=255)
                   gaagaGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACta                                                                                   <  1:5676417/101‑1 (MQ=255)
                   gaagaGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCGCTGCTAACTGCCTTAAAGCAGAAGCTATCCACta                                                                                   <  1:12796380/101‑1 (MQ=255)
                    aagtGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACtat                                                                                  <  1:7077241/101‑1 (MQ=255)
                    aagaGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACtat                                                                                  >  1:1332281/1‑101 (MQ=255)
                      gagaAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGGCTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTatat                                                                                <  1:4272895/101‑1 (MQ=255)
                       agaAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGTCTTAAAGCAGAAGCTATCCACTATATc                                                                               <  1:14486127/101‑1 (MQ=255)
                        gaAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCg                                                                              <  1:10068361/101‑1 (MQ=255)
                        gaAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCg                                                                              >  1:4593476/1‑101 (MQ=255)
                        gaAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCg                                                                              >  1:13424086/1‑101 (MQ=255)
                         aaaCCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCgg                                                                             >  1:8429633/1‑101 (MQ=255)
                           aCCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTg                                                                           <  1:6322505/101‑1 (MQ=255)
                             cAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGAt                                                                         <  1:6114251/101‑1 (MQ=255)
                               gagTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCt                                                                       >  1:2403855/1‑101 (MQ=255)
                                   tCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTa                                                                   <  1:9140897/101‑1 (MQ=255)
                                   tCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTa                                                                   >  1:5051072/1‑101 (MQ=255)
                                   tCGATCCGATCCTGCTGCGCCCTGTTGACGACCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTa                                                                   >  1:10215145/1‑101 (MQ=255)
                                    cgatccgatcCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTAc                                                                  <  1:1313729/101‑1 (MQ=255)
                                     gatccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACa                                                                 >  1:5520191/1‑101 (MQ=255)
                                       tccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGc                                                               <  1:3809686/101‑1 (MQ=255)
                                       tccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGc                                                               <  1:964167/101‑1 (MQ=255)
                                        ccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCg                                                              <  1:5425804/101‑1 (MQ=255)
                                        ccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCg                                                              >  1:2131180/1‑101 (MQ=255)
                                        ccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCg                                                              >  1:13864195/1‑101 (MQ=255)
                                        ccgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCg                                                              <  1:7695018/101‑1 (MQ=255)
                                         cgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGt                                                             <  1:7202610/101‑1 (MQ=255)
                                         cgatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGt                                                             >  1:10519772/1‑101 (MQ=255)
                                          gatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTa                                                            >  1:6485524/1‑101 (MQ=255)
                                          gatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTa                                                            >  1:6617706/1‑101 (MQ=255)
                                          gatccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTa                                                            >  1:6996468/1‑101 (MQ=255)
                                           atccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTAc                                                           >  1:7541987/1‑101 (MQ=255)
                                            tccTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTAcc                                                          >  1:4156909/1‑101 (MQ=255)
                                              ctgctgCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGa                                                        <  1:9312649/101‑1 (MQ=255)
                                              ctgctgCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGa                                                        >  1:7396097/1‑101 (MQ=255)
                                               tgctgcGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAg                                                       <  1:2934570/101‑1 (MQ=255)
                                                gctgcGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAgg                                                      >  1:197013/1‑101 (MQ=255)
                                                 ctgcGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGt                                                     >  1:9562168/1‑101 (MQ=255)
                                                 ctgcGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGt                                                     >  1:11447300/1‑101 (MQ=255)
                                                  tgcGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGtt                                                    >  1:13778560/1‑101 (MQ=255)
                                                   gcgcCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTg                                                   <  1:11173114/101‑1 (MQ=255)
                                                    cgcCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGa                                                  <  1:14152471/101‑1 (MQ=255)
                                                     gcCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAg                                                 <  1:8888382/101‑1 (MQ=255)
                                                     gcCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAg                                                 <  1:4484584/101‑1 (MQ=255)
                                                     gcCCTGTTGACGATCGGGAATTGACTGTCCGCGGTGCTAACTGCCTTAAAGCAGAGGCTACCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAg                                                 <  1:3575494/101‑1 (MQ=255)
                                                      cccTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGc                                                <  1:12893313/101‑1 (MQ=255)
                                                      cccTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGc                                                <  1:340241/101‑1 (MQ=255)
                                                          gTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCt                                            >  1:12572544/1‑101 (MQ=255)
                                                           ttGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCtt                                           <  1:5281966/101‑1 (MQ=255)
                                                           ttGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCtt                                           >  1:9431705/1‑101 (MQ=255)
                                                            tGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTa                                          >  1:14047169/1‑101 (MQ=255)
                                                            tGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTa                                          <  1:11976456/101‑1 (MQ=255)
                                                            tGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTa                                          <  1:507056/101‑1 (MQ=255)
                                                            tGACGATCGGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTa                                          <  1:3247620/101‑1 (MQ=255)
                                                             gACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTaa                                         <  1:1515974/101‑1 (MQ=255)
                                                                 aTCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACg                                     <  1:14111640/101‑1 (MQ=255)
                                                                 aTCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACg                                     <  1:13461565/101‑1 (MQ=255)
                                                                  tCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGc                                    <  1:14398462/101‑1 (MQ=255)
                                                                    tGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCt                                  <  1:11340089/101‑1 (MQ=255)
                                                                     ggAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTa                                 >  1:6002640/1‑101 (MQ=255)
                                                                        aTTGACTGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAAcc                              <  1:11037149/101‑1 (MQ=255)
                                                                           gACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTg                           <  1:9806673/101‑1 (MQ=255)
                                                                              tGTCCGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTa                        <  1:11921464/101‑1 (MQ=255)
                                                                               gTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaa                       <  1:13179744/101‑1 (MQ=255)
                                                                               gTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaa                       >  1:13590086/1‑101 (MQ=255)
                                                                                tCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaa                      <  1:4418035/101‑1 (MQ=255)
                                                                                tCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaa                      <  1:11675057/101‑1 (MQ=255)
                                                                                 ccGCTCTGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaaa                     <  1:10698303/101‑1 (MQ=255)
                                                                                 ccGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaaa                     >  1:2316639/1‑101 (MQ=255)
                                                                                  cGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaaaa                    >  1:9405387/1‑101 (MQ=255)
                                                                                  cGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaaaa                    >  1:3147411/1‑101 (MQ=255)
                                                                                   gCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTaaaaaa                   >  1:1804126/1‑101 (MQ=255)
                                                                                     tctGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAAtc                 <  1:826517/101‑1 (MQ=255)
                                                                                       tGCTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAAtctc               >  1:10367683/1‑101 (MQ=255)
                                                                                        gCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATctct              >  1:7308415/1‑101 (MQ=255)
                                                                                         cTAACTGCCTTGAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCtt             >  1:493323/1‑101 (MQ=255)
                                                                                         cTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCtt             <  1:8973912/101‑1 (MQ=255)
                                                                                          tAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTa            >  1:1207864/1‑101 (MQ=255)
                                                                                           aaCTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGACTAACCTTGGTAAAAAATCTCTTAc           >  1:4735132/1‑101 (MQ=255)
                                                                                            aCTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTACt          >  1:8237355/1‑101 (MQ=255)
                                                                                              tGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTACTga        >  1:7961364/1‑101 (MQ=255)
                                                                                              tGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTACTga        >  1:1564521/1‑101 (MQ=255)
                                                                                                ccTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAGAATCTCTTACTgaga      >  1:8686309/1‑101 (MQ=255)
                                                                                                 cTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTACTGAGAt     <  1:7576897/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ACGTCAGCCTGAAGTGAAAGAAGAGAAACCAGAGTTCGATCCGATCCTGCTGCGCCCTGTTGACGATCTGGAATTGACTGTCCGCTCTGCTAACTGCCTTAAAGCAGAAGCTATCCACTATATCGGTGATCTGGTACAGCGTACCGAGGTTGAGCTCCTTAAAACGCCTAACCTTGGTAAAAAATCTCTTACTGAGATTAA  >  NC_007779/4200097‑4200297

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.