Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,773 C→G 6.1% E1222Q (GAG→CAG)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,7730C→G 6.1% 34.2 ‑1.8 ‑0.1231E1222Q (GAG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (10/4):  ref base (106/111):  total (116/115)
Fisher's exact test strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.67e-01

Only 100 of 241 total aligned reads displayed.

TAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCACCG  >  NC_007779/3451674‑3451873
                                                                                                   |                                                                                                    
tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                                                                                                     >  1:1277053/1‑101 (MQ=255)
tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                                                                                                     >  1:3896560/1‑101 (MQ=255)
 aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                                                                                                    >  1:9326629/1‑101 (MQ=255)
   ccAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAAc                                                                                                  >  1:4798672/1‑101 (MQ=255)
     aGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTg                                                                                                >  1:4107294/1‑101 (MQ=255)
        cTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTc                                                                                             <  1:3035089/101‑1 (MQ=255)
        cTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTc                                                                                             <  1:11410735/101‑1 (MQ=255)
          gTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGTACGCTCGAACTGTTCAc                                                                                           >  1:5419256/1‑101 (MQ=255)
           tAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCAcc                                                                                          >  1:13330123/1‑101 (MQ=255)
            aaGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCa                                                                                         <  1:7852776/101‑1 (MQ=255)
            aaGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCa                                                                                         <  1:11723653/101‑1 (MQ=255)
              gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                                                                                       >  1:9636631/1‑101 (MQ=255)
                aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                                                                     <  1:11108021/101‑1 (MQ=255)
                aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                                                                     <  1:13948490/101‑1 (MQ=255)
                aCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGc                                                                                     >  1:12832610/1‑101 (MQ=255)
                 ccGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCg                                                                                    >  1:5559358/1‑101 (MQ=255)
                  cGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGa                                                                                   <  1:13373596/101‑1 (MQ=255)
                  cGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGa                                                                                   <  1:4529096/101‑1 (MQ=255)
                  cGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGa                                                                                   >  1:13590075/1‑101 (MQ=255)
                     tGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCa                                                                                >  1:5168308/1‑101 (MQ=255)
                       gAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTc                                                                              <  1:13696643/101‑1 (MQ=255)
                        aaCGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCg                                                                             <  1:280298/101‑1 (MQ=255)
                          cgcgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTa                                                                           >  1:6598426/1‑101 (MQ=255)
                            cgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACa                                                                         >  1:1926129/1‑101 (MQ=255)
                            agcgTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGCCCGTCGAACa                                                                         <  1:13874133/100‑1 (MQ=255)
                             gcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAg                                                                        >  1:8146066/1‑101 (MQ=255)
                             gcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAg                                                                        >  1:4722635/1‑101 (MQ=255)
                              cgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAgg                                                                       >  1:1413959/1‑101 (MQ=255)
                               gTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGc                                                                      >  1:11644876/1‑101 (MQ=255)
                               gTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGc                                                                      <  1:6387330/101‑1 (MQ=255)
                                tGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCg                                                                     <  1:6920682/101‑1 (MQ=255)
                                 gCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCgg                                                                    <  1:3558395/101‑1 (MQ=255)
                                  cATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGa                                                                   <  1:11654569/101‑1 (MQ=255)
                                   aTCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGTGACCATCGTACAGGCGGAt                                                                  >  1:2707213/1‑101 (MQ=255)
                                   aTCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGAt                                                                  >  1:1212251/1‑101 (MQ=255)
                                    tCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATc                                                                 <  1:11159980/101‑1 (MQ=255)
                                      ttGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTg                                                               <  1:6145525/101‑1 (MQ=255)
                                       tGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGa                                                              >  1:4075946/1‑101 (MQ=255)
                                         tcgtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGAcc                                                            <  1:7250471/101‑1 (MQ=255)
                                          cgtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCg                                                           >  1:11957309/1‑101 (MQ=255)
                                           gtcgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCgg                                                          >  1:8937271/1‑101 (MQ=255)
                                             cgACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGaa                                                        >  1:1398417/1‑101 (MQ=255)
                                               aCCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGt                                                      >  1:3975224/1‑101 (MQ=255)
                                                  aGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCgg                                                   >  1:11225595/1‑101 (MQ=255)
                                                  aGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCgg                                                   <  1:10199310/101‑1 (MQ=255)
                                                    ttggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGCCCGGAAGTCGGCa                                                 <  1:4385268/100‑1 (MQ=255)
                                                    gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa                                                 <  1:5262600/101‑1 (MQ=255)
                                                    gtggtTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCa                                                 <  1:10883591/101‑1 (MQ=255)
                                                        ttCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTc                                             <  1:2797717/101‑1 (MQ=255)
                                                        ttCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTc                                             >  1:3316833/1‑101 (MQ=255)
                                                        ttCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTc                                             <  1:7131524/101‑1 (MQ=255)
                                                         tCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCg                                            >  1:1983426/1‑101 (MQ=255)
                                                          cAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGc                                           <  1:8729205/101‑1 (MQ=255)
                                                            gTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGCCCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCa                                         <  1:848844/101‑1 (MQ=255)
                                                              ttCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAg                                       >  1:6089095/1‑101 (MQ=255)
                                                              ttCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAg                                       <  1:11322704/101‑1 (MQ=255)
                                                               tCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGt                                      <  1:3217199/101‑1 (MQ=255)
                                                                   catgtacatgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTca                                  <  1:8467457/101‑1 (MQ=255)
                                                                    atgtacatgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTcag                                 <  1:12519803/101‑1 (MQ=255)
                                                                        acatgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCagcagc                             >  1:10094040/1‑101 (MQ=255)
                                                                        acatgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCagcagc                             <  1:13142859/101‑1 (MQ=255)
                                                                        acatgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCagcagc                             >  1:13493628/1‑101 (MQ=255)
                                                                         catgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGct                            >  1:5632038/1‑101 (MQ=255)
                                                                         catgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGct                            <  1:7073055/101‑1 (MQ=255)
                                                                          atgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGctc                           >  1:11208329/1‑101 (MQ=255)
                                                                          atgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGctc                           >  1:5401282/1‑101 (MQ=255)
                                                                          atgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGGCGCCAAGTTTCAGCAGctc                           >  1:2663514/1‑101 (MQ=255)
                                                                           tgtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGctct                          <  1:7156468/101‑1 (MQ=255)
                                                                            gtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCtt                         <  1:7288873/101‑1 (MQ=255)
                                                                            gtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCtt                         <  1:2508713/101‑1 (MQ=255)
                                                                            gtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCtt                         <  1:3442475/101‑1 (MQ=255)
                                                                            gtaACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCACCAAGTTTCAGCAGCTCtt                         <  1:7685196/101‑1 (MQ=255)
                                                                              aaCCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTa                       >  1:4464902/1‑101 (MQ=255)
                                                                                 cAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAAtt                    <  1:13761348/101‑1 (MQ=255)
                                                                                 cAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAAtt                    >  1:6710127/1‑101 (MQ=255)
                                                                                  aaCGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAAttt                   >  1:4589099/1‑101 (MQ=255)
                                                                                  aaCGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAAttt                   <  1:1539638/101‑1 (MQ=255)
                                                                                  aaCGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTCACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAAttt                   >  1:9933694/1‑101 (MQ=255)
                                                                                    cGGCTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCTTACAGGCGGATCTGACCGGAAGTCGGCAGGTGGCCAAGTTTCAGCAGCTCTTTAATTTCt                 >  1:11384228/1‑101 (MQ=255)
                                                                                     ggTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTg                <  1:5858160/101‑1 (MQ=255)
                                                                                       ttACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGct              >  1:10506649/1‑101 (MQ=255)
                                                                                        tACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGctt             >  1:6458749/1‑101 (MQ=255)
                                                                                         aCCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGcttc            <  1:5220548/101‑1 (MQ=255)
                                                                                          ccGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGcttct           <  1:1838499/101‑1 (MQ=255)
                                                                                          ccGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGcttct           <  1:3697134/101‑1 (MQ=255)
                                                                                           cGGACGCTCGAACTGTTCACCAGTGCGCCCATCGTGCAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGcttctt          <  1:2117516/101‑1 (MQ=255)
                                                                                           cGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGcttctt          <  1:5521333/101‑1 (MQ=255)
                                                                                              aCGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGACGTCGGCTGGTCGCCCAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCg       >  1:3951575/1‑101 (MQ=255)
                                                                                              aCGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCg       >  1:10688879/1‑101 (MQ=255)
                                                                                              aCGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCg       <  1:10229102/101‑1 (MQ=255)
                                                                                               cGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGc      <  1:6247879/101‑1 (MQ=255)
                                                                                                gCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCa     >  1:3979759/1‑101 (MQ=255)
                                                                                                gCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCa     <  1:10725643/101‑1 (MQ=255)
                                                                                                 cTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCAc    >  1:13555406/1‑101 (MQ=255)
                                                                                                  tCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCAcc   >  1:12263838/1‑101 (MQ=255)
                                                                                                  tCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCAcc   <  1:4244417/101‑1 (MQ=255)
                                                                                                  tCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTACTTCTTTCGCAcc   >  1:2159728/1‑101 (MQ=255)
                                                                                                   cGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCACCg  >  1:9054807/1‑101 (MQ=255)
                                                                                                   cGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCACCg  >  1:10419489/1‑101 (MQ=255)
                                                                                                   cGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCACCg  <  1:9393457/101‑1 (MQ=255)
                                                                                                   |                                                                                                    
TAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGGATCTGACCGGAAGTCGGCAGGTCGCCAAGTTTCAGCAGCTCTTTAATTTCTGCTTCTTTCGCACCG  >  NC_007779/3451674‑3451873

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