Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,400 C→A 7.7% A413E (GCG→GAG)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_00777913,4000C→A 7.7% 47.8 ‑0.4 ‑0.2273A413E (GCG→GAG) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  new base (12/9):  ref base (129/123):  total (141/132)
Fisher's exact test strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.83e-01

Only 100 of 284 total aligned reads displayed.

GGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAA  >  NC_007779/13300‑13500
                                                                                                    |                                                                                                    
ggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATcgc                                                                                                      >  1:5988064/1‑101 (MQ=255)
    ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa                                                                                                  <  1:10095205/101‑1 (MQ=255)
       tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                               >  1:12811591/1‑101 (MQ=255)
       tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                               >  1:9091326/1‑101 (MQ=255)
       tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                               >  1:11729468/1‑101 (MQ=255)
        gACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAaca                                                                                              <  1:4606391/101‑1 (MQ=255)
         aCTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAacac                                                                                             >  1:10156010/1‑101 (MQ=255)
          cTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacc                                                                                            <  1:5940407/101‑1 (MQ=255)
          cTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacc                                                                                            <  1:7608581/101‑1 (MQ=255)
           tGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAacacaa                                                                                           >  1:6055943/1‑99 (MQ=255)
            ggTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAAcaccac                                                                                          <  1:283575/101‑1 (MQ=255)
             gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                         <  1:9194095/101‑1 (MQ=255)
              tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                        >  1:8183539/1‑101 (MQ=255)
              tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                        <  1:11924110/101‑1 (MQ=255)
               gACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                       >  1:6390587/1‑101 (MQ=255)
                aggtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATc                                                                                      >  1:7386556/9‑101 (MQ=255)
                aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACGATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                      >  1:8949511/1‑101 (MQ=255)
                aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                      <  1:13487717/101‑1 (MQ=255)
                aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                      >  1:11097324/1‑101 (MQ=255)
                 cGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAATCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATcc                                                                                     <  1:11325235/101‑1 (MQ=255)
                  gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                    >  1:12619079/1‑101 (MQ=255)
                  gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                    <  1:4865504/101‑1 (MQ=255)
                       aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                                                               >  1:6717284/1‑101 (MQ=255)
                         aCGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAg                                                                             <  1:10612801/101‑1 (MQ=255)
                          cGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGc                                                                            <  1:4140094/101‑1 (MQ=255)
                          cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                            <  1:6207860/101‑1 (MQ=255)
                           gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGca                                                                           >  1:13084181/1‑101 (MQ=255)
                           gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCCGATCGCGAAAAACACCACTATCCCGACCAAGca                                                                           <  1:3414473/101‑1 (MQ=255)
                             aCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcaca                                                                         <  1:10757105/101‑1 (MQ=255)
                             aCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcaca                                                                         <  1:9588983/101‑1 (MQ=255)
                             aCTGCTGCTGGACGTTACCCCGCCGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcaca                                                                         <  1:3343322/101‑1 (MQ=255)
                              ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                        <  1:8334975/101‑1 (MQ=255)
                              ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                        <  1:561047/101‑1 (MQ=255)
                              ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                        <  1:213903/101‑1 (MQ=255)
                                gctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGcc                                                                      <  1:10867515/101‑1 (MQ=255)
                                gctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcgg                                                                      >  1:2661050/1‑95 (MQ=255)
                                  tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCTACCAAGCACAGCCAg                                                                    >  1:9782811/1‑101 (MQ=255)
                                     tgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtg                                                                 <  1:12552224/101‑1 (MQ=255)
                                       gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAACACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                               >  1:1745378/1‑101 (MQ=255)
                                       gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                               >  1:8116005/1‑101 (MQ=255)
                                         cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                             <  1:13571633/101‑1 (MQ=255)
                                          gTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtctc                                                            <  1:3249066/101‑1 (MQ=255)
                                          gTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtctc                                                            >  1:7918375/1‑101 (MQ=255)
                                           ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                           <  1:4448612/101‑1 (MQ=255)
                                             aCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAc                                                         >  1:5892539/1‑101 (MQ=255)
                                               cccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCg                                                       >  1:12426289/1‑101 (MQ=255)
                                                ccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCTCTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGc                                                      <  1:13788114/101‑1 (MQ=255)
                                                 cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                     <  1:11420900/101‑1 (MQ=255)
                                                  gCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTg                                                    >  1:3565784/1‑101 (MQ=255)
                                                   cTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGa                                                   <  1:992135/101‑1 (MQ=255)
                                                    tGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGaa                                                  <  1:11174978/101‑1 (MQ=255)
                                                    tGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGaa                                                  <  1:6926556/101‑1 (MQ=255)
                                                     gTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAg                                                 >  1:5443340/1‑101 (MQ=255)
                                                       ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                               >  1:10008204/1‑101 (MQ=255)
                                                        tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa                                              >  1:3683747/1‑101 (MQ=255)
                                                        tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa                                              <  1:4415832/101‑1 (MQ=255)
                                                              tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                        >  1:8784478/1‑101 (MQ=255)
                                                              tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                        >  1:14222118/1‑101 (MQ=255)
                                                              tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                        >  1:2920849/1‑101 (MQ=255)
                                                               aTCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTc                                       >  1:3482195/1‑101 (MQ=255)
                                                                tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCt                                      >  1:12309838/1‑101 (MQ=255)
                                                                tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCt                                      <  1:19802/101‑1 (MQ=255)
                                                                 cGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTg                                     <  1:6283285/101‑1 (MQ=255)
                                                                   aaaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCg                                   >  1:5795352/1‑101 (MQ=255)
                                                                   aaaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCg                                   >  1:1343182/1‑101 (MQ=255)
                                                                     aCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGt                                 <  1:10285787/101‑1 (MQ=255)
                                                                       cATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTaa                               <  1:510269/101‑1 (MQ=255)
                                                                       cATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTaa                               <  1:11382693/101‑1 (MQ=255)
                                                                        aTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAc                              >  1:7824087/1‑101 (MQ=255)
                                                                        aTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAc                              >  1:7447172/1‑101 (MQ=255)
                                                                        aTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAc                              >  1:1952144/1‑101 (MQ=255)
                                                                          gggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcca                            >  1:2432190/1‑101 (MQ=255)
                                                                          gggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcca                            <  1:3817408/101‑1 (MQ=255)
                                                                             cGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcc                         >  1:4533685/1‑101 (MQ=255)
                                                                              ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                        >  1:8401101/1‑101 (MQ=255)
                                                                               gtgtGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatccat                       <  1:8057657/101‑1 (MQ=255)
                                                                               gtgtGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatccat                       <  1:3009537/101‑1 (MQ=255)
                                                                                    atgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgct                  <  1:6582530/101‑1 (MQ=255)
                                                                                    atgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgct                  <  1:6412789/101‑1 (MQ=255)
                                                                                    atgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgct                  <  1:4013320/101‑1 (MQ=255)
                                                                                     tgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctg                 <  1:7666907/101‑1 (MQ=255)
                                                                                     tgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctg                 >  1:14148804/1‑101 (MQ=255)
                                                                                         gacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGTGGTAACCATCCATGTGCTGCAgg             <  1:13230128/101‑1 (MQ=255)
                                                                                          acgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg            >  1:14268138/1‑101 (MQ=255)
                                                                                          acgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg            <  1:13716706/101‑1 (MQ=255)
                                                                                           cgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGt           >  1:9449770/1‑101 (MQ=255)
                                                                                            gCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg          <  1:12896748/101‑1 (MQ=255)
                                                                                            gCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg          <  1:54771/101‑1 (MQ=255)
                                                                                               gATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc       <  1:1535440/101‑1 (MQ=255)
                                                                                               gATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc       <  1:2953658/101‑1 (MQ=255)
                                                                                                aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg      <  1:11446258/101‑1 (MQ=255)
                                                                                                aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg      >  1:12130238/1‑101 (MQ=255)
                                                                                                aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg      <  1:7686751/101‑1 (MQ=255)
                                                                                                aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg      >  1:13888931/1‑101 (MQ=255)
                                                                                                 tCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGt     >  1:1446014/1‑101 (MQ=255)
                                                                                                 tCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGt     <  1:556515/101‑1 (MQ=255)
                                                                                                    cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAATGTaaa  <  1:8584918/101‑1 (MQ=255)
                                                                                                    cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa  <  1:9024462/101‑1 (MQ=255)
                                                                                                    cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa  <  1:5523807/101‑1 (MQ=255)
                                                                                                    cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa  >  1:10032539/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAA  >  NC_007779/13300‑13500

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.