Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,682 C→G 24.6% S1252T (AGC→ACC)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,6820C→G24.6% 274.4 ‑0.5 ‑0.1317S1252T (AGC→ACC) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (41/37):  ref base (116/123):  total (157/160)
Fisher's exact test strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.97e-01

Only 100 of 331 total aligned reads displayed.

GTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTC  >  NC_007779/3451582‑3451782
                                                                                                    |                                                                                                    
 tATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCt                                                                                                     <  1:6628393/101‑1 (MQ=255)
  aTGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTg                                                                                                    >  1:7047666/1‑101 (MQ=255)
  aTGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTg                                                                                                    <  1:8953693/101‑1 (MQ=255)
   tgctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGt                                                                                                   <  1:8457017/101‑1 (MQ=255)
    gctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTa                                                                                                  >  1:13892874/1‑101 (MQ=255)
    gctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTa                                                                                                  <  1:10341831/101‑1 (MQ=255)
     ctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTaa                                                                                                 >  1:6772732/1‑101 (MQ=255)
      tgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACTAGGCTGTAAg                                                                                                >  1:493969/1‑101 (MQ=255)
      tgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAg                                                                                                >  1:4317031/1‑101 (MQ=255)
        cgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGaa                                                                                              <  1:1069739/101‑1 (MQ=255)
            gTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCgg                                                                                          >  1:176965/1‑101 (MQ=255)
            gTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCgg                                                                                          <  1:9022672/101‑1 (MQ=255)
             tATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGt                                                                                         <  1:5573592/101‑1 (MQ=255)
              aTGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTg                                                                                        <  1:11384985/101‑1 (MQ=255)
              aTGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTg                                                                                        <  1:2697894/101‑1 (MQ=255)
                   tCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcg                                                                                   <  1:5784597/101‑1 (MQ=255)
                     cattgcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAAcgag                                                                                 <  1:342390/97‑3 (MQ=255)
                     cattgcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAAcgag                                                                                 <  1:9717329/97‑3 (MQ=255)
                     cAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcgcg                                                                                 <  1:3826113/101‑1 (MQ=255)
                       gcgcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcgcgcg                                                                               >  1:6962740/1‑101 (MQ=255)
                       gcgcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcgcgcg                                                                               <  1:1149819/101‑1 (MQ=255)
                        cgcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGt                                                                              <  1:2089599/101‑1 (MQ=255)
                        cgcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGt                                                                              <  1:13734448/101‑1 (MQ=255)
                         gcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTg                                                                             >  1:3134202/1‑101 (MQ=255)
                         gcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTg                                                                             >  1:7441142/1‑101 (MQ=255)
                          cccACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGc                                                                            <  1:4480201/101‑1 (MQ=255)
                              cacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCt                                                                        >  1:230809/1‑101 (MQ=255)
                              cacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGGGAATCt                                                                        >  1:12150527/1‑101 (MQ=255)
                               acTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCtt                                                                       <  1:13568404/101‑1 (MQ=255)
                                cctCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCttc                                                                      >  1:257768/3‑100 (MQ=255)
                                 ctCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCttct                                                                     <  1:11111831/100‑3 (MQ=255)
                                  tCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtc                                                                    <  1:11065400/101‑1 (MQ=255)
                                   ccATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcg                                                                   >  1:12997649/1‑101 (MQ=255)
                                     aTCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcgtc                                                                 <  1:11942692/101‑1 (MQ=255)
                                     aTCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcgtc                                                                 >  1:6359271/1‑101 (MQ=255)
                                      tctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGtcgtcg                                                                >  1:2948123/1‑101 (MQ=255)
                                      tctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGtcgtcg                                                                <  1:6314551/101‑1 (MQ=255)
                                       ctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGa                                                               >  1:404467/1‑101 (MQ=255)
                                       ctcACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGa                                                               >  1:3184077/1‑101 (MQ=255)
                                        tcaccGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGAc                                                              >  1:3293561/4‑101 (MQ=255)
                                         caccGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGAcc                                                             <  1:1609645/99‑1 (MQ=255)
                                            cGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAgg                                                          <  1:10925293/101‑1 (MQ=255)
                                             gAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGGCCAggt                                                         >  1:6077352/1‑101 (MQ=255)
                                             gAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggt                                                         >  1:13265252/1‑101 (MQ=255)
                                             gAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggt                                                         >  1:13152715/1‑101 (MQ=255)
                                             gAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggt                                                         <  1:6256202/101‑1 (MQ=255)
                                              aaaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtg                                                        >  1:7287453/1‑101 (MQ=255)
                                              aaaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtg                                                        >  1:1570429/1‑101 (MQ=255)
                                                gcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGAGGAACGAGCGTGAATCTTCTCGTCGACCAggtggt                                                      <  1:9201128/100‑1 (MQ=255)
                                                aCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtggt                                                      >  1:12718096/1‑101 (MQ=255)
                                                 cGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGtt                                                     >  1:9551890/1‑101 (MQ=255)
                                                 cGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtggtg                                                     <  1:1239049/101‑2 (MQ=255)
                                                  gCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTc                                                    >  1:533303/1‑101 (MQ=255)
                                                   cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                   >  1:1988223/1‑101 (MQ=255)
                                                    tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                  >  1:469732/1‑101 (MQ=255)
                                                    tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                  >  1:9167942/1‑101 (MQ=255)
                                                     gACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGt                                                 >  1:7545805/1‑101 (MQ=255)
                                                      accaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGtt                                                <  1:6462836/101‑1 (MQ=255)
                                                      accaccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAgct                                                <  1:2351096/101‑3 (MQ=255)
                                                        caccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTc                                              >  1:12576413/1‑101 (MQ=255)
                                                         accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                             <  1:12270506/101‑1 (MQ=255)
                                                         accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                             <  1:10066557/101‑1 (MQ=255)
                                                            gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                          >  1:11892772/1‑101 (MQ=255)
                                                            gAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGca                                          <  1:5994748/101‑1 (MQ=255)
                                                             aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                         <  1:8015771/101‑1 (MQ=255)
                                                             aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                         <  1:2760824/101‑1 (MQ=255)
                                                             aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                         >  1:8827308/1‑101 (MQ=255)
                                                             aaCTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcat                                         <  1:399451/101‑1 (MQ=255)
                                                                tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgta                                      <  1:7651485/101‑1 (MQ=255)
                                                                 gagCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtac                                     <  1:8720197/99‑1 (MQ=255)
                                                                  tgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtaca                                    >  1:787690/1‑101 (MQ=255)
                                                                     cTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatgt                                 >  1:6207111/1‑101 (MQ=255)
                                                                      ttACCACCCAGCTGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACACGCGTGCATCTTGTCGTCGACCAGGTGGTTCCGTTTCAGcatgtacatgta                                >  1:13844075/1‑101 (MQ=255)
                                                                      ttACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatgta                                >  1:13356586/1‑101 (MQ=255)
                                                                       tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                               >  1:9477206/1‑101 (MQ=255)
                                                                         ccaccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGcc                             >  1:179145/1‑101 (MQ=255)
                                                                         ccaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAAcc                             >  1:8296005/1‑101 (MQ=255)
                                                                           accCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccca                           <  1:14244795/101‑3 (MQ=255)
                                                                             ccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACg                         >  1:2877309/1‑101 (MQ=255)
                                                                              cAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCTTGTACATGTAACCAACgg                        >  1:8031205/1‑101 (MQ=255)
                                                                              cAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACgg                        <  1:13119946/101‑1 (MQ=255)
                                                                              cAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACgg                        <  1:2435556/101‑1 (MQ=255)
                                                                               aGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccga                       <  1:11250974/101‑4 (MQ=255)
                                                                               aGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGt                       >  1:2186104/1‑101 (MQ=255)
                                                                               aGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGt                       <  1:10508078/101‑1 (MQ=255)
                                                                                gCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGtt                      >  1:12405433/1‑101 (MQ=255)
                                                                                gCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGtt                      <  1:834280/101‑1 (MQ=255)
                                                                                  ggCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaac                    <  1:8185495/101‑7 (MQ=255)
                                                                                   gctggtGAGTACCCAGGCTGTAAGAACCGGTGGAACGCGCCTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAGCGGTTAcc                   <  1:12505066/101‑1 (MQ=255)
                                                                                    ctgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCg                  >  1:2674152/1‑101 (MQ=255)
                                                                                     tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg                 >  1:7219196/1‑101 (MQ=255)
                                                                                      gctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGa                >  1:5870036/1‑101 (MQ=255)
                                                                                       ctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGAc               >  1:2230520/1‑101 (MQ=255)
                                                                                       ctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGAc               >  1:526108/1‑101 (MQ=255)
                                                                                       ctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGAc               <  1:3617775/101‑1 (MQ=255)
                                                                                          aGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTGACCAACGGTTACCGGACGCt            <  1:11157402/101‑1 (MQ=255)
                                                                                            tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg          >  1:3896560/1‑101 (MQ=255)
                                                                                             aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa         >  1:13983432/1‑101 (MQ=255)
                                                                                                 aGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTg     >  1:4107294/1‑101 (MQ=255)
                                                                                                   gCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCCACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGTACGCTCGAACTGtt   <  1:2997053/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTC  >  NC_007779/3451582‑3451782

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.