Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,450,056 C→A 22.3% A426S (GCA→TCA)  rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,450,0560C→A22.3% 247.3 ‑0.1 ‑0.1363A426S (GCA→TCA) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  new base (42/39):  ref base (150/132):  total (192/171)
Fisher's exact test strand distribution p-value = 9.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.89e-01

Only 100 of 367 total aligned reads displayed.

GCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAATGGTG  >  NC_007779/3449956‑3450156
                                                                                                    |                                                                                                    
 cGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCa                                                                                                     >  1:4746718/1‑101 (MQ=255)
 cGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCa                                                                                                     >  1:1113356/1‑101 (MQ=255)
  gTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTCCCTTCGATCAGTACCGGTTCAAATGCCTGGATTCCCAGACGGTGCAGAGTCGGTGCAc                                                                                                    <  1:6576155/101‑1 (MQ=255)
  gTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCAc                                                                                                    <  1:6801218/101‑1 (MQ=255)
   ttATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACg                                                                                                   <  1:10126037/101‑1 (MQ=255)
   ttATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACg                                                                                                   >  1:7028907/1‑101 (MQ=255)
     atatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGt                                                                                                 <  1:8133326/101‑1 (MQ=255)
      tatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGtt                                                                                                <  1:10864420/101‑1 (MQ=255)
      tatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGtt                                                                                                >  1:10729059/1‑101 (MQ=255)
        tGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTca                                                                                              <  1:9220044/101‑1 (MQ=255)
         gCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcag                                                                                             >  1:10519647/1‑101 (MQ=255)
           cGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagca                                                                                           <  1:9163450/101‑1 (MQ=255)
            gCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGCGCAGAGTCGGTGCACGGTTcagcag                                                                                          >  1:5403993/1‑101 (MQ=255)
             cacaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGt                                                                                         >  1:12954426/1‑101 (MQ=255)
             cacaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGt                                                                                         >  1:13794180/1‑101 (MQ=255)
              acaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                                                                                        <  1:13904459/101‑1 (MQ=255)
              acaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                                                                                        >  1:13228360/1‑101 (MQ=255)
              acaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                                                                                        <  1:5906709/101‑1 (MQ=255)
               caAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAc                                                                                       <  1:695245/101‑1 (MQ=255)
                aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                                                                                      <  1:1105019/101‑1 (MQ=255)
                aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTAcc                                                                                      <  1:7411557/101‑1 (MQ=255)
                   ccAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCggg                                                                                   <  1:11926033/101‑1 (MQ=255)
                    cAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgt                                                                                  >  1:4377798/1‑101 (MQ=255)
                     aGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgtg                                                                                 >  1:172979/1‑101 (MQ=255)
                       cGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGtt                                                                               >  1:3576643/1‑101 (MQ=255)
                        gggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTc                                                                              >  1:777397/1‑101 (MQ=255)
                        gggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTc                                                                              >  1:13152709/1‑101 (MQ=255)
                          gTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcgc                                                                            >  1:10595252/1‑101 (MQ=255)
                            gCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCgg                                                                          >  1:5733952/1‑101 (MQ=255)
                              aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGAt                                                                        >  1:11285873/1‑101 (MQ=255)
                               gCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATa                                                                       <  1:4627218/101‑1 (MQ=255)
                               gCTGGATAGCTATACCTTCGATCAGGACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATa                                                                       <  1:2624144/101‑1 (MQ=255)
                                cTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATaa                                                                      >  1:10191927/1‑101 (MQ=255)
                                 tGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAAc                                                                     >  1:5786675/1‑101 (MQ=255)
                                 tGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAAc                                                                     >  1:7188188/1‑101 (MQ=255)
                                  ggATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACt                                                                    >  1:2571786/1‑101 (MQ=255)
                                   gATAGCTTTTCCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACtt                                                                   <  1:13055562/101‑1 (MQ=255)
                                    aTAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTc                                                                  <  1:11834384/101‑1 (MQ=255)
                                     tAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCg                                                                 >  1:2696737/1‑101 (MQ=255)
                                      aGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGt                                                                >  1:5616749/1‑101 (MQ=255)
                                      aGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGt                                                                >  1:2491434/1‑101 (MQ=255)
                                       gCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTc                                                               >  1:7263590/1‑101 (MQ=255)
                                       gCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTc                                                               <  1:8840256/101‑1 (MQ=255)
                                       gCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCACACGGCGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTc                                                               >  1:6185263/1‑101 (MQ=255)
                                        cTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTcc                                                              <  1:3743746/101‑1 (MQ=255)
                                           tACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAgg                                                           >  1:4566725/1‑101 (MQ=255)
                                           tACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAgg                                                           >  1:10298692/1‑101 (MQ=255)
                                           tACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAgg                                                           >  1:12637580/1‑101 (MQ=255)
                                           tACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTATTCGCGGATAACTTCGTCCAgg                                                           <  1:9940041/101‑1 (MQ=255)
                                            aCCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGa                                                          >  1:6671985/1‑101 (MQ=255)
                                               ttCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGatat                                                       >  1:12777530/1‑101 (MQ=255)
                                                tCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATc                                                      >  1:9571150/1‑101 (MQ=255)
                                                 cGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATcc                                                     <  1:9707214/101‑1 (MQ=255)
                                                   aTCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCa                                                   <  1:3810748/101‑1 (MQ=255)
                                                    tCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCaa                                                  <  1:14217347/101‑1 (MQ=255)
                                                     cAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCaaa                                                 >  1:8536345/1‑101 (MQ=255)
                                                      aGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAAc                                                <  1:1657612/101‑1 (MQ=255)
                                                      aGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAAc                                                <  1:12209727/101‑1 (MQ=255)
                                                      aGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAAc                                                >  1:3069186/1‑101 (MQ=255)
                                                      aGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAAc                                                <  1:9436607/101‑1 (MQ=255)
                                                       gTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACg                                               >  1:2345915/1‑101 (MQ=255)
                                                        tACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGa                                              >  1:4237776/1‑101 (MQ=255)
                                                        tACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGa                                              >  1:8400812/1‑101 (MQ=255)
                                                         aCCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGAc                                             >  1:1950557/1‑101 (MQ=255)
                                                         aCCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGAc                                             >  1:1490874/1‑101 (MQ=255)
                                                          ccGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACa                                            <  1:301795/101‑1 (MQ=255)
                                                            ggTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGc                                          <  1:10433525/101‑1 (MQ=255)
                                                                cAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGcttct                                      <  1:60965/101‑1 (MQ=255)
                                                                cAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGcttct                                      <  1:3644543/101‑1 (MQ=255)
                                                                   aTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcg                                   <  1:5286906/101‑1 (MQ=255)
                                                                   aTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcg                                   <  1:2893525/101‑1 (MQ=255)
                                                                    tGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcgc                                  <  1:3402915/101‑1 (MQ=255)
                                                                    tGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcgc                                  <  1:14083874/101‑1 (MQ=255)
                                                                      ccTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCgcgc                                >  1:7888641/1‑101 (MQ=255)
                                                                      ccTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCgcgc                                <  1:4046199/101‑1 (MQ=255)
                                                                       cTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCt                               >  1:14434777/1‑101 (MQ=255)
                                                                        tGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAACCGACAGCTTCTTCGCGCTc                              <  1:8948134/101‑1 (MQ=255)
                                                                        tGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTc                              <  1:7270630/101‑1 (MQ=255)
                                                                        tGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCACTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTc                              <  1:12407692/101‑1 (MQ=255)
                                                                         ggATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCa                             <  1:11583407/101‑1 (MQ=255)
                                                                          gATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCaa                            <  1:14159637/101‑1 (MQ=255)
                                                                            tACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAAcc                          <  1:7684328/101‑1 (MQ=255)
                                                                              cccAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAt                        >  1:4016764/1‑101 (MQ=255)
                                                                               ccAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAtt                       >  1:4802542/1‑101 (MQ=255)
                                                                               ccAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAtt                       <  1:10005442/101‑1 (MQ=255)
                                                                               ccAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAtt                       <  1:6409129/101‑1 (MQ=255)
                                                                                cAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAttt                      <  1:14369553/101‑1 (MQ=255)
                                                                                 aGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAtttt                     >  1:1860296/1‑101 (MQ=255)
                                                                                    cGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttctt                  >  1:3297327/1‑101 (MQ=255)
                                                                                     ggTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttcttc                 <  1:826101/101‑1 (MQ=255)
                                                                                       tGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGc               <  1:7471250/101‑1 (MQ=255)
                                                                                        gCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCa              <  1:5056476/101‑1 (MQ=255)
                                                                                          agagTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGc            >  1:11967136/1‑101 (MQ=255)
                                                                                           gagTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCt           <  1:9298150/101‑1 (MQ=255)
                                                                                           gagTCGGTGAACGGTTCAGCAGTACCAGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCt           <  1:7972271/101‑1 (MQ=255)
                                                                                            agTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCtt          <  1:13600237/101‑1 (MQ=255)
                                                                                               cGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTaa       >  1:1246144/1‑101 (MQ=255)
                                                                                                 gTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAAtg     >  1:10376937/1‑101 (MQ=255)
                                                                                                 gTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAAtg     <  1:1337915/101‑1 (MQ=255)
                                                                                                   gCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAAtggt   <  1:9827813/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAATGGTG  >  NC_007779/3449956‑3450156

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.