Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,400 C→A 13.3% 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→A13.3% 92.9 ‑2.7 ‑0.1254A413E (GCG→GAG) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  new base (21/13):  ref base (97/123):  total (118/136)
Fisher's exact test strand distribution p-value = 6.51e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.10e-01

Only 100 of 277 total aligned reads displayed.

GGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAA  >  NC_007779/13300‑13500
                                                                                                    |                                                                                                     
ggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATcgc                                                                                                       <  1:4850587/101‑1 (MQ=255)
ggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATcgc                                                                                                       <  1:11990100/101‑1 (MQ=255)
  tgtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGa                                                                                                     <  1:10859024/101‑1 (MQ=255)
   gtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaa                                                                                                    >  1:4410381/1‑101 (MQ=255)
   gtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaa                                                                                                    <  1:11253451/101‑1 (MQ=255)
    ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa                                                                                                   <  1:10055631/101‑1 (MQ=255)
    ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa                                                                                                   >  1:11153481/1‑101 (MQ=255)
     tCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaa                                                                                                  <  1:6179805/101‑1 (MQ=255)
       tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                                <  1:9297419/101‑1 (MQ=255)
        gACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAaca                                                                                               <  1:11052334/101‑1 (MQ=255)
          cTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacc                                                                                             >  1:6637810/1‑101 (MQ=255)
            ggTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcaccac                                                                                           <  1:1673565/101‑1 (MQ=255)
             gTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGGGCAACACCACt                                                                                          >  1:4047689/1‑101 (MQ=255)
             gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                          <  1:8011773/101‑1 (MQ=255)
              tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                         >  1:6585505/1‑101 (MQ=255)
              tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                         >  1:3531127/1‑101 (MQ=255)
              tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                         <  1:10919644/101‑1 (MQ=255)
               gaggtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTAt                                                                                        <  1:13045255/92‑1 (MQ=255)
                aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                       >  1:7795509/1‑101 (MQ=255)
                 ggtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATcc                                                                                      <  1:8090745/94‑1 (MQ=255)
                  gtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATccc                                                                                     <  1:12624311/95‑1 (MQ=255)
                  gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                     >  1:1152671/1‑101 (MQ=255)
                   tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                    <  1:4147142/101‑1 (MQ=255)
                   tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                    >  1:1049176/1‑101 (MQ=255)
                    aaaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGa                                                                                   >  1:9031682/1‑101 (MQ=255)
                     aaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGAc                                                                                  >  1:14880736/4‑101 (MQ=255)
                      aggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGAcc                                                                                 >  1:4584375/3‑101 (MQ=255)
                      aaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGAcc                                                                                 <  1:14467244/101‑1 (MQ=255)
                       aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                                                                <  1:5047557/101‑1 (MQ=255)
                       aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                                                                <  1:7740420/101‑1 (MQ=255)
                        gACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCaa                                                                               <  1:3387621/101‑1 (MQ=255)
                         aCGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAg                                                                              >  1:13380766/1‑101 (MQ=255)
                          cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                            <  1:5901343/101‑1 (MQ=255)
                           gtgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcg                                                                           >  1:2721964/4‑100 (MQ=255)
                             gctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCCAGGGCAACACCCCTATCCCGACCCAGcgtt                                                                         >  1:15028004/2‑98 (MQ=255)
                               tgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGc                                                                       <  1:12186114/101‑1 (MQ=255)
                               tgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcg                                                                       >  1:696547/1‑96 (MQ=255)
                                gctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGcc                                                                      >  1:6352528/1‑101 (MQ=255)
                                gctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGcc                                                                      <  1:14441456/101‑1 (MQ=255)
                                 ctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCa                                                                     <  1:7272181/101‑1 (MQ=255)
                                 ctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCa                                                                     >  1:14437417/1‑101 (MQ=255)
                                 ctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcgga                                                                     >  1:6255058/1‑94 (MQ=255)
                                  tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAg                                                                    >  1:14557492/1‑101 (MQ=255)
                                  tgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttctgag                                                                    >  1:8464595/1‑93 (MQ=255)
                                   gctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcggaga                                                                   >  1:9172838/1‑92 (MQ=255)
                                    ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                                                                  <  1:8896873/101‑1 (MQ=255)
                                    ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                                                                  <  1:1330249/101‑1 (MQ=255)
                                      ggACGTTGCCCCGCTGTCTCTGGGTATCGAAACCAGGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtgt                                                                <  1:3945291/101‑1 (MQ=255)
                                       gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                               >  1:9535726/1‑101 (MQ=255)
                                       gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                               >  1:11803643/1‑101 (MQ=255)
                                       gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                               >  1:14762747/1‑101 (MQ=255)
                                        aCGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGGGGGTTGACGCCGCTGACGCGAAAAACCCCAATATCCCCACCCAGCCACACCCAGGTGTtc                                                              >  1:5223523/1‑101 (MQ=255)
                                         cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                             <  1:13941451/101‑1 (MQ=255)
                                         cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                             <  1:8956241/101‑1 (MQ=255)
                                           ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                           <  1:4113451/101‑1 (MQ=255)
                                             aCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAc                                                         <  1:9736830/101‑1 (MQ=255)
                                              ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAcc                                                        >  1:13741995/1‑101 (MQ=255)
                                              ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAcc                                                        <  1:12355812/101‑1 (MQ=255)
                                               cccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCg                                                       >  1:3556963/1‑101 (MQ=255)
                                                 cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                     <  1:2232176/101‑1 (MQ=255)
                                                 cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                     >  1:10225626/1‑101 (MQ=255)
                                                 cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                     <  1:3205873/101‑1 (MQ=255)
                                                  gCTGTCTCTGGGTATCGAAACCTTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTg                                                    <  1:3786242/101‑1 (MQ=255)
                                                    tGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGaa                                                  >  1:5614891/1‑101 (MQ=255)
                                                      tctcTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACGCCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGa                                                <  1:10108929/101‑1 (MQ=255)
                                                       ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                               <  1:5122354/101‑1 (MQ=255)
                                                        tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa                                              >  1:5394877/1‑101 (MQ=255)
                                                        tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa                                              >  1:10173480/1‑101 (MQ=255)
                                                          tGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAAc                                            <  1:7163828/101‑1 (MQ=255)
                                                              tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                        <  1:3437946/101‑1 (MQ=255)
                                                                 cGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTg                                     <  1:5028949/101‑1 (MQ=255)
                                                                  gAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGc                                    <  1:1619110/101‑1 (MQ=255)
                                                                  gAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGc                                    <  1:2390300/101‑1 (MQ=255)
                                                                  gAAACCATGGGCGGTGTGAAGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGc                                    >  1:9283192/1‑101 (MQ=255)
                                                                      ccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTa                                >  1:797230/1‑101 (MQ=255)
                                                                      ccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTa                                <  1:6100003/101‑1 (MQ=255)
                                                                           ggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccat                           <  1:1282302/101‑1 (MQ=255)
                                                                            gCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatc                          <  1:6833397/101‑1 (MQ=255)
                                                                            gCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatc                          >  1:14149282/1‑101 (MQ=255)
                                                                            gCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatc                          >  1:15171932/1‑101 (MQ=255)
                                                                              ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                        <  1:224652/101‑1 (MQ=255)
                                                                                tgtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtg                      >  1:2431582/1‑101 (MQ=255)
                                                                                  tgatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtgtg                    <  1:12437910/101‑1 (MQ=255)
                                                                                   gatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgc                   <  1:8288275/101‑1 (MQ=255)
                                                                                     tgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctg                 >  1:2579362/1‑101 (MQ=255)
                                                                                      gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                <  1:1948599/101‑1 (MQ=255)
                                                                                      gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                <  1:4336352/101‑1 (MQ=255)
                                                                                       acgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCa               <  1:7500614/101‑1 (MQ=255)
                                                                                        cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg              <  1:588362/101‑1 (MQ=255)
                                                                                          acgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg            <  1:7338457/101‑1 (MQ=255)
                                                                                             cTGATCACGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGa         <  1:4019211/101‑1 (MQ=255)
                                                                                              tGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGaa        >  1:3625139/1‑101 (MQ=255)
                                                                                              tGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGaa        >  1:9419460/1‑101 (MQ=255)
                                                                                               gATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc       >  1:3073305/1‑101 (MQ=255)
                                                                                               gATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc       <  1:11273128/101‑1 (MQ=255)
                                                                                                aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg      <  1:8806758/101‑1 (MQ=255)
                                                                                                  cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTa    >  1:8108645/1‑101 (MQ=255)
                                                                                                  cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTa    <  1:10545602/101‑1 (MQ=255)
                                                                                                   gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa   >  1:7579121/1‑101 (MQ=255)
                                                                                                   gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa   >  1:5602902/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAA  >  NC_007779/13300‑13500

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.