Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,322 A→G 15.5% Q164R (CAG→CGG)  hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,3220A→G15.5% 105.7 ‑3.0 ‑0.1247Q164R (CAG→CGG) hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (13/25):  ref base (109/99):  total (123/124)
Fisher's exact test strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.21e-01

Only 100 of 253 total aligned reads displayed.

CCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTG  >  NC_007779/1485222‑1485422
                                                                                                    |                                                                                                      
ccGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCa                                                                                                        >  1:8527398/1‑101 (MQ=255)
  gggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGc                                                                                                      >  1:3890214/1‑101 (MQ=255)
  gggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGc                                                                                                      >  1:8146656/1‑101 (MQ=255)
     cGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAg                                                                                                   >  1:9399921/1‑101 (MQ=255)
       gTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAc                                                                                                 <  1:8765506/101‑1 (MQ=255)
        ttGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAcc                                                                                                <  1:7993831/101‑1 (MQ=255)
        ttGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAcc                                                                                                >  1:8860901/1‑101 (MQ=255)
            aTCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGAATCCTGCTGGCGGAGATCCAGCAAGACCGCct                                                                                            >  1:13021008/1‑101 (MQ=255)
               ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                         >  1:3243785/1‑101 (MQ=255)
                gTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgctg                                                                                        >  1:4897936/1‑101 (MQ=255)
                gTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgctg                                                                                        >  1:2204757/1‑101 (MQ=255)
                   aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATg                                                                                   <  1:7421544/101‑1 (MQ=255)
                     aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                 >  1:10236725/1‑101 (MQ=255)
                      aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                <  1:11921126/101‑1 (MQ=255)
                       aGTGCGTTTCAGCGATCCCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                               <  1:357722/101‑1 (MQ=255)
                       aGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                               >  1:5537429/1‑101 (MQ=255)
                         tGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTAc                                                                             <  1:5216310/101‑1 (MQ=255)
                         tGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTAc                                                                             <  1:421506/101‑1 (MQ=255)
                           cGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGa                                                                           >  1:100778/1‑101 (MQ=255)
                            gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGac                                                                          >  1:661919/1‑101 (MQ=255)
                             tttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGaca                                                                         >  1:1401726/1‑101 (MQ=255)
                             cttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGaca                                                                         <  1:11278164/100‑1 (MQ=255)
                              ttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGacac                                                                        >  1:3183072/1‑101 (MQ=255)
                                   cGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCa                                                                   >  1:4447729/1‑101 (MQ=255)
                                   cGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCa                                                                   >  1:1253779/1‑101 (MQ=255)
                                    gATCACGTAAGTGATAACACGATGGTCAAGCTGACGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                  >  1:395564/1‑101 (MQ=255)
                                      tCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatta                                                                >  1:2575538/1‑101 (MQ=255)
                                        aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                              <  1:13276779/101‑1 (MQ=255)
                                         cGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTg                                                             <  1:4122125/101‑1 (MQ=255)
                                          gTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGa                                                            >  1:8938612/1‑101 (MQ=255)
                                          gTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGa                                                            <  1:10816801/101‑1 (MQ=255)
                                             aGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACCACACTATCATTATTGACGa                                                         >  1:14332944/1‑101 (MQ=255)
                                               tGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAg                                                       <  1:13562920/101‑1 (MQ=255)
                                                gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                      <  1:9795165/101‑1 (MQ=255)
                                                gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                      <  1:14026046/101‑1 (MQ=255)
                                                 aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcg                                                     <  1:4046499/101‑1 (MQ=255)
                                                 aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcg                                                     >  1:2966711/1‑101 (MQ=255)
                                                   aaCACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCa                                                   >  1:12351308/1‑101 (MQ=255)
                                                   aaCACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCa                                                   >  1:13880213/1‑101 (MQ=255)
                                                    acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                  <  1:762810/101‑1 (MQ=255)
                                                    acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                  <  1:1719701/101‑1 (MQ=255)
                                                     cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                 <  1:10026515/101‑1 (MQ=255)
                                                      acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGa                                                >  1:5889756/1‑101 (MQ=255)
                                                       cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                               >  1:6307052/1‑101 (MQ=255)
                                                       cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                               <  1:7741977/101‑1 (MQ=255)
                                                       cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                               >  1:9949630/1‑101 (MQ=255)
                                                        gATGGTCAAGCTGGTGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                              >  1:3386493/1‑101 (MQ=255)
                                                        gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                              <  1:1670083/101‑1 (MQ=255)
                                                          tGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                            >  1:1099335/1‑101 (MQ=255)
                                                           cattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCa                                           <  1:1408366/94‑1 (MQ=255)
                                                             tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                         >  1:7476425/1‑101 (MQ=255)
                                                             tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                         <  1:1999844/101‑1 (MQ=255)
                                                               aaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCt                                       <  1:6325777/98‑1 (MQ=255)
                                                               aaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCt                                       >  1:13636630/4‑101 (MQ=255)
                                                               aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                       <  1:10479829/101‑1 (MQ=255)
                                                                 gCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGa                                     >  1:11789124/1‑101 (MQ=255)
                                                                  cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                    <  1:12313827/101‑1 (MQ=255)
                                                                  cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                    <  1:7060855/101‑1 (MQ=255)
                                                                  cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                    >  1:6234903/1‑101 (MQ=255)
                                                                  cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                    >  1:9741732/1‑101 (MQ=255)
                                                                  cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                    >  1:13123256/1‑101 (MQ=255)
                                                                  cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                    >  1:9486380/1‑101 (MQ=255)
                                                                   tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                   <  1:6131085/101‑1 (MQ=255)
                                                                   tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                   >  1:4788245/1‑101 (MQ=255)
                                                                   tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                   <  1:5787573/101‑2 (MQ=255)
                                                                    gatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaaca                                  <  1:6467691/101‑3 (MQ=255)
                                                                     atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAatat                                 <  1:78715/101‑1 (MQ=255)
                                                                     atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAatat                                 >  1:6838311/1‑101 (MQ=255)
                                                                     atgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAacat                                 <  1:5504919/101‑1 (MQ=255)
                                                                      tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                <  1:8633877/101‑1 (MQ=255)
                                                                      tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                <  1:6391924/101‑1 (MQ=255)
                                                                      tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                <  1:4986333/101‑1 (MQ=255)
                                                                       gaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCg                               >  1:11885461/1‑101 (MQ=255)
                                                                       gaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCg                               <  1:2079879/101‑1 (MQ=255)
                                                                       gaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCg                               <  1:6760842/101‑1 (MQ=255)
                                                                         ccGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGac                             <  1:9217276/101‑2 (MQ=255)
                                                                         ccGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGac                             <  1:5294527/101‑2 (MQ=255)
                                                                          cGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtt                            <  1:2291697/101‑1 (MQ=255)
                                                                             cGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAttttt                         <  1:10695431/101‑1 (MQ=255)
                                                                              ggTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttttt                        >  1:6367562/1‑101 (MQ=255)
                                                                              ggTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttttt                        >  1:1589900/1‑101 (MQ=255)
                                                                              ggTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGActtcct                        <  1:13785570/101‑4 (MQ=255)
                                                                              ggTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGActtcct                        <  1:9300453/101‑4 (MQ=255)
                                                                                 aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                     <  1:2870609/101‑1 (MQ=255)
                                                                                 aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                     >  1:8401873/1‑101 (MQ=255)
                                                                                 aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                     >  1:13346858/1‑101 (MQ=255)
                                                                                  tCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                    <  1:11185267/101‑1 (MQ=255)
                                                                                   ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                   >  1:3006430/1‑101 (MQ=255)
                                                                                        tgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAt              <  1:8197228/101‑1 (MQ=255)
                                                                                          gCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt            >  1:11622384/1‑101 (MQ=255)
                                                                                           cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           <  1:3840893/101‑1 (MQ=255)
                                                                                           cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           <  1:3088886/101‑1 (MQ=255)
                                                                                           cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           >  1:8044942/1‑101 (MQ=255)
                                                                                           cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           >  1:7439870/1‑101 (MQ=255)
                                                                                           cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           <  1:6048684/101‑1 (MQ=255)
                                                                                              agaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaaa        <  1:9632573/101‑1 (MQ=255)
                                                                                               gaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag       <  1:5659244/101‑1 (MQ=255)
                                                                                                aTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAaga      >  1:9689324/1‑101 (MQ=255)
                                                                                                 tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag     <  1:10941179/101‑1 (MQ=255)
                                                                                                 tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag     >  1:7685877/1‑101 (MQ=255)
                                                                                                    |                                                                                                      
CCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTG  >  NC_007779/1485222‑1485422

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