Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,726 G→T 4.4% H1237Q (CAC→CAA)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,7260G→T 4.4% 25.6 ‑0.5 ‑0.5296H1237Q (CAC→CAA) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (8/5):  ref base (149/134):  total (157/139)
Fisher's exact test strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.39e-01

Only 100 of 331 total aligned reads displayed.

CGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCA  >  NC_007779/3451626‑3451826
                                                                                                    |                                                                                                    
cGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAgg                                                                                                      <  1:1899180/101‑1 (MQ=255)
cGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAgg                                                                                                      <  1:13051294/101‑1 (MQ=255)
  aagcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtg                                                                                                    <  1:2516382/98‑1 (MQ=255)
  aaaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtg                                                                                                    >  1:11313218/1‑101 (MQ=255)
   agcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtgg                                                                                                   >  1:1827591/3‑101 (MQ=255)
   aaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtgg                                                                                                   >  1:6658316/1‑101 (MQ=255)
    gcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtggt                                                                                                  <  1:6279260/100‑1 (MQ=255)
       cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                                                               >  1:10191905/1‑101 (MQ=255)
       cTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCa                                                                                               >  1:4966175/1‑101 (MQ=255)
       cTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCa                                                                                               <  1:11821085/101‑1 (MQ=255)
        tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                                              <  1:3661027/101‑1 (MQ=255)
        tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                                              >  1:14739163/1‑101 (MQ=255)
         gACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGt                                                                                             <  1:7076978/101‑1 (MQ=255)
          accaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGtt                                                                                            >  1:8586959/1‑101 (MQ=255)
          accaccGAACTGTGCCTTACCACCCAGCGGCTGCCGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGtt                                                                                            >  1:1274315/1‑101 (MQ=255)
          accaccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGGAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAgct                                                                                            >  1:11144337/1‑99 (MQ=255)
           ccaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCAGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGttt                                                                                           <  1:7029954/101‑1 (MQ=255)
            caccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTc                                                                                          >  1:217255/1‑101 (MQ=255)
            caccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTc                                                                                          >  1:4873505/1‑101 (MQ=255)
            caccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTc                                                                                          >  1:520775/1‑101 (MQ=255)
            caccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTc                                                                                          >  1:73352/1‑101 (MQ=255)
             accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                                                                         >  1:8320416/1‑101 (MQ=255)
               cGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGc                                                                                       <  1:7394985/101‑1 (MQ=255)
                gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                                      <  1:13156732/101‑1 (MQ=255)
                gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                                      >  1:8392202/1‑101 (MQ=255)
                gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                                      >  1:2810568/1‑101 (MQ=255)
                 aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                                                                     <  1:5052691/101‑1 (MQ=255)
                 aaCTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcat                                                                                     <  1:14133201/101‑1 (MQ=255)
                  aCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatg                                                                                    <  1:14824831/101‑1 (MQ=255)
                  aCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatg                                                                                    >  1:6204967/1‑101 (MQ=255)
                     gagCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtac                                                                                 >  1:6366345/3‑101 (MQ=255)
                       gCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacat                                                                               <  1:13649281/101‑1 (MQ=255)
                        ccTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatg                                                                              <  1:7827935/101‑1 (MQ=255)
                         cTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgt                                                                             >  1:8876052/1‑101 (MQ=255)
                          ttACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgta                                                                            >  1:6116082/1‑101 (MQ=255)
                          ttACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgta                                                                            <  1:9887495/101‑1 (MQ=255)
                          ttACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgta                                                                            <  1:14630916/101‑1 (MQ=255)
                           tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                                           >  1:13573798/1‑101 (MQ=255)
                           tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                                           <  1:15080073/101‑1 (MQ=255)
                            accaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAAc                                                                          >  1:14172669/1‑101 (MQ=255)
                             ccaccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGcc                                                                         <  1:12248223/101‑1 (MQ=255)
                              caccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccc                                                                        <  1:6482258/101‑2 (MQ=255)
                              caccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccc                                                                        >  1:7676872/1‑100 (MQ=255)
                                cccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGATCGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAAc                                                                      <  1:10884200/101‑1 (MQ=255)
                                 ccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACg                                                                     >  1:1285489/1‑101 (MQ=255)
                                    gCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaa                                                                  <  1:1892088/101‑5 (MQ=255)
                                     gggCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaa                                                                 <  1:15146871/100‑6 (MQ=255)
                                     cGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTa                                                                 <  1:1694293/101‑1 (MQ=255)
                                      ggCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaac                                                                <  1:13196794/101‑7 (MQ=255)
                                      ggCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAc                                                                >  1:11480398/1‑101 (MQ=255)
                                      ggCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAc                                                                >  1:5655321/1‑101 (MQ=255)
                                       gctgctGAGTAACCAGGCTGTAAGAACCGGTGGAGCGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAcc                                                               <  1:8093015/101‑1 (MQ=255)
                                         tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg                                                             >  1:1578919/1‑101 (MQ=255)
                                          gctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGa                                                            >  1:3193716/1‑101 (MQ=255)
                                          gctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGa                                                            <  1:981973/101‑1 (MQ=255)
                                           ctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGAc                                                           >  1:6391604/1‑101 (MQ=255)
                                           ctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGAc                                                           <  1:9378192/101‑1 (MQ=255)
                                            tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg                                                          >  1:10421470/1‑101 (MQ=255)
                                                 aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                                                     >  1:11665858/1‑101 (MQ=255)
                                                   ccAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAAc                                                   <  1:9473918/101‑1 (MQ=255)
                                                    cAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACt                                                  >  1:341333/1‑101 (MQ=255)
                                                     aGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTg                                                 <  1:721286/101‑1 (MQ=255)
                                                     aGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTg                                                 <  1:6032887/101‑1 (MQ=255)
                                                         tGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCa                                             >  1:1610701/1‑101 (MQ=255)
                                                           tAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCAcc                                           >  1:3089121/1‑101 (MQ=255)
                                                            aaGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCa                                          >  1:13699774/1‑101 (MQ=255)
                                                             aGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAg                                         <  1:5419598/101‑1 (MQ=255)
                                                              gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                                        >  1:12012881/1‑101 (MQ=255)
                                                              gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                                        >  1:3884852/1‑101 (MQ=255)
                                                              gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                                        >  1:4051189/1‑101 (MQ=255)
                                                                aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                      <  1:5094526/101‑1 (MQ=255)
                                                                aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                      <  1:8480102/101‑1 (MQ=255)
                                                                 ccGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCg                                     >  1:12555388/1‑101 (MQ=255)
                                                                 ccGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCg                                     >  1:9688375/1‑101 (MQ=255)
                                                                    gTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGAcc                                  <  1:12984496/101‑1 (MQ=255)
                                                                     tGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccg                                 >  1:9937658/1‑100 (MQ=255)
                                                                     tGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccg                                 >  1:4212530/1‑100 (MQ=255)
                                                                       gAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATc                               >  1:10770124/1‑101 (MQ=255)
                                                                         aCGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCga                             >  1:4130545/1‑100 (MQ=255)
                                                                           gcgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTAc                           >  1:12386711/1‑101 (MQ=255)
                                                                           gcgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTAc                           >  1:6961627/1‑101 (MQ=255)
                                                                            cgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACa                          >  1:874606/1‑101 (MQ=255)
                                                                            cgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACa                          >  1:10804504/1‑101 (MQ=255)
                                                                             gcgTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAg                         >  1:14518620/1‑101 (MQ=255)
                                                                              cgTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgc                        >  1:14895327/1‑100 (MQ=255)
                                                                               ttGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGc                       >  1:4599958/2‑101 (MQ=255)
                                                                               gTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGc                       <  1:9152658/101‑1 (MQ=255)
                                                                                tgaaTCTTCTCGTCGACCAGGTGGTGCATCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgcgt                      <  1:10246661/98‑4 (MQ=255)
                                                                                 gCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCgg                     <  1:813583/101‑1 (MQ=255)
                                                                                 gCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCgg                     <  1:441455/101‑1 (MQ=255)
                                                                                     cTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAgcgtggcct                 >  1:8529704/1‑93 (MQ=255)
                                                                                         tcgtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGAcc             <  1:2501232/101‑1 (MQ=255)
                                                                                          cgtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCg            >  1:4701792/1‑101 (MQ=255)
                                                                                          cgtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCg            >  1:7950709/1‑101 (MQ=255)
                                                                                           gtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCgg           >  1:10536403/1‑101 (MQ=255)
                                                                                            tcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGa          >  1:5868438/1‑101 (MQ=255)
                                                                                              gACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAg        <  1:2393874/101‑1 (MQ=255)
                                                                                                 cAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCg     <  1:4354433/101‑1 (MQ=255)
                                                                                                    gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa  <  1:1940705/101‑1 (MQ=255)
                                                                                                    gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa  >  1:13334516/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCA  >  NC_007779/3451626‑3451826

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