Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,332 C→A 14.9% R167R (CGC→CGA)  hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,3320C→A14.9% 96.5 ‑3.9 ‑0.4255R167R (CGC→CGA) hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (13/25):  ref base (121/96):  total (134/121)
Fisher's exact test strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.23e-01

Only 100 of 259 total aligned reads displayed.

ATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGC  >  NC_007779/1485234‑1485432
                                                                                                  |                                                                                                      
aTCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGAATCCTGCTGGCGGAGATCCAGCAAGACCGCct                                                                                                      >  1:13021008/1‑101 (MQ=255)
 tCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctg                                                                                                     >  1:9301524/1‑101 (MQ=255)
   ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                                   <  1:9090210/101‑1 (MQ=255)
   ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                                   <  1:4660275/101‑1 (MQ=255)
       aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑Tg                                                                                             <  1:7421544/101‑1 (MQ=255)
        cAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGc                                                                                            <  1:7750607/101‑1 (MQ=255)
         aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCa                                                                                           >  1:10236725/1‑101 (MQ=255)
         aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCa                                                                                           >  1:7642363/1‑101 (MQ=255)
          aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAg                                                                                          <  1:1522024/101‑1 (MQ=255)
          aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAg                                                                                          <  1:5197667/101‑1 (MQ=255)
             tGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTAc                                                                                       <  1:13622458/101‑1 (MQ=255)
               cGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGa                                                                                     >  1:100778/1‑101 (MQ=255)
                gtgtCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGac                                                                                    <  1:11955556/98‑1 (MQ=255)
                gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGTTCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGac                                                                                    <  1:8906981/101‑1 (MQ=255)
                gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGac                                                                                    <  1:3487222/101‑1 (MQ=255)
                 tttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGaca                                                                                   >  1:1401726/1‑101 (MQ=255)
                  ttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGacac                                                                                  >  1:3183072/1‑101 (MQ=255)
                   tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACt                                                                                 >  1:13499961/1‑101 (MQ=255)
                   tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACt                                                                                 >  1:9846945/1‑101 (MQ=255)
                      gCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATc                                                                              >  1:6380570/1‑101 (MQ=255)
                        gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCat                                                                            >  1:4884349/1‑101 (MQ=255)
                        gATCACGTAAGTGATAACACGATGGTCAAGCTGACGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCat                                                                            >  1:395564/1‑101 (MQ=255)
                          tCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCatta                                                                          >  1:2575538/1‑101 (MQ=255)
                           cACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCattat                                                                         >  1:11886769/1‑101 (MQ=255)
                            aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCattatt                                                                        >  1:14056225/1‑101 (MQ=255)
                            aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCattatt                                                                        <  1:13276779/101‑1 (MQ=255)
                            aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCattatt                                                                        >  1:8634934/1‑101 (MQ=255)
                            aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCattatt                                                                        <  1:9792955/101‑1 (MQ=255)
                               tAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGAc                                                                     >  1:425465/1‑101 (MQ=255)
                                aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACg                                                                    >  1:7726089/1‑101 (MQ=255)
                                 aGTGATAACACGATGGTCAAGCAGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGa                                                                   >  1:12329696/1‑101 (MQ=255)
                                    gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGTAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAgc                                                                >  1:11746237/1‑101 (MQ=255)
                                     aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAgcg                                                               <  1:4046499/101‑1 (MQ=255)
                                     aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAgcg                                                               <  1:6094406/101‑1 (MQ=255)
                                       aaCACGATGGCCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCCGCAAGACCGCCTGCTG‑A‑TGCAGGACGACACTATCATTATTGACGAAGCGCa                                                             >  1:7929880/1‑101 (MQ=255)
                                        acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                            >  1:8586732/1‑101 (MQ=255)
                                        acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                            <  1:1719701/101‑1 (MQ=255)
                                         cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                           <  1:10026515/101‑1 (MQ=255)
                                          acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGa                                                          >  1:1092001/1‑101 (MQ=255)
                                          acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGa                                                          >  1:5889756/1‑101 (MQ=255)
                                           cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                                         >  1:9949630/1‑101 (MQ=255)
                                           cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                                         >  1:6307052/1‑101 (MQ=255)
                                            gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                        <  1:1670083/101‑1 (MQ=255)
                                            ctgcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑C‑GTGAGTACGACACCATCATCATCGATGAGGCGCACGAAc                                                        >  1:965618/11‑101 (MQ=255)
                                             aTGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACg                                                       <  1:322369/101‑1 (MQ=255)
                                              tGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                                      >  1:1099335/1‑101 (MQ=255)
                                              tGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                                      >  1:9718466/1‑101 (MQ=255)
                                              gcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCG‑TG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGc                                                      <  1:5355154/93‑1 (MQ=255)
                                                 tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                                   <  1:1999844/101‑1 (MQ=255)
                                                   aaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑C‑GTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCt                                                 >  1:308764/4‑101 (MQ=255)
                                                   aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                                 >  1:8023807/1‑101 (MQ=255)
                                                   aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                                 <  1:10479829/101‑1 (MQ=255)
                                                     gCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGa                                               <  1:13869394/101‑1 (MQ=255)
                                                     gCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGa                                               >  1:2914896/1‑101 (MQ=255)
                                                      cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                              >  1:2243244/1‑101 (MQ=255)
                                                      cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑C‑GTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                              >  1:13123256/1‑101 (MQ=255)
                                                       tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                             <  1:6131085/101‑1 (MQ=255)
                                                       tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑CGTG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                             <  1:5787573/101‑2 (MQ=255)
                                                       tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑C‑GTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                             >  1:4967340/1‑100 (MQ=255)
                                                        gatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑C‑GTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaaca                                            <  1:6467691/101‑3 (MQ=255)
                                                         atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAatat                                           >  1:6838311/1‑101 (MQ=255)
                                                          tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCG‑TG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                          <  1:6391924/101‑1 (MQ=255)
                                                          tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCG‑TG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                          <  1:6332820/101‑1 (MQ=255)
                                                           gaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCG‑TG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCg                                         <  1:2079879/101‑1 (MQ=255)
                                                           gaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCG‑TG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCg                                         <  1:6760842/101‑1 (MQ=255)
                                                            aCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGa                                        >  1:4144650/1‑101 (MQ=255)
                                                             ccGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑CGTG‑AGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGac                                       <  1:9217276/101‑2 (MQ=255)
                                                               gACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAttt                                     >  1:5287739/1‑101 (MQ=255)
                                                                 cGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAttttt                                   <  1:10695431/101‑1 (MQ=255)
                                                                  ggTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTG‑C‑GTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGActtcct                                  <  1:11865980/101‑4 (MQ=255)
                                                                   gTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTg                                 >  1:3651969/1‑101 (MQ=255)
                                                                     aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                               >  1:8401873/1‑101 (MQ=255)
                                                                     aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                               >  1:1248278/1‑101 (MQ=255)
                                                                     aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                               <  1:2870609/101‑1 (MQ=255)
                                                                     aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                               >  1:6609943/1‑101 (MQ=255)
                                                                      tCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                              <  1:11185267/101‑1 (MQ=255)
                                                                      tCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                              >  1:13747336/1‑101 (MQ=255)
                                                                       ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                             >  1:3006430/1‑101 (MQ=255)
                                                                        ctgctgGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                            <  1:552005/101‑1 (MQ=255)
                                                                        ctgctgGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                            <  1:11440066/101‑1 (MQ=255)
                                                                           ctgGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTa                         <  1:756037/101‑1 (MQ=255)
                                                                               cGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg                     <  1:3840893/101‑1 (MQ=255)
                                                                               cGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg                     <  1:8217490/101‑1 (MQ=255)
                                                                                ggAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGTATATCGATTTTTTGCTCGGCTATTTGa                    <  1:2292028/101‑1 (MQ=255)
                                                                                 gagaTCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaa                   <  1:8146920/101‑1 (MQ=255)
                                                                                 gagaTCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaa                   <  1:14131094/101‑1 (MQ=255)
                                                                                   gaTCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag                 <  1:5659244/101‑1 (MQ=255)
                                                                                   gaTCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag                 <  1:7105528/101‑1 (MQ=255)
                                                                                     tCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag               <  1:6236081/101‑1 (MQ=255)
                                                                                       cAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt             <  1:7456293/101‑1 (MQ=255)
                                                                                        aGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtg            >  1:9539460/1‑101 (MQ=255)
                                                                                          cAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgct          >  1:12642924/1‑101 (MQ=255)
                                                                                          cAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgct          >  1:11339144/1‑101 (MQ=255)
                                                                                           aaGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgctg         >  1:13946378/1‑101 (MQ=255)
                                                                                             gACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGcc       >  1:10796120/1‑101 (MQ=255)
                                                                                              aCCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcg      >  1:7013567/1‑101 (MQ=255)
                                                                                              aCCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcg      <  1:12215007/101‑1 (MQ=255)
                                                                                              aCCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcg      >  1:1613302/1‑101 (MQ=255)
                                                                                               ccGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcgc     >  1:3895108/1‑101 (MQ=255)
                                                                                                 gCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg   <  1:9363025/101‑1 (MQ=255)
                                                                                                  |                                                                                                      
ATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTG‑A‑TGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGC  >  NC_007779/1485234‑1485432

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.