Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,702 C→A 26.2% A1245A (GCG→GCT)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,7020C→A26.2% 225.0 ‑2.8 ‑0.0289A1245A (GCG→GCT) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (45/31):  ref base (98/115):  total (143/146)
Fisher's exact test strand distribution p-value = 6.11e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

Only 100 of 309 total aligned reads displayed.

CCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTA  >  NC_007779/3451602‑3451800
                                                                                                    |                                                                                                  
 cAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcgcg                                                                                                   <  1:71064/101‑1 (MQ=255)
   ttgcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAAcgagcg                                                                                                 >  1:13981927/3‑101 (MQ=255)
   ttgcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAAcgagcg                                                                                                 >  1:12235822/3‑101 (MQ=255)
   ttgcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAAcgagcg                                                                                                 <  1:8963679/99‑1 (MQ=255)
     gcCTACACTTCCATCTCCCCGAAACGCTGACCAGCGAACTGTGGCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTg                                                                                               <  1:4179132/101‑1 (MQ=255)
     gcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTg                                                                                               <  1:12292559/101‑1 (MQ=255)
     gcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTg                                                                                               >  1:1134610/1‑101 (MQ=255)
        cacaCTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCAt                                                                                            <  1:10917561/101‑1 (MQ=255)
        cacaCCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAAt                                                                                            >  1:6287020/1‑101 (MQ=255)
          cacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCt                                                                                          >  1:1604873/1‑101 (MQ=255)
           acTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCtt                                                                                         >  1:464604/1‑101 (MQ=255)
             ttCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgt                                                                                       >  1:7269905/1‑101 (MQ=255)
              tCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtc                                                                                      <  1:1298045/101‑1 (MQ=255)
              tCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTctc                                                                                      <  1:393885/101‑1 (MQ=255)
               ccATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcg                                                                                     >  1:5349500/1‑101 (MQ=255)
                 aTCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTctcgtc                                                                                   >  1:10242954/1‑101 (MQ=255)
                   ctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGa                                                                                 >  1:10748261/1‑101 (MQ=255)
                   ctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGCGGAACGCGCGTGCATCTTGTCGTCGa                                                                                 >  1:11465097/1‑101 (MQ=255)
                   ctcACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGa                                                                                 >  1:2681102/1‑101 (MQ=255)
                    tcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGAc                                                                                >  1:1247/1‑101 (MQ=255)
                     ccccGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGAcc                                                                               >  1:14010071/1‑101 (MQ=255)
                      cccGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCa                                                                              <  1:11577214/101‑1 (MQ=255)
                        cGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAgg                                                                            >  1:9613389/1‑101 (MQ=255)
                         gAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggt                                                                           <  1:11235216/101‑1 (MQ=255)
                         gAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggt                                                                           <  1:4020991/101‑1 (MQ=255)
                           agcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGCGGAACGAGCGTGAATCTTCTCGTCGACCAggtgg                                                                         >  1:1439023/3‑101 (MQ=255)
                           aaCGCTGACCACCGAACTGTGCCTTACCACCCACCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtgg                                                                         >  1:475684/1‑101 (MQ=255)
                            gcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtggt                                                                        >  1:7433112/2‑101 (MQ=255)
                             cGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGtt                                                                       <  1:8002396/101‑1 (MQ=255)
                               cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                                     <  1:1787216/101‑1 (MQ=255)
                               cTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCa                                                                     <  1:10086752/101‑1 (MQ=255)
                                tGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAg                                                                    <  1:3781170/101‑1 (MQ=255)
                                tGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAg                                                                    >  1:7528262/1‑101 (MQ=255)
                                tGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAg                                                                    <  1:815107/101‑1 (MQ=255)
                                tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                    <  1:14249938/101‑1 (MQ=255)
                                 gACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAgc                                                                   >  1:6604200/1‑100 (MQ=255)
                                   ccaccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGctt                                                                 <  1:14185281/101‑1 (MQ=255)
                                    caccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGGGCAGCTTc                                                                >  1:11920290/1‑101 (MQ=255)
                                     accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                                               >  1:1923510/1‑101 (MQ=255)
                                      ccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAg                                                              <  1:3314354/101‑1 (MQ=255)
                                      ccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTCGTTCAGTTTCAg                                                              <  1:8519388/101‑1 (MQ=255)
                                      ccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAg                                                              <  1:3980810/101‑1 (MQ=255)
                                       cGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGc                                                             <  1:12688840/101‑1 (MQ=255)
                                        gAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGca                                                            >  1:2256902/1‑101 (MQ=255)
                                          aCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatg                                                          <  1:5097291/101‑1 (MQ=255)
                                          aCTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatg                                                          <  1:2122096/101‑1 (MQ=255)
                                          aCTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatg                                                          <  1:10484858/101‑1 (MQ=255)
                                            tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgta                                                        <  1:1690037/101‑1 (MQ=255)
                                            tgagCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgta                                                        <  1:13727088/98‑1 (MQ=255)
                                             gtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtac                                                       >  1:4960520/1‑101 (MQ=255)
                                              tgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtaca                                                      <  1:4417195/101‑1 (MQ=255)
                                              tgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGGCGACCAGGTGGTTCAGTTTCAGcatgtaca                                                      <  1:6226719/101‑1 (MQ=255)
                                               gCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacat                                                     <  1:9819563/101‑1 (MQ=255)
                                               gCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacat                                                     <  1:10581515/101‑1 (MQ=255)
                                                ccTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatg                                                    <  1:2209161/101‑1 (MQ=255)
                                                    accaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAAc                                                >  1:9054136/1‑101 (MQ=255)
                                                    accaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAAc                                                <  1:6555346/101‑1 (MQ=255)
                                                      caccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccc                                              >  1:4243213/1‑100 (MQ=255)
                                                       accCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCaa                                             <  1:1732967/101‑1 (MQ=255)
                                                        cccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAAc                                            <  1:11848960/101‑1 (MQ=255)
                                                         ccAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAcc                                           >  1:5932626/1‑100 (MQ=255)
                                                          cAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACgg                                          >  1:7926995/1‑101 (MQ=255)
                                                           aGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccga                                         >  1:7741964/1‑98 (MQ=255)
                                                           aGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGt                                         >  1:3420870/1‑101 (MQ=255)
                                                             cGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTa                                       <  1:10013960/101‑1 (MQ=255)
                                                             cGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTa                                       >  1:13181590/1‑101 (MQ=255)
                                                               gctgctGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaact                                     <  1:3670418/101‑8 (MQ=255)
                                                               gctgctGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaacc                                     <  1:8229951/101‑8 (MQ=255)
                                                               gctgctGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAcc                                     <  1:5486589/101‑1 (MQ=255)
                                                                ctgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCg                                    <  1:13641699/101‑1 (MQ=255)
                                                                 tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg                                   <  1:6737811/101‑1 (MQ=255)
                                                                  gctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGa                                  >  1:8734832/1‑101 (MQ=255)
                                                                  gctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGa                                  >  1:12783301/1‑101 (MQ=255)
                                                                   ctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGAc                                 <  1:6476347/101‑1 (MQ=255)
                                                                    tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg                                >  1:10185353/1‑101 (MQ=255)
                                                                    tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg                                <  1:6965133/101‑1 (MQ=255)
                                                                       gTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTc                             >  1:7651068/1‑101 (MQ=255)
                                                                       gTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTc                             <  1:577646/101‑1 (MQ=255)
                                                                        tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                            <  1:6032842/101‑1 (MQ=255)
                                                                        tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                            >  1:8914955/1‑101 (MQ=255)
                                                                        tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                            >  1:13464533/1‑101 (MQ=255)
                                                                         aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                           >  1:11037892/1‑101 (MQ=255)
                                                                           ccAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAAc                         >  1:2498214/1‑101 (MQ=255)
                                                                              ggCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGt                      >  1:13057368/1‑101 (MQ=255)
                                                                              ggCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGt                      >  1:13902832/1‑101 (MQ=255)
                                                                                 tGTAAGAACCGGTGGAACGCGTGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCa                   >  1:8637725/1‑101 (MQ=255)
                                                                                 tGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCa                   <  1:8345737/101‑1 (MQ=255)
                                                                                 tGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCa                   <  1:4522686/101‑1 (MQ=255)
                                                                                  gTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCAc                  >  1:4502356/1‑101 (MQ=255)
                                                                                   tAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCAcc                 >  1:7806075/1‑101 (MQ=255)
                                                                                     aGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAg               <  1:13273779/101‑1 (MQ=255)
                                                                                     aGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAg               <  1:9531566/101‑1 (MQ=255)
                                                                                      gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt              <  1:10434322/101‑1 (MQ=255)
                                                                                        aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc            <  1:5211072/101‑1 (MQ=255)
                                                                                         ccGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCg           >  1:9685997/1‑101 (MQ=255)
                                                                                          cGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCgg          >  1:12729398/1‑100 (MQ=255)
                                                                                            gTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGAcc        >  1:9414734/1‑101 (MQ=255)
                                                                                              ggAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCAt      <  1:9200005/101‑1 (MQ=255)
                                                                                              ggAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccgt      >  1:3404286/1‑99 (MQ=255)
                                                                                                 aCGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGt   <  1:11172596/101‑1 (MQ=255)
                                                                                                    |                                                                                                  
CCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTA  >  NC_007779/3451602‑3451800

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.