Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,336 C→T 10.9% L169L (CTG→TTG)  hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,3360C→T10.9% 63.7 ‑0.2 ‑0.0248L169L (CTG→TTG) hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (11/16):  ref base (99/121):  total (111/137)
Fisher's exact test strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

Only 100 of 268 total aligned reads displayed.

CGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTC  >  NC_007779/1485236‑1485435
                                                                                                    |                                                                                                    
cGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGATCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgc                                                                                                      >  1:1197248/1‑101 (MQ=255)
 ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                                     >  1:11788407/1‑101 (MQ=255)
 ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                                     >  1:7253840/1‑101 (MQ=255)
  gTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgctg                                                                                                    <  1:13644910/101‑1 (MQ=255)
     aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATg                                                                                                <  1:1764270/101‑1 (MQ=255)
     aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATg                                                                                                <  1:6345501/101‑1 (MQ=255)
      cAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGc                                                                                               >  1:2172154/1‑101 (MQ=255)
       aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                              >  1:12245344/1‑101 (MQ=255)
       aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                              >  1:5051853/1‑101 (MQ=255)
        aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                             <  1:6888212/101‑1 (MQ=255)
        aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                             >  1:8883658/1‑101 (MQ=255)
         aGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                                            >  1:1005031/1‑101 (MQ=255)
         aGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                                            <  1:11226277/101‑1 (MQ=255)
         aGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                                            <  1:5642434/101‑1 (MQ=255)
          gTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTa                                                                                           <  1:5833192/101‑1 (MQ=255)
          gTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTa                                                                                           >  1:12694585/1‑101 (MQ=255)
          gTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTa                                                                                           >  1:12793672/1‑101 (MQ=255)
              gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGac                                                                                       >  1:5522881/1‑101 (MQ=255)
                ttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGacac                                                                                     <  1:14129773/101‑1 (MQ=255)
                  cAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTa                                                                                   >  1:2641439/1‑101 (MQ=255)
                    gCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATc                                                                                 >  1:1603257/1‑101 (MQ=255)
                     cGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCa                                                                                >  1:1657823/1‑101 (MQ=255)
                      gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                               >  1:3676848/1‑101 (MQ=255)
                      gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                               >  1:13441968/1‑101 (MQ=255)
                      gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                               <  1:7642373/101‑1 (MQ=255)
                       aTCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatt                                                                              <  1:9141936/101‑1 (MQ=255)
                       aTCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatt                                                                              >  1:8367022/1‑101 (MQ=255)
                          aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                                           <  1:3466101/101‑1 (MQ=255)
                          aCGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattatt                                                                           <  1:12922075/101‑1 (MQ=255)
                            gTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGa                                                                         >  1:10890369/1‑101 (MQ=255)
                             tAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGAc                                                                        >  1:3944681/1‑101 (MQ=255)
                             tAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGAc                                                                        >  1:790613/1‑101 (MQ=255)
                              aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACg                                                                       >  1:13722217/1‑101 (MQ=255)
                                gTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGaa                                                                     <  1:11282295/101‑1 (MQ=255)
                                      acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                               <  1:9642305/101‑1 (MQ=255)
                                       cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                              >  1:4337186/1‑101 (MQ=255)
                                       cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                              <  1:7227648/101‑1 (MQ=255)
                                       cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                              <  1:7821079/101‑1 (MQ=255)
                                        acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGa                                                             <  1:12341544/101‑1 (MQ=255)
                                          gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                           >  1:8128554/1‑101 (MQ=255)
                                          ctgcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAAc                                                           <  1:1399495/91‑1 (MQ=255)
                                           tgcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACg                                                          >  1:1901017/10‑101 (MQ=255)
                                           tgcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACg                                                          <  1:4978027/92‑1 (MQ=255)
                                           aTGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACg                                                          <  1:8972207/101‑1 (MQ=255)
                                            tGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                                         <  1:9075611/101‑1 (MQ=255)
                                            gcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGc                                                         >  1:5420356/9‑101 (MQ=255)
                                                cAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGcc                                                     >  1:7040777/1‑101 (MQ=255)
                                                    cTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGaa                                                 <  1:2687897/101‑1 (MQ=255)
                                                     tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                                >  1:289535/1‑101 (MQ=255)
                                                     tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                                >  1:5136135/1‑101 (MQ=255)
                                                     tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                                >  1:5735005/1‑101 (MQ=255)
                                                      gatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAata                                               >  1:3858077/1‑101 (MQ=255)
                                                      gatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAata                                               >  1:10618268/1‑101 (MQ=255)
                                                      gatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaaca                                               <  1:12507595/101‑3 (MQ=255)
                                                      gatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaaca                                               <  1:458712/101‑3 (MQ=255)
                                                       atgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAacat                                              <  1:3309664/101‑1 (MQ=255)
                                                        tgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATc                                             <  1:9525527/101‑1 (MQ=255)
                                                         gaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCg                                            <  1:2594173/101‑1 (MQ=255)
                                                          aCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGa                                           <  1:7854827/101‑1 (MQ=255)
                                                            cGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtt                                         >  1:2383482/1‑101 (MQ=255)
                                                            cGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGact                                         <  1:11173074/101‑3 (MQ=255)
                                                               cGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAttttt                                      <  1:9771593/101‑1 (MQ=255)
                                                                ggTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttttt                                     <  1:8530847/101‑1 (MQ=255)
                                                                ggTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGActtcct                                     <  1:8459031/101‑4 (MQ=255)
                                                                 gTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTg                                    <  1:4046325/101‑1 (MQ=255)
                                                                 gTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTg                                    <  1:12625329/101‑1 (MQ=255)
                                                                   aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                                  >  1:3487590/1‑101 (MQ=255)
                                                                    tCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                                 >  1:2000548/1‑101 (MQ=255)
                                                                     ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                                <  1:1142858/101‑1 (MQ=255)
                                                                     ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                                <  1:2805308/101‑1 (MQ=255)
                                                                     ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                                >  1:13381677/1‑101 (MQ=255)
                                                                      ctgctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                               <  1:1062501/101‑1 (MQ=255)
                                                                      ctgctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACCACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                               <  1:14328955/101‑1 (MQ=255)
                                                                         ctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTa                            <  1:2103716/101‑1 (MQ=255)
                                                                         ctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTa                            <  1:6953120/101‑1 (MQ=255)
                                                                          tgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAt                           >  1:5550871/1‑101 (MQ=255)
                                                                           ggCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAtt                          >  1:1072308/1‑101 (MQ=255)
                                                                              ggAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGa                       <  1:6513313/101‑1 (MQ=255)
                                                                                agaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaaa                     <  1:4732761/101‑1 (MQ=255)
                                                                                 gaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag                    >  1:12785349/1‑101 (MQ=255)
                                                                                 gaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag                    <  1:5297316/101‑1 (MQ=255)
                                                                                    ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                 <  1:2405848/101‑1 (MQ=255)
                                                                                     cAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt                <  1:4484096/101‑1 (MQ=255)
                                                                                     cAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt                >  1:11457651/1‑101 (MQ=255)
                                                                                     cAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt                >  1:13859085/1‑101 (MQ=255)
                                                                                      aGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTAACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtg               >  1:6021222/1‑101 (MQ=255)
                                                                                        cAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgct             >  1:7085609/1‑101 (MQ=255)
                                                                                        cAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgct             <  1:10786216/101‑1 (MQ=255)
                                                                                            aCCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcg         <  1:6573505/101‑1 (MQ=255)
                                                                                             ccGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcgc        <  1:8742801/101‑1 (MQ=255)
                                                                                              cGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcgcg       <  1:9930643/101‑1 (MQ=255)
                                                                                               gCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg      >  1:4185144/1‑101 (MQ=255)
                                                                                                ccggctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTTAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggc     <  1:389220/98‑1 (MQ=255)
                                                                                                ccTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggc     <  1:11002397/101‑1 (MQ=255)
                                                                                                 ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg    >  1:2499697/1‑101 (MQ=255)
                                                                                                 ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg    >  1:10359411/1‑101 (MQ=255)
                                                                                                 ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg    <  1:13226619/101‑1 (MQ=255)
                                                                                                  tgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGt   >  1:10156236/1‑101 (MQ=255)
                                                                                                   gctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc  >  1:3210556/1‑101 (MQ=255)
                                                                                                   gctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc  >  1:4858426/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTC  >  NC_007779/1485236‑1485435

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.