Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,453,849 T→G 6.3% I530L (ATC→CTC)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,453,8490T→G 6.3% 41.9 ‑0.2 ‑0.9223I530L (ATC→CTC) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (8/6):  ref base (108/101):  total (116/107)
Fisher's exact test strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.20e-01

Only 100 of 240 total aligned reads displayed.

GTTTCGATTGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATCGGC  >  NC_007779/3453749‑3453949
                                                                                                    |                                                                                                    
gTTTCGATTGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGAt                                                                                                      <  1:13280995/101‑1 (MQ=255)
 tttCGATTGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATa                                                                                                     <  1:1811491/101‑1 (MQ=255)
  ttCGATTGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATac                                                                                                    >  1:11899227/1‑101 (MQ=255)
     gATTGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATacgac                                                                                                 >  1:13670792/1‑101 (MQ=255)
       ttGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGt                                                                                               >  1:4756774/1‑101 (MQ=255)
        ttgACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGtt                                                                                              <  1:2030358/99‑1 (MQ=255)
        tGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACTACGtt                                                                                              >  1:8473056/1‑101 (MQ=255)
        tGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGtt                                                                                              <  1:10073692/101‑1 (MQ=255)
          gACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTtg                                                                                            >  1:11460272/1‑101 (MQ=255)
           aCATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTtgt                                                                                           >  1:21173/1‑101 (MQ=255)
             aTACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGc                                                                                         >  1:6506025/1‑101 (MQ=255)
              tACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCg                                                                                        >  1:6795611/1‑101 (MQ=255)
              tACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCg                                                                                        <  1:691010/101‑1 (MQ=255)
               aCGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGt                                                                                       <  1:13345171/101‑1 (MQ=255)
                cgcgACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTa                                                                                      >  1:12685074/1‑101 (MQ=255)
                 gcgACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCCAGTGCGGAGATACGACGTTTGTGCGTaa                                                                                     <  1:12964265/101‑1 (MQ=255)
                  cgACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAAt                                                                                    >  1:3580172/1‑101 (MQ=255)
                   gACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAAtc                                                                                   >  1:11107167/1‑101 (MQ=255)
                   gACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAAtc                                                                                   >  1:672062/1‑101 (MQ=255)
                     ccGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAAtctc                                                                                 <  1:12635279/101‑1 (MQ=255)
                      ggTAGTGAGTGGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCa                                                                                <  1:7222988/100‑1 (MQ=255)
                       gTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAg                                                                               >  1:8917779/1‑101 (MQ=255)
                         aGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGAc                                                                             <  1:5285638/101‑1 (MQ=255)
                           tGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAg                                                                           >  1:9271213/1‑101 (MQ=255)
                            gAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGc                                                                          <  1:8156940/101‑1 (MQ=255)
                             aGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCg                                                                         <  1:855621/101‑1 (MQ=255)
                              gTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCgg                                                                        <  1:3419904/101‑1 (MQ=255)
                               tCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCggg                                                                       >  1:12038197/1‑101 (MQ=255)
                                cGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGgt                                                                      <  1:4362118/101‑1 (MQ=255)
                                    tgtACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGgttgtt                                                                  <  1:12148715/101‑1 (MQ=255)
                                        cGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTgg                                                              >  1:9464483/1‑101 (MQ=255)
                                         gTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGt                                                             <  1:9234419/101‑1 (MQ=255)
                                          tctcGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCGGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTc                                                            >  1:14250174/1‑101 (MQ=255)
                                           ctcGAACTTCGAAGCCTGCACGTTCACGGGTCAGGCCGCCCGGGCCTTGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTcc                                                           <  1:9320018/101‑1 (MQ=255)
                                           ctcGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTcc                                                           <  1:7104813/101‑1 (MQ=255)
                                            tcGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCa                                                          >  1:14071025/1‑101 (MQ=255)
                                             cGAACTTCGAAGCCTGCCCGTTCCCGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCAt                                                         <  1:1642626/101‑1 (MQ=255)
                                             cGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCAt                                                         <  1:1848717/101‑1 (MQ=255)
                                             cGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCAt                                                         >  1:10265681/1‑101 (MQ=255)
                                              gAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATa                                                        >  1:6987485/1‑101 (MQ=255)
                                              gAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATa                                                        >  1:11343704/1‑101 (MQ=255)
                                               aaCTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATaa                                                       >  1:1088474/1‑101 (MQ=255)
                                               aaCTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATaa                                                       >  1:2852404/1‑101 (MQ=255)
                                                aCTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATaaa                                                      <  1:11570466/101‑1 (MQ=255)
                                                      agaccagCACGGTCACGGGTCAGACCACCCGGACCCAGTGCGGACAGACGACGCTTGTGGGTCAGACCCGACAGCGAGTTGTTCTGGTCCATGAACTggga                                                <  1:1062351/95‑1 (MQ=255)
                                                      agaccagCACGGTCACGGGTCAGACCACCCGGACCCAGTGCGGACAGACGACGCTTGTGGGTCAGACCCGACAGCGAGTTGTTCTGGTCCATGAACTggga                                                <  1:13632926/95‑1 (MQ=255)
                                                      aaGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTgaga                                                >  1:10100013/1‑101 (MQ=255)
                                                      aaGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTgaga                                                <  1:12953151/101‑1 (MQ=255)
                                                      aaGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTgaga                                                >  1:13010790/1‑101 (MQ=255)
                                                         ccTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAg                                             <  1:14295783/101‑1 (MQ=255)
                                                          cTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAgc                                            <  1:2206887/101‑1 (MQ=255)
                                                           tGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAgct                                           >  1:4991455/1‑101 (MQ=255)
                                                            gCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGTAGATACGACGTTTGTGCGTAATCTCAGACAGCTGGTTGTTCTGGTCCATAAACTGAGACAgctg                                          <  1:4811072/101‑1 (MQ=255)
                                                            gCACGGTCACGGGTCAGACCACCCGGACCCAGTGCGGACAGACGACGCTTGTGGGTCAGACCCGACAGCGAGTTGTTCTGGTCCATGAACTGGGACAgctg                                          <  1:14111075/101‑1 (MQ=255)
                                                             cACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAgctgg                                         >  1:13853834/1‑101 (MQ=255)
                                                             cACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAgctgg                                         <  1:7778035/101‑1 (MQ=255)
                                                               cGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAgctggct                                       >  1:1002107/1‑101 (MQ=255)
                                                                 ttCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGctggctgg                                     <  1:10735336/101‑1 (MQ=255)
                                                                 ttCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGctggctgg                                     <  1:3599542/101‑1 (MQ=255)
                                                                   cACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGaa                                   <  1:13535131/101‑1 (MQ=255)
                                                                     cGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAAcc                                 <  1:2281752/101‑1 (MQ=255)
                                                                     cGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAAcc                                 <  1:2733129/101‑1 (MQ=255)
                                                                     cGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAAcc                                 <  1:5098215/101‑1 (MQ=255)
                                                                      gggTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCg                                >  1:2296510/1‑101 (MQ=255)
                                                                        gTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCgaa                              <  1:6132271/101‑1 (MQ=255)
                                                                          cAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCgaaga                            >  1:7891785/1‑101 (MQ=255)
                                                                          cAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCgaaga                            >  1:4554735/1‑101 (MQ=255)
                                                                            gACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAAc                          >  1:7460257/1‑101 (MQ=255)
                                                                             aCCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAAct                         >  1:12557202/1‑101 (MQ=255)
                                                                              ccGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAActc                        >  1:5515234/1‑101 (MQ=255)
                                                                              ccGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAActc                        <  1:2603274/101‑1 (MQ=255)
                                                                              ccGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAActc                        >  1:1944703/1‑101 (MQ=255)
                                                                               cGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAActct                       >  1:5444144/1‑101 (MQ=255)
                                                                               cGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAActct                       >  1:6099242/1‑101 (MQ=255)
                                                                                gCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCtt                      <  1:4371693/101‑1 (MQ=255)
                                                                                gCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCtt                      <  1:12216940/101‑1 (MQ=255)
                                                                                gCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCtt                      <  1:7592818/101‑1 (MQ=255)
                                                                                 ccTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCttt                     >  1:4486393/1‑101 (MQ=255)
                                                                                   tGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCa                   >  1:8137266/1‑101 (MQ=255)
                                                                                   tGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCa                   >  1:8799654/1‑101 (MQ=255)
                                                                                     ggCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCAct                 <  1:6436970/101‑1 (MQ=255)
                                                                                     ggCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCAct                 >  1:4273731/1‑101 (MQ=255)
                                                                                      gCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCActg                >  1:5949269/1‑101 (MQ=255)
                                                                                        cGAGTGCGGGGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTTActgct              <  1:3039357/101‑1 (MQ=255)
                                                                                        cGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCActgct              >  1:9434649/1‑101 (MQ=255)
                                                                                        cGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCActgct              <  1:9556340/101‑1 (MQ=255)
                                                                                         gAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCActgctg             <  1:13930589/101‑1 (MQ=255)
                                                                                         caGTGCGGAGATATGCCGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGCCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCGctgctg             <  1:2325550/100‑1 (MQ=255)
                                                                                          aGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACtgctgc            <  1:10978834/101‑1 (MQ=255)
                                                                                           gTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCg           <  1:5555574/101‑1 (MQ=255)
                                                                                            tGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCgg          >  1:14153643/1‑101 (MQ=255)
                                                                                            tGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCgg          >  1:5833264/1‑101 (MQ=255)
                                                                                            tGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCgg          <  1:9088417/101‑1 (MQ=255)
                                                                                              cGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGaa        >  1:10375707/1‑101 (MQ=255)
                                                                                                gagaTACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAAt      >  1:8244483/1‑101 (MQ=255)
                                                                                                 agaTACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATc     >  1:973454/1‑101 (MQ=255)
                                                                                                 agaTACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATc     <  1:5206323/101‑1 (MQ=255)
                                                                                                   aTACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATCgg   >  1:3252562/1‑101 (MQ=255)
                                                                                                   aTACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATCgg   >  1:3170167/1‑101 (MQ=255)
                                                                                                    tACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATCGGc  >  1:3789935/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GTTTCGATTGGACATACGCGACCGTAGTGAGTCGGGTGTACGTCTCGAACTTCGAAGCCTGCACGTTCACGGGTCAGACCGCCTGGGCCGAGTGCGGAGATACGACGTTTGTGCGTAATCTCAGACAGCGGGTTGTTCTGGTCCATAAACTGAGACAGCTGGCTGGAACCGAAGAACTCTTTCACTGCTGCGGAAATCGGC  >  NC_007779/3453749‑3453949

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.