Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,449,991 G→A 8.1% I447I (ATC→ATT rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,449,9910G→A 8.1% 46.5 ‑0.7 ‑0.0257I447I (ATC→ATTrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (110/125):  total (118/139)
Fisher's exact test strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.59e-01

Only 100 of 272 total aligned reads displayed.

GCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTT  >  NC_007779/3449891‑3450091
                                                                                                    |                                                                                                    
gcgcTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTgg                                                                                                      >  1:5816198/1‑101 (MQ=255)
gcgcTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCAAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTgg                                                                                                      >  1:13695921/1‑101 (MQ=255)
 cgcTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGa                                                                                                     >  1:13371164/1‑101 (MQ=255)
 cgcTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGa                                                                                                     <  1:5608273/101‑1 (MQ=255)
 cgcTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGa                                                                                                     <  1:4669326/101‑1 (MQ=255)
   cTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATa                                                                                                   <  1:5087988/101‑1 (MQ=255)
   cTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATa                                                                                                   <  1:10808374/101‑1 (MQ=255)
     tCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGc                                                                                                 >  1:1893723/1‑101 (MQ=255)
      ccAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCCCCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCt                                                                                                <  1:12729477/101‑1 (MQ=255)
      ccAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCt                                                                                                >  1:244139/1‑101 (MQ=255)
       cAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCtt                                                                                               >  1:679361/1‑101 (MQ=255)
        aGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCttt                                                                                              >  1:9257788/1‑101 (MQ=255)
          cTGGGCTTCCAGCGTCAGCGGTGCGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTAc                                                                                            <  1:12101667/101‑1 (MQ=255)
          cTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTAc                                                                                            >  1:4857174/1‑101 (MQ=255)
          cTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTAc                                                                                            >  1:4192181/1‑101 (MQ=255)
               cTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCg                                                                                       <  1:4940918/101‑1 (MQ=255)
                ttCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGa                                                                                      >  1:10316808/1‑101 (MQ=255)
                 tCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCACCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCACCTGGATAGCTTTACCTTCGAt                                                                                     >  1:8954488/1‑101 (MQ=255)
                  ccAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATc                                                                                    <  1:1864039/101‑1 (MQ=255)
                   cAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTCCCGTCGATCa                                                                                   <  1:2830988/101‑1 (MQ=255)
                    aGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAg                                                                                  >  1:11972721/1‑101 (MQ=255)
                    aGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAg                                                                                  >  1:3642860/1‑101 (MQ=255)
                       gTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTAc                                                                               >  1:13892496/1‑101 (MQ=255)
                          aGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCgg                                                                            >  1:6412814/1‑101 (MQ=255)
                          aGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCgg                                                                            <  1:13499386/101‑1 (MQ=255)
                           gCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGt                                                                           >  1:11889890/1‑101 (MQ=255)
                            ccgTAAGTGAACAGCCATCTGGTCCCCATCGAAGTCGGCGTTATATGCCGCAGAAACCAGCGGGTGCAGCTGGATAGCTTTCCCTTCGATCAGTCCCGGtt                                                                          <  1:42144/99‑1 (MQ=255)
                            cGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGtt                                                                          >  1:6256353/1‑101 (MQ=255)
                            cGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGtt                                                                          <  1:6467438/101‑1 (MQ=255)
                            cGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGtt                                                                          >  1:3164741/1‑101 (MQ=255)
                             ggTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTc                                                                         <  1:10748277/101‑1 (MQ=255)
                              gTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCa                                                                        <  1:3750364/101‑1 (MQ=255)
                                     aaCAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCt                                                                 <  1:3406548/101‑1 (MQ=255)
                                     aaCAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCt                                                                 >  1:12810551/1‑101 (MQ=255)
                                       cAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTgg                                                               <  1:695088/101‑1 (MQ=255)
                                         gCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGAt                                                             >  1:12590713/1‑101 (MQ=255)
                                           cATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATAc                                                           <  1:1107121/101‑1 (MQ=255)
                                             tCTGTTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATAccc                                                         <  1:2069942/101‑1 (MQ=255)
                                             tCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATAccc                                                         >  1:2990101/1‑101 (MQ=255)
                                              cTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTTGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCa                                                        >  1:13935089/1‑101 (MQ=255)
                                              cTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCa                                                        >  1:7157771/1‑101 (MQ=255)
                                              cTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCa                                                        <  1:10391151/101‑1 (MQ=255)
                                               tGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAg                                                       >  1:13311091/1‑101 (MQ=255)
                                               tGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAg                                                       >  1:11264019/1‑101 (MQ=255)
                                                  tCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACg                                                    <  1:1509505/101‑1 (MQ=255)
                                                   cACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACgg                                                   >  1:5411251/1‑101 (MQ=255)
                                                    aCCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGt                                                  <  1:9344035/101‑1 (MQ=255)
                                                    aCCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGt                                                  >  1:11456424/1‑101 (MQ=255)
                                                    aCCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGt                                                  <  1:8086315/101‑1 (MQ=255)
                                                    aCCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGt                                                  >  1:7982824/1‑101 (MQ=255)
                                                     ccATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTg                                                 <  1:7564272/101‑1 (MQ=255)
                                                      cATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGc                                                <  1:9830563/101‑1 (MQ=255)
                                                       aTCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCa                                               >  1:1824039/1‑101 (MQ=255)
                                                       aTCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCa                                               >  1:3088556/1‑101 (MQ=255)
                                                       aTCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCa                                               <  1:6095901/101‑1 (MQ=255)
                                                           aaGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGt                                           <  1:4758152/101‑1 (MQ=255)
                                                            aGTCGGTGTTAGACGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGGTCCGGGTAAAATGCATAGATAAACAGACGGTGCAGAGTc                                          <  1:2112245/101‑1 (MQ=255)
                                                            aGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTc                                          <  1:4281865/101‑1 (MQ=255)
                                                             gTCGGCGTTATATGCCGCACAAACCAGCGGCTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGATACCCAGACCGTGCAGAGTCgt                                        >  1:14122538/1‑100 (MQ=255)
                                                              tCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCgg                                        <  1:271124/101‑1 (MQ=255)
                                                               cGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGt                                       <  1:6907474/101‑1 (MQ=255)
                                                                  cGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCa                                    <  1:38694/101‑1 (MQ=255)
                                                                  cGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCa                                    >  1:13549577/1‑101 (MQ=255)
                                                                      atatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGt                                >  1:7338217/1‑101 (MQ=255)
                                                                       tatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGtt                               <  1:2901419/101‑1 (MQ=255)
                                                                         tGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTca                             <  1:5885687/101‑1 (MQ=255)
                                                                           cctcacatgccAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTcagc                           <  1:8145173/92‑1 (MQ=255)
                                                                           ccGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagc                           >  1:11790261/1‑101 (MQ=255)
                                                                            cGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagca                          >  1:3495715/1‑101 (MQ=255)
                                                                              cacaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGt                        <  1:1975341/101‑1 (MQ=255)
                                                                               gctatccAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                       >  1:989531/6‑101 (MQ=255)
                                                                               acaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                       >  1:1650133/1‑101 (MQ=255)
                                                                                 aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                     <  1:4997267/101‑1 (MQ=255)
                                                                                 aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                     >  1:4084234/1‑101 (MQ=255)
                                                                                 aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                     <  1:8358056/101‑1 (MQ=255)
                                                                                 aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                     <  1:8186866/101‑1 (MQ=255)
                                                                                 aaaCAAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                     >  1:3907486/1‑101 (MQ=255)
                                                                                     cAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGgt                 <  1:10601924/101‑1 (MQ=255)
                                                                                     cAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCCAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCgggg                 >  1:12247308/1‑100 (MQ=255)
                                                                                     cAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCCAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGgt                 >  1:10898976/1‑101 (MQ=255)
                                                                                       gCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgtgt               >  1:770625/1‑101 (MQ=255)
                                                                                        cGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGtt              <  1:14056222/101‑1 (MQ=255)
                                                                                         gggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTc             >  1:86551/1‑101 (MQ=255)
                                                                                          ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcg            >  1:10254218/1‑101 (MQ=255)
                                                                                          ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcg            <  1:2172206/101‑1 (MQ=255)
                                                                                          ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcg            <  1:4573416/101‑1 (MQ=255)
                                                                                           gTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcgc           >  1:1435224/1‑101 (MQ=255)
                                                                                           gTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATGCCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGGTGTCcgg           <  1:14159685/101‑2 (MQ=255)
                                                                                             gCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGGTGTCcggta         <  1:13054887/101‑4 (MQ=255)
                                                                                              cAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGa        <  1:7509390/101‑1 (MQ=255)
                                                                                              cAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGa        <  1:9633913/101‑1 (MQ=255)
                                                                                               aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGAt       >  1:11706587/1‑101 (MQ=255)
                                                                                               aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGAt       >  1:11490872/1‑101 (MQ=255)
                                                                                               aGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGGTGTCcggtaat       >  1:8529042/1‑96 (MQ=255)
                                                                                                gCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATa      <  1:7068719/101‑1 (MQ=255)
                                                                                                gCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCCGAATACCCAGGCGGTGCAGTGTAGGTGCACGGTTCAGCAGGACCGGGGGTCcggtaatg      >  1:11859196/1‑95 (MQ=255)
                                                                                                 cTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATaa     >  1:12517571/1‑101 (MQ=255)
                                                                                                 cTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATaa     >  1:1234390/1‑101 (MQ=255)
                                                                                                   ggATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACt   >  1:12700765/1‑101 (MQ=255)
                                                                                                    gATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACtt  <  1:261335/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTT  >  NC_007779/3449891‑3450091

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