Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,690,713 A→G 8.4% noncoding (2265/2904 nt) rrlC ← 23S ribosomal RNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,690,7130A→G 8.4% 25.4 ‑1.0 ‑0.2143noncoding (2265/2904 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  new base (3/9):  ref base (54/77):  total (57/86)
Fisher's exact test strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.22e-01

Only 100 of 148 total aligned reads displayed.

TGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAG  >  NC_007779/3690667‑3690813
                                              |                                                                                                    
tGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATc                                                <  1:11351805/101‑1 (MQ=11)
tGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATc                                                <  1:10603409/101‑1 (MQ=11)
 gTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCa                                               <  1:7788057/101‑1 (MQ=31)
 gTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCACCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCa                                               >  1:698724/1‑101 (MQ=14)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                                              <  1:605973/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                                              <  1:12263441/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                                              <  1:551568/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                                              <  1:11761384/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                                              <  1:13197022/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                                              <  1:5899676/101‑1 (MQ=32)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                                              <  1:11083907/97‑1 (MQ=2)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                                              <  1:14284747/97‑1 (MQ=2)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                                              <  1:2875789/97‑1 (MQ=2)
    caaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACgg                                            >  1:4760728/3‑101 (MQ=2)
    cGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACgg                                            >  1:11600471/1‑101 (MQ=32)
     gACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACggg                                           <  1:623883/101‑1 (MQ=32)
     aaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACggg                                           <  1:2132112/100‑1 (MQ=2)
     aaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACggg                                           <  1:9429745/100‑1 (MQ=2)
      aCCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGt                                          <  1:11468777/101‑1 (MQ=32)
       ccAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACggggc                                         >  1:9625322/1‑99 (MQ=17)
        cAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTcc                                        <  1:9055766/101‑1 (MQ=33)
        cAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGACCGCAACCCGGATCACGGGTcc                                        >  1:5672622/1‑101 (MQ=18)
         aGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTccc                                       >  1:6736287/1‑100 (MQ=18)
         aGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCa                                       <  1:1977420/101‑1 (MQ=34)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                                      <  1:13247630/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                                      <  1:11669485/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                                      <  1:10806971/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                                      >  1:9909284/1‑101 (MQ=35)
           gATTAGCCACCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTCCCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACg                                     <  1:11321058/101‑1 (MQ=255)
           gATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACg                                     >  1:4534304/1‑101 (MQ=35)
            aTTAGCCAACCTTCGTGCTCCTCCGTTCCTCTTTAGGAGGAGACCGCCCCAGTCAAACTTCCCCCCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt                                    <  1:2697860/101‑1 (MQ=255)
            aTTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt                                    >  1:11242962/1‑101 (MQ=35)
            aTTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt                                    >  1:13494670/1‑101 (MQ=35)
            aTTAGCCAACCTTCGTGCTCCTCCGTTACTCTTGAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt                                    <  1:7643929/101‑1 (MQ=25)
             ttAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGtt                                   <  1:11483906/101‑1 (MQ=37)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGGCCGACCCAGTCAAACTACCCACCAGACACTATCTGCAACCCGGATCACGGGACCACGTTa                                  >  1:13813698/1‑101 (MQ=11)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTa                                  >  1:4334454/1‑101 (MQ=35)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTa                                  <  1:3182578/101‑1 (MQ=35)
                gCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGa                                >  1:13593/1‑101 (MQ=35)
                gCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGa                                <  1:5451100/101‑1 (MQ=37)
                gCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGa                                <  1:10030795/101‑1 (MQ=37)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa                               >  1:6193739/1‑101 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa                               <  1:7874073/101‑1 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa                               <  1:8079234/101‑1 (MQ=37)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa                               >  1:11590181/1‑101 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa                               >  1:10748849/1‑101 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa                               <  1:10159659/101‑1 (MQ=37)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc                              >  1:13301946/1‑101 (MQ=35)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc                              <  1:12812611/101‑1 (MQ=37)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc                              >  1:7899105/1‑101 (MQ=35)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc                              <  1:8674723/101‑1 (MQ=37)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa                             >  1:4326087/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa                             >  1:3652769/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa                             >  1:2875458/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa                             >  1:14035885/1‑101 (MQ=35)
                    aCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACAt                            >  1:2881165/1‑101 (MQ=35)
                    aCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACAt                            <  1:11974329/101‑1 (MQ=37)
                    aCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACAt                            <  1:12776279/101‑1 (MQ=37)
                     ccTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATc                           >  1:12104119/1‑101 (MQ=35)
                     ccTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATc                           <  1:902233/101‑1 (MQ=38)
                      cTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCa                          >  1:4886280/1‑101 (MQ=35)
                       ttCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaa                         >  1:7391111/1‑101 (MQ=35)
                        tCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaaa                        <  1:5149445/101‑1 (MQ=38)
                         cGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAc                       >  1:1855895/1‑101 (MQ=35)
                         cGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAc                       <  1:397650/101‑1 (MQ=35)
                          gTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACa                      <  1:12138271/101‑1 (MQ=38)
                          gTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACa                      <  1:9172778/101‑1 (MQ=38)
                           tGCTCCTCCGTTACTCTTTAGGAGGAGGCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACAt                     <  1:507380/101‑1 (MQ=37)
                           tGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACAt                     <  1:1313409/101‑1 (MQ=35)
                           tGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACAt                     <  1:6427863/101‑1 (MQ=37)
                           tGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACAt                     <  1:14324521/101‑1 (MQ=35)
                              tcctccGTTTCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaa                  <  1:10315056/101‑1 (MQ=25)
                              tcctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaa                  <  1:4806705/101‑1 (MQ=38)
                              tcctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaa                  <  1:2275095/101‑1 (MQ=37)
                              tcctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaa                  <  1:13708530/101‑1 (MQ=38)
                               cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTGCCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaaa                 <  1:3394639/101‑1 (MQ=25)
                               cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaaa                 >  1:7619067/1‑101 (MQ=35)
                               cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaaa                 >  1:3518477/1‑101 (MQ=35)
                               cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaaa                 >  1:5774311/1‑101 (MQ=35)
                               cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaaa                 >  1:3479038/1‑101 (MQ=35)
                               cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTaaa                 <  1:11926115/101‑1 (MQ=35)
                                 tccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAgg               >  1:2659798/1‑101 (MQ=35)
                                 tccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAgg               >  1:7286301/1‑101 (MQ=35)
                                   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggt             <  1:6946964/101‑1 (MQ=37)
                                   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggt             <  1:2053202/101‑1 (MQ=35)
                                   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggt             >  1:5307556/1‑101 (MQ=35)
                                    gTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtg            >  1:12555810/1‑101 (MQ=35)
                                    gTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtg            <  1:10238713/101‑1 (MQ=38)
                                    gTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtg            <  1:4847172/101‑1 (MQ=38)
                                     ttACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtgg           >  1:12861085/1‑101 (MQ=35)
                                     ttACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtgg           >  1:8819242/1‑101 (MQ=35)
                                     ttACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACAGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtgg           >  1:5581237/1‑101 (MQ=25)
                                      tACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGggtggt          >  1:4940399/1‑101 (MQ=35)
                                       aCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTa         <  1:10021071/101‑1 (MQ=35)
                                       aCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTa         >  1:2440240/1‑101 (MQ=35)
                                         tctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTAtt       >  1:13870298/1‑101 (MQ=35)
                                          ctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTAttt      >  1:1050228/1‑101 (MQ=35)
                                            ttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTATTTCa    <  1:1994992/101‑1 (MQ=35)
                                            ttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTATTTCa    <  1:2573595/101‑1 (MQ=35)
                                             tAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTATTTCaa   >  1:605285/1‑101 (MQ=35)
                                              |                                                                                                    
TGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAG  >  NC_007779/3690667‑3690813

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