Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,323 G→C 15.3% Q164H (CAG→CAC hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,3230G→C15.3% 120.1 ‑3.0 ‑0.0248Q164H (CAG→CAChrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (13/25):  ref base (111/99):  total (124/124)
Fisher's exact test strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.49e-01

Only 100 of 250 total aligned reads displayed.

GGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTG  >  NC_007779/1485224‑1485422
                                                                                                   |                                                                                                     
gggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGc                                                                                                      >  1:3890214/1‑101 (MQ=255)
gggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGc                                                                                                      >  1:8146656/1‑101 (MQ=255)
gggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGc                                                                                                      >  1:2011955/1‑101 (MQ=255)
   cGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAg                                                                                                   >  1:9399921/1‑101 (MQ=255)
    ggTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGa                                                                                                  >  1:13391463/1‑101 (MQ=255)
    ggTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGa                                                                                                  >  1:5541495/1‑101 (MQ=255)
     gTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAc                                                                                                 <  1:11224726/101‑1 (MQ=255)
      ttGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAcc                                                                                                <  1:7993831/101‑1 (MQ=255)
      ttGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAcc                                                                                                >  1:8860901/1‑101 (MQ=255)
          aTCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCct                                                                                            >  1:1689559/1‑101 (MQ=255)
          aTCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGAATCCTGCTGGCGGAGATCCAGCAAGACCGCct                                                                                            >  1:13021008/1‑101 (MQ=255)
           tCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctg                                                                                           >  1:9301524/1‑101 (MQ=255)
             ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                         <  1:9090210/101‑1 (MQ=255)
             ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                         >  1:3243785/1‑101 (MQ=255)
              gTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgctg                                                                                        >  1:4897936/1‑101 (MQ=255)
              gTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgctg                                                                                        >  1:2204757/1‑101 (MQ=255)
               ttACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGa                                                                                      <  1:12365058/101‑1 (MQ=255)
                 aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATg                                                                                   <  1:7421544/101‑1 (MQ=255)
                  cAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGc                                                                                  >  1:5309/1‑101 (MQ=255)
                   aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                 >  1:10236725/1‑101 (MQ=255)
                    aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                <  1:1522024/101‑1 (MQ=255)
                     aGTGCGTTTCAGCGATCCCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                               <  1:357722/101‑1 (MQ=255)
                      gTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTa                                                                              <  1:8222548/101‑1 (MQ=255)
                       tGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTAc                                                                             <  1:5216310/101‑1 (MQ=255)
                       tGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTAc                                                                             >  1:8791045/1‑101 (MQ=255)
                         cGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGa                                                                           >  1:100778/1‑101 (MQ=255)
                          gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGac                                                                          >  1:661919/1‑101 (MQ=255)
                           tttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGaca                                                                         >  1:1401726/1‑101 (MQ=255)
                            ttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGacac                                                                        >  1:3183072/1‑101 (MQ=255)
                             tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACt                                                                       >  1:9846945/1‑101 (MQ=255)
                                 cGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCa                                                                   >  1:4447729/1‑101 (MQ=255)
                                  gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                  >  1:4884349/1‑101 (MQ=255)
                                  gATCACGTAAGTGATAACACGATGGTCAAGCTGACGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                  >  1:395564/1‑101 (MQ=255)
                                    tCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatta                                                                >  1:2575538/1‑101 (MQ=255)
                                     cACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattat                                                               >  1:11886769/1‑101 (MQ=255)
                                      aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                              >  1:14056225/1‑101 (MQ=255)
                                      aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                              <  1:13276779/101‑1 (MQ=255)
                                       cGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTg                                                             <  1:4122125/101‑1 (MQ=255)
                                        gTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGa                                                            <  1:10816801/101‑1 (MQ=255)
                                         tAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGAc                                                           >  1:425465/1‑101 (MQ=255)
                                          aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACg                                                          >  1:5136344/1‑101 (MQ=255)
                                          aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACg                                                          <  1:11798361/101‑1 (MQ=255)
                                           aGTGATAACACGATGGTCAAGCAGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGa                                                         >  1:12329696/1‑101 (MQ=255)
                                             tGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAg                                                       <  1:13562920/101‑1 (MQ=255)
                                              gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGTAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                      >  1:11746237/1‑101 (MQ=255)
                                               aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcg                                                     >  1:2966711/1‑101 (MQ=255)
                                               aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcg                                                     <  1:4046499/101‑1 (MQ=255)
                                               aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcg                                                     <  1:6094406/101‑1 (MQ=255)
                                                 aaCACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCa                                                   >  1:1193898/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:9949630/1‑101 (MQ=255)
                                                     cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                               <  1:7741977/101‑1 (MQ=255)
                                                     cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                               >  1:6307052/1‑101 (MQ=255)
                                                      gATGGTCAAGCTGGTGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                              >  1:3386493/1‑101 (MQ=255)
                                                        tGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                            >  1:1099335/1‑101 (MQ=255)
                                                        tGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                            >  1:9718466/1‑101 (MQ=255)
                                                         cattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCa                                           <  1:1408366/94‑1 (MQ=255)
                                                           ttaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGc                                         >  1:13895876/6‑101 (MQ=255)
                                                           tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                         >  1:7476425/1‑101 (MQ=255)
                                                             aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                       <  1:10479829/101‑1 (MQ=255)
                                                               acTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCtca                                     <  1:12991051/100‑3 (MQ=255)
                                                                cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                    <  1:7060855/101‑1 (MQ=255)
                                                                cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                    >  1:9741732/1‑101 (MQ=255)
                                                                 tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                   <  1:6131085/101‑1 (MQ=255)
                                                                  gatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaaca                                  <  1:6467691/101‑3 (MQ=255)
                                                                   atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAatat                                 <  1:78715/101‑1 (MQ=255)
                                                                   atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAatat                                 >  1:6838311/1‑101 (MQ=255)
                                                                   atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAatat                                 <  1:5841324/101‑1 (MQ=255)
                                                                   atgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAacat                                 <  1:5504919/101‑1 (MQ=255)
                                                                    tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                <  1:6391924/101‑1 (MQ=255)
                                                                    tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                <  1:6332820/101‑1 (MQ=255)
                                                                    tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                <  1:8633877/101‑1 (MQ=255)
                                                                     gaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCg                               <  1:6760842/101‑1 (MQ=255)
                                                                     gaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCg                               <  1:2079879/101‑1 (MQ=255)
                                                                       ccGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGac                             <  1:9217276/101‑2 (MQ=255)
                                                                        cGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtt                            <  1:2291697/101‑1 (MQ=255)
                                                                        cGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtt                            >  1:8556761/1‑101 (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)
                                                                               aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                     <  1:2870609/101‑1 (MQ=255)
                                                                                 ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                   >  1:3006430/1‑101 (MQ=255)
                                                                                  ctgctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                  <  1:11440066/101‑1 (MQ=255)
                                                                                      tgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAt              <  1:8197228/101‑1 (MQ=255)
                                                                                        gCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt            <  1:13182672/101‑1 (MQ=255)
                                                                                         cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           >  1:8044942/1‑101 (MQ=255)
                                                                                         cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           >  1:7439870/1‑101 (MQ=255)
                                                                                         cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           <  1:3088886/101‑1 (MQ=255)
                                                                                         cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg           <  1:3840893/101‑1 (MQ=255)
                                                                                            agaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaaa        <  1:9632573/101‑1 (MQ=255)
                                                                                             gaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag       <  1:5659244/101‑1 (MQ=255)
                                                                                               tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag     <  1:10941179/101‑1 (MQ=255)
                                                                                               tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag     >  1:7685877/1‑101 (MQ=255)
                                                                                               tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag     <  1:6236081/101‑1 (MQ=255)
                                                                                                ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt    >  1:13676248/1‑101 (MQ=255)
                                                                                                  aGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtg  >  1:9539460/1‑101 (MQ=255)
                                                                                                   |                                                                                                     
GGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTG  >  NC_007779/1485224‑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.