Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,450,056 C→A 23.6% 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→A23.6% 159.3 ‑0.8 ‑1.9233A426S (GCA→TCA) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  new base (27/28):  ref base (100/78):  total (127/106)
Fisher's exact test strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.33e-02

Only 100 of 235 total aligned reads displayed.

GCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAATG  >  NC_007779/3449956‑3450153
                                                                                                    |                                                                                                 
gCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGc                                                                                                   <  1:9645551/101‑1 (MQ=255)
  gTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCAc                                                                                                 >  1:5807318/1‑101 (MQ=255)
      tatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGtt                                                                                             >  1:8860742/1‑101 (MQ=255)
       atGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTc                                                                                            >  1:5157714/1‑101 (MQ=255)
         gCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcag                                                                                          <  1:1866208/101‑1 (MQ=255)
         gCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGGTcag                                                                                          <  1:3283582/101‑1 (MQ=255)
          ccGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagc                                                                                         >  1:3775787/1‑101 (MQ=255)
           cGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagca                                                                                        <  1:1064282/101‑1 (MQ=255)
           cGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagca                                                                                        >  1:9320046/1‑101 (MQ=255)
           cGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagca                                                                                        >  1:3674724/1‑101 (MQ=255)
               caAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTAc                                                                                    >  1:4932317/1‑101 (MQ=255)
                aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                                                                                   >  1:2564833/1‑101 (MQ=255)
                aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                                                                                   <  1:3706928/101‑1 (MQ=255)
                   ccAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCggg                                                                                <  1:1832003/101‑1 (MQ=255)
                    cAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgt                                                                               >  1:1305161/1‑101 (MQ=255)
                     aGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgtg                                                                              >  1:1669288/1‑101 (MQ=255)
                     aGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGgtg                                                                              >  1:5781919/1‑101 (MQ=255)
                      gCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgtgt                                                                             >  1:4404683/1‑101 (MQ=255)
                       cGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGtt                                                                            <  1:7997131/101‑1 (MQ=255)
                       cGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGtt                                                                            >  1:1007959/1‑101 (MQ=255)
                         ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcg                                                                          >  1:5546564/1‑101 (MQ=255)
                         ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTcg                                                                          >  1:9020894/1‑101 (MQ=255)
                          gTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTcgc                                                                         <  1:8523408/101‑1 (MQ=255)
                            gCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCgg                                                                       >  1:9383359/1‑101 (MQ=255)
                             cAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGa                                                                      <  1:8189082/101‑1 (MQ=255)
                              aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGAt                                                                     >  1:3019290/1‑101 (MQ=255)
                              aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGAt                                                                     >  1:5445224/1‑101 (MQ=255)
                              aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGAt                                                                     <  1:7951987/101‑1 (MQ=255)
                               gCTGGATAGCTTTTCCTTCGATCAGTACAGGTTCAAATGCCTGGATGCCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATa                                                                    <  1:3859046/101‑1 (MQ=255)
                               gCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATa                                                                    <  1:4740110/101‑1 (MQ=255)
                                   gATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACtt                                                                >  1:7353767/1‑101 (MQ=255)
                                       gCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTc                                                            <  1:6722112/101‑1 (MQ=255)
                                        cTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTcc                                                           >  1:7171833/1‑101 (MQ=255)
                                         tttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCa                                                          <  1:7308980/101‑1 (MQ=255)
                                         tttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCa                                                          >  1:1506137/1‑101 (MQ=255)
                                         tttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCa                                                          >  1:3169837/1‑101 (MQ=255)
                                          ttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAg                                                         >  1:169751/1‑101 (MQ=255)
                                          ttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAg                                                         <  1:8908533/101‑1 (MQ=255)
                                             ccTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGat                                                      >  1:5689798/1‑101 (MQ=255)
                                              cTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGata                                                     <  1:2914026/101‑1 (MQ=255)
                                              cTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGata                                                     <  1:6866164/101‑1 (MQ=255)
                                               ttCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGatat                                                    <  1:8446176/101‑1 (MQ=255)
                                               ttCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGatat                                                    >  1:7687932/1‑101 (MQ=255)
                                                tCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATc                                                   >  1:8415809/1‑101 (MQ=255)
                                                 cGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATcc                                                  >  1:960625/1‑101 (MQ=255)
                                                 cGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATcc                                                  >  1:3372073/1‑101 (MQ=255)
                                                 cGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGCGTTCGCGGATAACTTCGTCCAGGATATcc                                                  >  1:4792941/1‑101 (MQ=255)
                                                  gATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATccc                                                 <  1:3991039/101‑1 (MQ=255)
                                                  gATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATccc                                                 >  1:570464/1‑101 (MQ=255)
                                                   aTCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCa                                                >  1:1637830/1‑101 (MQ=255)
                                                      aGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAAc                                             >  1:2378324/1‑101 (MQ=255)
                                                      aGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAAc                                             >  1:5047047/1‑101 (MQ=255)
                                                        tACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGa                                           >  1:2416353/1‑101 (MQ=255)
                                                         aCCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGAc                                          >  1:5474208/1‑101 (MQ=255)
                                                         aCCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGAc                                          <  1:1736291/101‑1 (MQ=255)
                                                          ccGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACa                                         >  1:6032840/1‑101 (MQ=255)
                                                            ggTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGc                                       >  1:4993002/1‑101 (MQ=255)
                                                             gTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGct                                      >  1:8626091/1‑101 (MQ=255)
                                                             gTTCAAATGCCTGGATACCCAGACGGGGCAGAGTCGGTGAACGGTTCAGCATTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGct                                      >  1:1165347/1‑101 (MQ=255)
                                                              ttCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTACACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGctt                                     <  1:6190212/101‑1 (MQ=255)
                                                               tCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGcttt                                    >  1:2259379/1‑100 (MQ=255)
                                                                 aaaTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGTTGTTCGCGGATAACATCGTCCAGGATATCCCAAACGACAGcttctt                                  >  1:6916426/1‑101 (MQ=255)
                                                                 aaaTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGcttctt                                  >  1:3115817/1‑101 (MQ=255)
                                                                 aaaTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGcttctt                                  <  1:8934631/101‑1 (MQ=255)
                                                                  aaTGCCTGGATACCGAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCttcttc                                 >  1:1374486/1‑101 (MQ=255)
                                                                  aaTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCttcttc                                 <  1:8293418/101‑1 (MQ=255)
                                                                  aaTGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCttcttc                                 <  1:1211527/101‑1 (MQ=255)
                                                                  aaTGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCttcttc                                 <  1:2147200/101‑1 (MQ=255)
                                                                   aTGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGTGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcg                                <  1:4044689/101‑1 (MQ=255)
                                                                    tGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcgc                               <  1:8719507/101‑1 (MQ=255)
                                                                    tGCCTGGATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcgc                               <  1:6281741/101‑1 (MQ=255)
                                                                     gCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTcgcg                              <  1:8095440/101‑1 (MQ=255)
                                                                      ccTGGATGCCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCgcgc                             <  1:29026/101‑1 (MQ=255)
                                                                      ccTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCgcgc                             >  1:7233073/1‑101 (MQ=255)
                                                                        tGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTc                           >  1:4163697/1‑101 (MQ=255)
                                                                          gATACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGCTATCCCAAACGACAGCTTCTTCGCGCTCaa                         <  1:9304673/101‑1 (MQ=255)
                                                                            tACCCAGACGGTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAAcc                       >  1:7870119/1‑101 (MQ=255)
                                                                             aCCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCTGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCa                      <  1:8602760/101‑1 (MQ=255)
                                                                             aCCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCa                      >  1:239571/1‑101 (MQ=255)
                                                                              cccAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAt                     >  1:2675722/1‑101 (MQ=255)
                                                                               ccAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAtt                    <  1:5336272/101‑1 (MQ=255)
                                                                                cAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAttt                   <  1:3077492/101‑1 (MQ=255)
                                                                                 aGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCAtttt                  <  1:945211/101‑1 (MQ=255)
                                                                                  gACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttc                 <  1:4306568/101‑1 (MQ=255)
                                                                                   aCGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttct                <  1:4629130/101‑1 (MQ=255)
                                                                                    cGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttctt               >  1:5967585/1‑101 (MQ=255)
                                                                                    cGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttctt               <  1:9260078/101‑1 (MQ=255)
                                                                                     ggTGCAGAGTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTttcttc              >  1:1563510/1‑101 (MQ=255)
                                                                                       tGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGc            >  1:6624052/1‑101 (MQ=255)
                                                                                       tGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGAAAGCTTCTTCGCGCTCAACCATTttcttaga            >  1:342972/1‑98 (MQ=255)
                                                                                        gCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCa           <  1:2067329/101‑1 (MQ=255)
                                                                                          agagTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGc         <  1:5892192/101‑1 (MQ=255)
                                                                                           gagTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCt        <  1:394435/101‑1 (MQ=255)
                                                                                           gagTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCt        >  1:56017/1‑101 (MQ=255)
                                                                                            agTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCtt       >  1:639571/1‑101 (MQ=255)
                                                                                             gTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCttt      <  1:5221518/101‑1 (MQ=255)
                                                                                             gTCGGTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCttt      >  1:7752646/1‑101 (MQ=255)
                                                                                                ggTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAAt   >  1:2076508/1‑101 (MQ=255)
                                                                                                ggTGAACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAAt   <  1:7396295/101‑1 (MQ=255)
                                                                                                 gTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGGCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAAtg  >  1:1997572/1‑101 (MQ=255)
                                                                                                    |                                                                                                 
GCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAAACGACAGCTTCTTCGCGCTCAACCATTTTCTTCGCAGCTTTAATG  >  NC_007779/3449956‑3450153

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