Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,690,692 C→T 13.8% noncoding (2286/2904 nt) rrlC ← 23S ribosomal RNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,690,6920C→T13.8% 25.2 ‑1.1 ‑0.087noncoding (2286/2904 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  new base (3/9):  ref base (32/43):  total (35/52)
Fisher's exact test strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.43e-01

TGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAC  >  NC_007779/3690667‑3690792
                         |                                                                                                    
tGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATc                           <  1:10603409/101‑1 (MQ=11)
tGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATc                           <  1:11351805/101‑1 (MQ=11)
 tgccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACTCGGATCa                          >  1:624377/6‑101 (MQ=2)
 gTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCa                          <  1:7788057/101‑1 (MQ=31)
 gTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCACCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCa                          >  1:698724/1‑101 (MQ=14)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:5899676/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:605973/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:551568/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         >  1:10701177/1‑101 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:12263441/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:13197022/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:11761384/101‑1 (MQ=32)
  tCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCAc                         <  1:12692953/101‑1 (MQ=32)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                         <  1:4801809/97‑1 (MQ=2)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                         <  1:11083907/97‑1 (MQ=2)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                         <  1:14284747/97‑1 (MQ=2)
  gccaaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCAc                         <  1:2875789/97‑1 (MQ=2)
   ccGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACGCTGTCCGCAACCCGGATCACg                        >  1:12230546/1‑101 (MQ=17)
   ccGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACg                        <  1:13814883/101‑1 (MQ=32)
   ccGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACg                        <  1:275279/101‑1 (MQ=32)
    caaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACgg                       >  1:4760728/3‑101 (MQ=2)
    cGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACgg                       >  1:11600471/1‑101 (MQ=32)
    cGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACgg                       >  1:7428972/1‑101 (MQ=32)
     gACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACggg                      <  1:623883/101‑1 (MQ=32)
     gACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACggg                      <  1:14184396/101‑1 (MQ=32)
     aaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACggg                      <  1:2132112/100‑1 (MQ=2)
     aaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACggg                      <  1:11392411/100‑1 (MQ=2)
     aaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACggg                      <  1:11282016/100‑1 (MQ=2)
     aaCCAGGCTGAGGGAACCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCAGGCACTGTCCCTAACCCGGATCACggg                      <  1:4993576/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)
       ccAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTc                    >  1:2580015/1‑101 (MQ=32)
        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:24118/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                 <  1:13247630/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                 >  1:9909284/1‑101 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                 <  1:11669485/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                 <  1:11484198/101‑1 (MQ=35)
          ggATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCAc                 <  1:10806971/101‑1 (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:13494670/1‑101 (MQ=35)
            aTTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt               <  1:8916525/101‑1 (MQ=35)
            aTTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt               >  1:11242962/1‑101 (MQ=35)
            aTTAGCCAACCTTCGTGCTCCTCCGTTACTCTTGAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGt               <  1:7643929/101‑1 (MQ=25)
             ttAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGtt              <  1:11483906/101‑1 (MQ=37)
             tgagggaaCCTTTGGGCGCCTCCGTTACTCTTTGGGAGGCAACCGCCCCAGTTAAACTACCCACCACGCACTGTCCCTAACCCGGACCACGGCCCCAGGtt              >  1:3560106/7‑101 (MQ=2)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGGCCGACCCAGTCAAACTACCCACCAGACACTATCTGCAACCCGGATCACGGGACCACGTTa             >  1:13813698/1‑101 (MQ=11)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTa             <  1:3182578/101‑1 (MQ=35)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTa             >  1:2087600/1‑101 (MQ=35)
              tAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTa             >  1:4334454/1‑101 (MQ=35)
                gCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGa           <  1:5451100/101‑1 (MQ=37)
                gCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGa           <  1:10030795/101‑1 (MQ=37)
                gCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGa           >  1:13593/1‑101 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          <  1:10159659/101‑1 (MQ=37)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          >  1:10748849/1‑101 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          >  1:11590181/1‑101 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          <  1:11841460/101‑1 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          <  1:8079234/101‑1 (MQ=37)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          <  1:7874073/101‑1 (MQ=35)
                 ccAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGaa          >  1:6193739/1‑101 (MQ=35)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc         <  1:5190063/101‑1 (MQ=35)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc         <  1:8674723/101‑1 (MQ=37)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc         >  1:7899105/1‑101 (MQ=35)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc         <  1:12812611/101‑1 (MQ=37)
                  cAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAAc         >  1:13301946/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa        >  1:14035885/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa        >  1:4326087/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa        >  1:2875458/1‑101 (MQ=35)
                   aaCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACa        >  1:3652769/1‑101 (MQ=35)
                    aCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACAt       <  1:12776279/101‑1 (MQ=37)
                    aCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACAt       <  1:11974329/101‑1 (MQ=37)
                    aCCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACAt       >  1:2881165/1‑101 (MQ=35)
                     ccTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATc      >  1:13045590/1‑101 (MQ=35)
                     ccTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATc      >  1:12104119/1‑101 (MQ=35)
                     ccTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATc      <  1:902233/101‑1 (MQ=38)
                      cTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCa     <  1:7763592/101‑1 (MQ=37)
                      cTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCa     >  1:4886280/1‑101 (MQ=35)
                       ttCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaa    <  1:13365226/101‑1 (MQ=38)
                       ttCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaa    >  1:7391111/1‑101 (MQ=35)
                        tCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaaa   >  1:9998156/1‑101 (MQ=35)
                        tCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaaa   <  1:5149445/101‑1 (MQ=38)
                        tCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCaaa   <  1:5807047/101‑1 (MQ=35)
                         cGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAc  >  1:1855895/1‑101 (MQ=35)
                         cGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAc  <  1:397650/101‑1 (MQ=35)
                         cGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAc  <  1:6937137/101‑1 (MQ=35)
                         |                                                                                                    
TGTCCGACCAGGATTAGCCAACCTTCGTGCTCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATCACGGGTCCACGTTAGAACATCAAAC  >  NC_007779/3690667‑3690792

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