Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,486,012 C→G 3.6% S394C (TCC→TGC)  hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,486,0120C→G 3.6% 14.6 ‑0.7 ‑0.0223S394C (TCC→TGC) hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (3/5):  ref base (113/102):  total (116/107)
Fisher's exact test strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.84e-01

Only 100 of 232 total aligned reads displayed.

GTGCAGCGTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCGCTGGGGCTGGGC  >  NC_007779/1485912‑1486112
                                                                                                    |                                                                                                    
 tGCAGCGTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTcc                                                                                                     <  1:13236693/101‑1 (MQ=255)
   cAGCGTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGa                                                                                                   <  1:2223146/101‑1 (MQ=255)
     gCGTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAg                                                                                                 <  1:7406796/101‑1 (MQ=255)
       gTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGAc                                                                                               <  1:1859774/101‑1 (MQ=255)
       gTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGAc                                                                                               <  1:7008561/101‑1 (MQ=255)
        tttGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACg                                                                                              <  1:264485/101‑1 (MQ=255)
        tttGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACg                                                                                              <  1:10045718/101‑1 (MQ=255)
          tGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGAt                                                                                            <  1:8811549/101‑1 (MQ=255)
              gATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTcc                                                                                        >  1:9470964/1‑101 (MQ=255)
               aTTGAGCCGATTTCCCATGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCct                                                                                       <  1:14867784/101‑1 (MQ=255)
               aTTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCct                                                                                       <  1:10271617/101‑1 (MQ=255)
               aTTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCct                                                                                       <  1:1239049/101‑1 (MQ=255)
                ttGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCctc                                                                                      >  1:14085801/1‑101 (MQ=255)
                  gAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCtctc                                                                                    <  1:545143/101‑1 (MQ=255)
                   aGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTcg                                                                                   >  1:8906298/1‑101 (MQ=255)
                    gCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTcgc                                                                                  >  1:8224943/1‑101 (MQ=255)
                    gCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTcgc                                                                                  >  1:10936404/1‑101 (MQ=255)
                     ccGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTcgcg                                                                                 >  1:6619745/1‑101 (MQ=255)
                     ccGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTcgcg                                                                                 <  1:939869/101‑1 (MQ=255)
                       gATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGcc                                                                               >  1:6398185/1‑101 (MQ=255)
                       gATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGcc                                                                               >  1:6072734/1‑101 (MQ=255)
                         tttCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCg                                                                             <  1:4022796/101‑1 (MQ=255)
                            cccAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAg                                                                          <  1:6816501/101‑1 (MQ=255)
                             ccAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGt                                                                         >  1:3279439/1‑101 (MQ=255)
                              cAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGtt                                                                        >  1:1179902/1‑101 (MQ=255)
                                  cGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGAGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTAcc                                                                    >  1:12071385/1‑101 (MQ=255)
                                    tCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGa                                                                  >  1:734181/1‑101 (MQ=255)
                                    tCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGa                                                                  <  1:13638294/101‑1 (MQ=255)
                                     cTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGAt                                                                 >  1:4183196/1‑101 (MQ=255)
                                       gCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATcc                                                               <  1:10202336/101‑1 (MQ=255)
                                        ccAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCg                                                              <  1:5713885/101‑1 (MQ=255)
                                          aaTCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGga                                                            <  1:14704961/101‑1 (MQ=255)
                                             cAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGAt                                                         >  1:567533/1‑101 (MQ=255)
                                              aGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGAtt                                                        >  1:9838933/1‑101 (MQ=255)
                                                 gTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTg                                                     >  1:6458858/1‑101 (MQ=255)
                                                 gTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTg                                                     >  1:2080042/1‑101 (MQ=255)
                                                 gTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTg                                                     <  1:1814652/101‑1 (MQ=255)
                                                   aaaGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCg                                                   >  1:11357203/1‑101 (MQ=255)
                                                   aaaGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCg                                                   >  1:7242961/1‑101 (MQ=255)
                                                    aaGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGt                                                  <  1:12691204/101‑1 (MQ=255)
                                                     aGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTa                                                 <  1:13849030/101‑1 (MQ=255)
                                                     aGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTa                                                 >  1:4684664/1‑101 (MQ=255)
                                                        ccGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCa                                              >  1:5099143/1‑101 (MQ=255)
                                                         cGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCaa                                             <  1:14068448/101‑1 (MQ=255)
                                                         cGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCaa                                             <  1:13271447/101‑1 (MQ=255)
                                                          gCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAAc                                            >  1:8770604/1‑101 (MQ=255)
                                                          gCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAAc                                            <  1:9038042/101‑1 (MQ=255)
                                                          gCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAAc                                            >  1:271670/1‑101 (MQ=255)
                                                          gCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAAc                                            <  1:13542241/101‑1 (MQ=255)
                                                           cTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAAcc                                           <  1:12884387/101‑1 (MQ=255)
                                                           cTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAAcc                                           <  1:12299945/101‑1 (MQ=255)
                                                              tgGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTgg                                        <  1:7830790/101‑1 (MQ=255)
                                                               ggTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGc                                       >  1:9249346/1‑101 (MQ=255)
                                                                gTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGcc                                      >  1:10405676/1‑101 (MQ=255)
                                                                  tgtgtgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTc                                    <  1:9655597/100‑1 (MQ=255)
                                                                  cGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTc                                    >  1:4463496/1‑101 (MQ=255)
                                                                   gtgtgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCg                                   <  1:12254096/101‑1 (MQ=255)
                                                                   gtgtgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCg                                   >  1:84948/1‑101 (MQ=255)
                                                                    tgtgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCgg                                  >  1:14219625/1‑101 (MQ=255)
                                                                    tgtgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCgg                                  >  1:10648583/1‑101 (MQ=255)
                                                                     gtgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGt                                 <  1:12974631/101‑1 (MQ=255)
                                                                      tgtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGtt                                >  1:160632/1‑101 (MQ=255)
                                                                       gtCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTa                               >  1:14989973/1‑101 (MQ=255)
                                                                        tCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTAt                              >  1:7621732/1‑101 (MQ=255)
                                                                        tCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTAt                              >  1:4353382/1‑101 (MQ=255)
                                                                          cGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTAttt                            <  1:809191/101‑1 (MQ=255)
                                                                           gAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTAtttt                           >  1:6238112/1‑101 (MQ=255)
                                                                           gAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTAtttt                           <  1:3389544/101‑1 (MQ=255)
                                                                           gAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTAtttt                           <  1:3548202/101‑1 (MQ=255)
                                                                            aaGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTg                          <  1:10225789/101‑1 (MQ=255)
                                                                            aaGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTg                          <  1:14557775/101‑1 (MQ=255)
                                                                              gggATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCa                        <  1:11542621/101‑1 (MQ=255)
                                                                               ggATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAg                       >  1:11595016/1‑101 (MQ=255)
                                                                                gATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGa                      >  1:12449820/1‑101 (MQ=255)
                                                                                gATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGa                      <  1:884240/101‑1 (MQ=255)
                                                                                 aTCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGAt                     <  1:9162319/101‑1 (MQ=255)
                                                                                 aTCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGAt                     >  1:7918641/1‑101 (MQ=255)
                                                                                  tCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATg                    >  1:13985614/1‑101 (MQ=255)
                                                                                  tCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATg                    <  1:6128838/101‑1 (MQ=255)
                                                                                   cTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGa                   <  1:150353/101‑1 (MQ=255)
                                                                                   cTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGa                   >  1:7976323/1‑101 (MQ=255)
                                                                                   cTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGa                   >  1:3060633/1‑101 (MQ=255)
                                                                                     gTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGAcc                 <  1:1273533/101‑1 (MQ=255)
                                                                                      tATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACcg                <  1:7739120/101‑1 (MQ=255)
                                                                                      tATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACcg                >  1:414663/1‑101 (MQ=255)
                                                                                      gaTTCGTCTCTACAGCGAGGAGGAATTTCTCTCTCGCCCCGAGTTTACGGATCCGGAGATTCTGCGCACCAATCTCGCTTCGGTTGTACTGCAAATGattt                <  1:2562552/100‑4 (MQ=255)
                                                                                        ttCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACcgcg              <  1:3840792/101‑1 (MQ=255)
                                                                                        ttCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACcgcg              <  1:4907307/101‑1 (MQ=255)
                                                                                         tCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCgcgc             <  1:5895908/101‑1 (MQ=255)
                                                                                          cGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgct            <  1:13991315/101‑1 (MQ=255)
                                                                                           gTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctg           <  1:10730710/101‑1 (MQ=255)
                                                                                           gTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCCTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctg           <  1:11972910/101‑1 (MQ=255)
                                                                                            tctcTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctgg          >  1:7419424/1‑101 (MQ=255)
                                                                                              tctATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctgggg        >  1:9488161/1‑101 (MQ=255)
                                                                                                 aTTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctggggctg     >  1:11145654/1‑101 (MQ=255)
                                                                                                  ttCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctggggctgg    <  1:5629388/101‑1 (MQ=255)
                                                                                                  ttCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctggggctgg    >  1:3439515/1‑101 (MQ=255)
                                                                                                   tCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctggggctggg   >  1:8323997/1‑101 (MQ=255)
                                                                                                   tCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctggggctggg   <  1:1699491/101‑1 (MQ=255)
                                                                                                   tCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCgctggggctggg   >  1:9939846/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GTGCAGCGTTTGCCGATTGAGCCGATTTCCCAGGCGTCTGCCAATCAGCGTAAAGGCCGCTGTGGTCGTGTGTCCGAAGGGATCTGTATTCGTCTCTATTCCGAAGACGATTTCCTCTCGCGCCCGGAGTTTACCGATCCGGAGATTCTGCGTACCAACCTGGCCTCGGTTATTTTGCAGATGACCGCGCTGGGGCTGGGC  >  NC_007779/1485912‑1486112

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