Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,101 T→G 11.8% A230A (GCA→GCC ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,1010T→G11.8% 114.4 ‑1.3 ‑0.6279A230A (GCA→GCCftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (13/20):  ref base (125/121):  total (138/141)
Fisher's exact test strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.48e-01

Only 100 of 293 total aligned reads displayed.

GACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTAC  >  NC_007779/3326001‑3326201
                                                                                                     |                                                                                                    
gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                                                                      <  1:10076991/101‑1 (MQ=255)
 aCCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATg                                                                                                     >  1:8950572/1‑101 (MQ=255)
 aCCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATg                                                                                                     <  1:6968747/101‑1 (MQ=255)
   cAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCa                                                                                                   >  1:6819310/1‑101 (MQ=255)
   cAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCAACGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCa                                                                                                   >  1:13899819/1‑101 (MQ=255)
    aGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAc                                                                                                  <  1:2717363/101‑1 (MQ=255)
      cgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAccc                                                                                                <  1:7466254/101‑1 (MQ=255)
      cgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAccc                                                                                                <  1:516215/101‑1 (MQ=255)
       gcCACGCTGGCGGCCTCCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCCCCca                                                                                               <  1:4467038/101‑1 (MQ=255)
        ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                                                                              <  1:2384702/101‑1 (MQ=255)
         cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                                                             <  1:9901728/101‑1 (MQ=255)
         cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                                                             <  1:11874903/101‑1 (MQ=255)
          aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                                                                            <  1:1161842/101‑1 (MQ=255)
            gCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcg                                                                                          <  1:11086068/101‑1 (MQ=255)
              tGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgac                                                                                        <  1:2587568/101‑1 (MQ=255)
               ggcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacg                                                                                       <  1:7259775/101‑1 (MQ=255)
               ggcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacg                                                                                       <  1:11358308/101‑1 (MQ=255)
                 cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                                                                     >  1:10941099/1‑101 (MQ=255)
                  ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                                                    <  1:6015184/101‑1 (MQ=255)
                    ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                                                                                  >  1:524815/1‑101 (MQ=255)
                     cgaCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACAtt                                                                                 <  1:6896897/99‑1 (MQ=255)
                     cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                                                                                 >  1:8573872/1‑101 (MQ=255)
                     cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                                                                                 >  1:7351312/1‑101 (MQ=255)
                      tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                                                                                <  1:3570722/101‑1 (MQ=255)
                       aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                                                                               >  1:5697835/1‑101 (MQ=255)
                       aCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTc                                                                               <  1:6857395/101‑1 (MQ=255)
                       aCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTc                                                                               >  1:4272655/1‑101 (MQ=255)
                        cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                                                              <  1:5889734/101‑1 (MQ=255)
                        cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                                                              >  1:7141262/1‑101 (MQ=255)
                        cGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcg                                                                              >  1:4142879/1‑100 (MQ=255)
                         ggCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTcga                                                                             <  1:7656244/101‑3 (MQ=255)
                          gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                                            <  1:10099382/101‑1 (MQ=255)
                          gCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTAc                                                                            <  1:2224804/101‑1 (MQ=255)
                           cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                                                                           >  1:9498850/1‑101 (MQ=255)
                           cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                                                                           <  1:8126731/101‑1 (MQ=255)
                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                                          >  1:8783185/1‑101 (MQ=255)
                            gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                                                          >  1:2756230/1‑101 (MQ=255)
                              cGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTACTACGAAg                                                                        >  1:4678655/1‑101 (MQ=255)
                                atctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTc                                                                      >  1:2975680/7‑101 (MQ=255)
                                aTTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGCACACGGGATGCACCCACACCGACGAACATTTCTACGAGGTc                                                                      >  1:3919495/1‑101 (MQ=255)
                                 tctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTcc                                                                     <  1:12367150/96‑2 (MQ=255)
                                  ttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAg                                                                    >  1:10692157/1‑101 (MQ=255)
                                  ctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTccg                                                                    >  1:10809357/5‑99 (MQ=255)
                                   tcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGa                                                                   <  1:9777376/98‑1 (MQ=255)
                                   tcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGa                                                                   <  1:8485110/98‑1 (MQ=255)
                                   tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                                                   >  1:1501699/1‑101 (MQ=255)
                                   tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                                                   >  1:6494560/1‑101 (MQ=255)
                                   tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                                                   >  1:2009993/1‑101 (MQ=255)
                                    cgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGaa                                                                  >  1:13365062/3‑101 (MQ=255)
                                      tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                                                                <  1:7066761/101‑1 (MQ=255)
                                       cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                                                               <  1:11266775/101‑1 (MQ=255)
                                        gATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAAcccg                                                              >  1:6158798/1‑99 (MQ=255)
                                        gATAAAGATGATGCACGGTCGCGCTTTCTTCGCCTGTTCGAACACCTCACCAACACGGGATGCACCCACACCGAAGAACATTTCTACGAAGTCAGAACCgg                                                              <  1:10518501/101‑1 (MQ=255)
                                           aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                                                           >  1:5311716/1‑101 (MQ=255)
                                           aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                                                           <  1:1218520/101‑1 (MQ=255)
                                              gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                                                        >  1:9644037/1‑101 (MQ=255)
                                              gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                                                        <  1:13215589/101‑1 (MQ=255)
                                               atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt                                                       <  1:474111/101‑1 (MQ=255)
                                                tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                                                      <  1:4922043/101‑1 (MQ=255)
                                                tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                                                      >  1:3294840/1‑101 (MQ=255)
                                                tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                                                      <  1:7905063/101‑1 (MQ=255)
                                                    gCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGa                                                  >  1:10783844/1‑101 (MQ=255)
                                                     cACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGaa                                                 >  1:13577232/1‑101 (MQ=255)
                                                     cACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGaa                                                 <  1:2127664/101‑1 (MQ=255)
                                                      aCGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGACAGTAAAGAAc                                                <  1:9190972/101‑1 (MQ=255)
                                                       cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg                                               <  1:2831453/101‑1 (MQ=255)
                                                       cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg                                               >  1:6379414/1‑101 (MQ=255)
                                                        ggtccctttttCTTCGCCTGTTCGAACATGTCACGAACACGGGATCCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACgg                                              <  1:12115457/93‑1 (MQ=255)
                                                        ggTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACgg                                              >  1:11491632/1‑101 (MQ=255)
                                                           gCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAAc                                           <  1:8013814/101‑1 (MQ=255)
                                                            ccGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACt                                          <  1:10346872/101‑1 (MQ=255)
                                                             cGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACtt                                         <  1:6667799/101‑1 (MQ=255)
                                                             cGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACtt                                         <  1:299567/101‑1 (MQ=255)
                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                                        >  1:2529521/1‑101 (MQ=255)
                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                                        <  1:8433475/101‑1 (MQ=255)
                                                              gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt                                        >  1:7811835/1‑101 (MQ=255)
                                                              ccTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCGCCCACACCGACGAACATTTCTACGAAGGCAGAACCGGAGATAGTAAAGAACGGAACttt                                        <  1:1605874/100‑1 (MQ=255)
                                                               cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc                                       <  1:6237428/101‑1 (MQ=255)
                                                               cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc                                       <  1:1846293/101‑1 (MQ=255)
                                                               cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc                                       <  1:12407921/101‑1 (MQ=255)
                                                                tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg                                      <  1:12765039/101‑1 (MQ=255)
                                                                 ttcttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgc                                     <  1:3769083/101‑1 (MQ=255)
                                                                    cttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgctt                                   <  1:5045788/101‑1 (MQ=255)
                                                                    cttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcactt                                   >  1:4361974/1‑101 (MQ=255)
                                                                     ttcGCCTGTTCGAACATGTCACGAACAGGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgcttc                                  >  1:6575895/1‑101 (MQ=255)
                                                                         ccTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCt                              <  1:3115736/101‑1 (MQ=255)
                                                                           tGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCNGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGc                            >  1:1926932/1‑101 (MQ=255)
                                                                            gTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCa                           >  1:2617998/1‑101 (MQ=255)
                                                                                gAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCg                       >  1:6129375/1‑101 (MQ=255)
                                                                                gAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCg                       >  1:4848241/1‑101 (MQ=255)
                                                                                   cATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCtt                    <  1:5591407/101‑1 (MQ=255)
                                                                                      gTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCg                 >  1:12968956/1‑101 (MQ=255)
                                                                                      gTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCg                 >  1:9121203/1‑101 (MQ=255)
                                                                                         aCGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCca              <  1:13457037/101‑1 (MQ=255)
                                                                                            aaCACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCcagca           >  1:7586339/1‑101 (MQ=255)
                                                                                            aaCACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCcagca           <  1:3634003/101‑1 (MQ=255)
                                                                                               acGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCg        <  1:12318634/101‑1 (MQ=255)
                                                                                                cGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGt       >  1:14227816/1‑101 (MQ=255)
                                                                                                cGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGt       >  1:1374607/1‑101 (MQ=255)
                                                                                                 gggATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGtt      >  1:3490590/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
GACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTAC  >  NC_007779/3326001‑3326201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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.