Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,340 T→G 8.0% M170R (ATG→AGG)  hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,3400T→G 8.0% 61.0 ‑1.0 ‑0.2261M170R (ATG→AGG) hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (13/8):  ref base (123/117):  total (136/125)
Fisher's exact test strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.64e-01

Only 100 of 266 total aligned reads displayed.

TACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGAC  >  NC_007779/1485240‑1485440
                                                                                                      |                                                                                                    
 aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATg                                                                                                     <  1:7421544/101‑1 (MQ=255)
  cAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGc                                                                                                    <  1:7750607/101‑1 (MQ=255)
  cAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCGGATGc                                                                                                    >  1:2983928/1‑101 (MQ=255)
   aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                                   >  1:7642363/1‑101 (MQ=255)
   aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                                   >  1:10236725/1‑101 (MQ=255)
    aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                                  <  1:1522024/101‑1 (MQ=255)
    aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                                  <  1:3870689/101‑1 (MQ=255)
    aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                                  <  1:5197667/101‑1 (MQ=255)
      gTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTa                                                                                                <  1:8222548/101‑1 (MQ=255)
         cGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGa                                                                                             >  1:100778/1‑101 (MQ=255)
          gtgtCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGac                                                                                            <  1:11955556/98‑1 (MQ=255)
          gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGac                                                                                            <  1:3487222/101‑1 (MQ=255)
           tttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGaca                                                                                           >  1:1401726/1‑101 (MQ=255)
           cttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGaca                                                                                           <  1:11278164/100‑1 (MQ=255)
             tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACt                                                                                         >  1:13499961/1‑101 (MQ=255)
             tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACt                                                                                         >  1:9846945/1‑101 (MQ=255)
                 cGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCa                                                                                     <  1:5878446/101‑1 (MQ=255)
                 cGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCa                                                                                     >  1:1253779/1‑101 (MQ=255)
                  gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                                    >  1:4884349/1‑101 (MQ=255)
                  gATCACGTAAGTGATAACACGATGGTCAAGCTGACGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                                    >  1:395564/1‑101 (MQ=255)
                   aTCACGTAAGTGATAACACGATGGTCAAGCTGATGACTGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatt                                                                                   >  1:3027883/1‑101 (MQ=255)
                    tCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatta                                                                                  >  1:2575538/1‑101 (MQ=255)
                     cACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattat                                                                                 >  1:11886769/1‑101 (MQ=255)
                      aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                                                <  1:13276779/101‑1 (MQ=255)
                      aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                                                <  1:9792955/101‑1 (MQ=255)
                        gTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGa                                                                              >  1:8938612/1‑101 (MQ=255)
                         tAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGAc                                                                             >  1:425465/1‑101 (MQ=255)
                          aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACATTATCATTATTGACg                                                                            <  1:10570558/101‑1 (MQ=255)
                          aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACg                                                                            >  1:7726089/1‑101 (MQ=255)
                          aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACg                                                                            <  1:11798361/101‑1 (MQ=255)
                           aGTGATAACACGATGGTCAAGCAGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGa                                                                           >  1:12329696/1‑101 (MQ=255)
                              gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGTAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                                        >  1:11746237/1‑101 (MQ=255)
                              gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                                        <  1:14026046/101‑1 (MQ=255)
                               aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcg                                                                       <  1:4046499/101‑1 (MQ=255)
                                 aaCACGATGGCCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCCGCAAGACCGCCTGCTGATGCAGGACGACACTATCATTATTGACGAAGCGCa                                                                     >  1:7929880/1‑101 (MQ=255)
                                  acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                                    >  1:8129983/1‑101 (MQ=255)
                                   cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                                   >  1:681239/1‑101 (MQ=255)
                                   cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                                   <  1:10026515/101‑1 (MQ=255)
                                     cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                                                 >  1:6307052/1‑101 (MQ=255)
                                     cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGaa                                                                 >  1:9949630/1‑101 (MQ=255)
                                      gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                <  1:9423439/101‑1 (MQ=255)
                                      gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                <  1:1670083/101‑1 (MQ=255)
                                       aTGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACg                                                               <  1:322369/101‑1 (MQ=255)
                                        gcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGc                                                              <  1:5355154/93‑1 (MQ=255)
                                         ggTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCa                                                             <  1:8483950/101‑1 (MQ=255)
                                         cattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCa                                                             <  1:1408366/94‑1 (MQ=255)
                                           ttaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGc                                                           >  1:13895876/6‑101 (MQ=255)
                                           tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                                           >  1:7476425/1‑101 (MQ=255)
                                           tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                                           <  1:1999844/101‑1 (MQ=255)
                                             aaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCt                                                         <  1:6325777/98‑1 (MQ=255)
                                             aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                                         <  1:172197/101‑1 (MQ=255)
                                             aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                                         >  1:8023807/1‑101 (MQ=255)
                                             aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                                         <  1:10479829/101‑1 (MQ=255)
                                              aGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTg                                                        <  1:10873455/101‑1 (MQ=255)
                                                cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                                      >  1:2243244/1‑101 (MQ=255)
                                                cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                                      >  1:6234903/1‑101 (MQ=255)
                                                cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                                      >  1:13123256/1‑101 (MQ=255)
                                                cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                                      >  1:7050752/1‑101 (MQ=255)
                                                cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                                      >  1:9741732/1‑101 (MQ=255)
                                                 tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                                     >  1:4788245/1‑101 (MQ=255)
                                                 tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                                     <  1:5787573/101‑2 (MQ=255)
                                                 tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                                     >  1:4967340/1‑100 (MQ=255)
                                                   atgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAacat                                                   <  1:8841200/101‑1 (MQ=255)
                                                    tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                                  <  1:8633877/101‑1 (MQ=255)
                                                    tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                                  <  1:6391924/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)
                                                            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)
                                                               aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                                       >  1:6609943/1‑101 (MQ=255)
                                                                tCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                                      <  1:11185267/101‑1 (MQ=255)
                                                                  ctgctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                                    <  1:11440066/101‑1 (MQ=255)
                                                                  ctgctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                                    <  1:552005/101‑1 (MQ=255)
                                                                     ctgGCGGACATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTAAATATCGATTTTTTGCTCGGCTa                                 <  1:6121179/101‑1 (MQ=255)
                                                                         cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg                             <  1:3840893/101‑1 (MQ=255)
                                                                         cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg                             <  1:6048684/101‑1 (MQ=255)
                                                                           gagaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaa                           >  1:933012/1‑101 (MQ=255)
                                                                            agaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaaa                          <  1:9632573/101‑1 (MQ=255)
                                                                             gaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag                         <  1:7105528/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)
                                                                                ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                      >  1:13676248/1‑101 (MQ=255)
                                                                                 cAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt                     >  1:3110569/1‑101 (MQ=255)
                                                                                 cAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt                     <  1:7456293/101‑1 (MQ=255)
                                                                                  aGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtg                    >  1:9539460/1‑101 (MQ=255)
                                                                                     aaGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgctg                 >  1:4132329/1‑101 (MQ=255)
                                                                                      aGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgctgc                <  1:12698040/101‑1 (MQ=255)
                                                                                       gACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGcc               >  1:10796120/1‑101 (MQ=255)
                                                                                        aCCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcg              <  1:12215007/101‑1 (MQ=255)
                                                                                        aCCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcg              >  1:13598382/1‑101 (MQ=255)
                                                                                          cGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcgcg            <  1:6438843/101‑1 (MQ=255)
                                                                                          cGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcgcg            >  1:6845397/1‑101 (MQ=255)
                                                                                           gCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg           >  1:2311683/1‑101 (MQ=255)
                                                                                             ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg         <  1:5047408/101‑1 (MQ=255)
                                                                                                  gATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTg    >  1:8183926/1‑101 (MQ=255)
                                                                                                      tGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGAc  >  1:3610486/1‑101 (MQ=255)
                                                                                                      |                                                                                                    
TACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGAC  >  NC_007779/1485240‑1485440

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.